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The Temperature Effect of Extensive Air Showers 


By A. L. HODSON 
The Physical Laboratories, The University, Manchester 


Communicated by P. M. S. Blackett; MS. received 3rd Fuly 1951 


ABSTRACT. Extensive air showers have been recorded at sea level over a period of 
nine months. The data show a fairly good correlation between the shower rate and the 
air temperature at about one cascade unit above the apparatus, e.g. at the 950 mb. level. A 
temperature coefficient of —(0°38+0-11)% per °c. at this level is indicated. This may be 
explained by the dependence of the lateral spread of an extensive shower on the air density. 


§1. INTRODUCTION 

N a homogeneous atmosphere the root mean square lateral spread of a 
] cascade shower is proportional to the length of the cascade unit (Blatt 1949) 

which in turn is inversely proportional to the air density p. The dependence 
of the lateral spread on p, and thus on the air temperature, leads us to expect a 
seasonal variation of the rate of extensive showers at sea level. 

The shower rate observed with a fixed arrangement of counter trays depends 
on p in three different ways: 

(a) The effective area of the counter trays is proportional to p*. Since 
varying the areas of the counter trays by a factor f produces a change of f” in 
the shower rate (where y is the exponent of the integral density spectrum of 
extensive showers), a dependence of the order of p?” is thus introduced into the 
shower rate. 

(0) The effective separation of the counter trays is proportional top. Now the 
shower rate recorded by counter trays a distance 7 apart varies approximately 
as7—“ where « lies between 0-1 and 0-5, depending on the experimental arrangement 
and on the altitude (Auger et al. 1939, Daudin and Daudin 1949, Kraybill 1949). 
The shower rate is therefore modified by a further factor p-”. 

(c) A counter tray is liable to be discharged by showers whose cores fall 
within a distance proportional to the root mean square lateral spread, i.e. within 
an area proportional to p~’. 

The variation of the lateral spread with air density is therefore expected to 
produce a total dependence of the shower rate on air density of the 
ovin po”. 

Now Janossy (1948, p. 323) has shown that most of the lateral spread of a 
shower observed at a given atmospheric depth is due to scattering in the previous 
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one or two cascade units. The lateral spread at a given depth in the earth’s — 
atmosphere is therefore more correctly given by taking p equal to the density at — 
one cascade unit above the apparatus. 

At constant pressure, the shower rate is thus expected to depend on the 
absolute temperature T at one cascade unit above the apparatus roughly as 
T~@-2—®, i.e, to have a temperature coefficient of the order of 


1 /dN 100 : 
xan) See 
ere ; 
" eee RPE 9 (1) 


The existence of a temperature effect in extensive air showers was suggested 
by Euler (1940) and by Daudin and Daudin (1949), but no observational evidence 
has previously been reported. 


§2. EXPERIMENTAL RESULTS 

Threefold coincidences were recorded between counter trays A, B, C 
(Figure 1) during the period April 1950-January 1951, the apparatus being 
situated in a heated laboratory under a light roof. The average hourly shower 
rate (corrected to 76 cm. Hg) for each monthly period is given in the Table and. 
plotted in Figure 2. Also shown in Figure 2 is the average air temperature at the 
950 mb. level measured at Liverpool during the same period. (We are indebted 
to the Air Ministry for the temperature data.) 
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Figure 1. 
Records ery: Average rate Average air 
Period dice (ROreS) Eka e per hour in temperature. 
1950 period. at 950 mb. (°r.) 

April-May 688 8038 11-68+ 0-13 43-8 
May-June AT 8196 ikea sede? 52-6 
June-July 688 7714 11-202 0-13 54:1 
July—Aug. 556 6241 11-22+0-14 DSi, 
Aug.—Sept. 530 6449 11:12+0-14 51:9 
Sept.—Oct. 606 6735 11-:11+0-14 47-9 
Oct.—Nov. 538 6265 11:64+0-15 40-3 
Noy.—Dec. 300 3448 11:-49+ 0-20 SS 
Dec.—Jan. 751. B52 4073 Wile 7/ap Oi ke! 31:9 


The results indicate a fairly definite correlation of the shower rate with air 
temperature. An analysis of the data gives the temperature coefficient of extensive 
showers at sea level as —(0:21+0-06)% per °F. or —(0-38+0-11)% per °c. 
(The errors quoted here and elsewhere are standard deviations.) 

The choice of pressure level for the correlation is based on the theoretical 
considerations above. If these are correct the shower data should give the 
highest correlation for a pressure level of about 950 mb. Unfortunately, the 
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and because of the limited statistical accuracy of the shower data it is not possible 

to distinguish between various levels. For example, the seasonal variation of the 
air temperature at the 500 mb. level differed from that at the 950 mb. level by 
less than 10%. ‘The correlation and regression coefficients for these two levels 
are thus not significantly different. 
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Figure 2. The seasonal variation of extensive air showers. 


fe From experiments to be published shortly, the value of y for the showers 
recorded by the arrangement of Figure 1 was found to be y=1-445+0-014. The 
value of « applicable to the present apparatus is probably of the same order as 
the value found by Kraybill (1949) at 720 m. with a threefold arrangement of 
counters («~~ 0-1). From expression (1) we thus expect a temperature coefficient 
of the order of —0-3°% per °c. This is quite consistent with the observed 
value. 


ry 


§3. OTHER VARIATIONS 

Since data were available for the number of showers recorded during each 
hourly period, a search was also made for a possible solar daily variation and for 
a sidereal variation. 

No significant solar daily variation was found, a Fourier analysis of the data 
giving the amplitude of the first harmonic as (0-30 + 0-61)%. ‘The amplitude of 
the average daily variation of air temperature during the period of the experiment 

was 1-2°r. (Compare this with the amplitude of the seasonal temperature 
variation, Figure 2 (of the order of 11°F.).) Assuming the temperature coefficient 
found above a solar diurnal variation of amplitude 0-25 % is expected. 
69-2 
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The analysis of the sidereal data gave the amplitude of the first harmonic o! 
the sidereal variation as (1:15+0-61)% with the ‘maximum’ at 2330 hours 
sidereal time. The average primary energy of the showers recorded was probably’ 
of the order of 5 x 10! ev. It is clear that even for primary cosmic rays of thisi 
energy no large sidereal variation exists. Since the average energy of the showers 
recorded by an arrangement of counters increases with the total counter sprea 
a much larger effect might be observed with counter trays, say, 50 m. apart. 
Daudin and Daudin (1949) claim to have observed a sidereal variation of 1-5°,, 
(for counter spreads 5—70 m.) with a maximum at approximately 1800 hours: 
sidereal time. However, their results were obtained over a rather short period 
and are subject to large statistical errors. Considerably more data are needed: 
to establish whether or not a sidereal variation exists in the very high energy: 
primaries. 
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The Disintegration of Light Nuclei by Meson Capture 


By A. C. CLARK anp 8S. N. RUDDLESDEN 


Department of Theoretical Physics, Liverpool University 
Communicated by H. Fréhlich; MS. received 19th Fune 1951 


ABSTRACT. The disintegration by 7~-meson capture of the elements helium and 
carbon is discussed. ‘T'wo types of meson theory, scalar and pseudo-scalar, are considered. 
The relative probabilities of the various modes of disintegration of helium are calculated 
and the energy spectra of the most likely one, viz. proton emission, are shown. 'The 
disintegration of carbon into two a-particles and one singly charged particle is considered, 
using a simple a-particle model. The energy spectra and angular correlations of the 
a-particles obtained by a direct transition and by a transition via an excited state of ®Be 
are compared. 


§1. INTRODUCTION 

HE production of ‘stars’ by the capture of negative 7-mesons in nuclei 

| is a well-known phenomenon which it is hoped may provide some 
information about the form of the interaction between mesons and 
nucleons. The stars formed by the disintegration of heavy nuclei do not provide’ 
such information since these have been explained satisfactorily by Le Couteur | 
(1950) using the evaporation theory of Bohr and Kalckar (1937), which is based | 
on rather general considerations and is independent of the meson—nucleon | 


interaction. The criteria for the applicability of evaporation theory, that 
: | 
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» (a) the excitation energy given to the nucleus must be much less than its binding 
{ energy, and (5) the number of particles in the nucleus must be great enough for 


approximate thermodynamic equilibrium to be attained between emission of 


q particles, are not obeyed when a light nucleus captures a meson. In this case 
| the concept of a compound nucleus has only a limited meaning and the 


disintegration must be considered as a direct process, as indeed must the 
emission of one high-energy particle from a heavy nucleus (Tamor 1950). This 
direct transition requires the use of first order perturbation theory and the 
results of applying this theory may be expected to show some dependence on 


} the meson-nucleon interaction. Considerations favouring the study of very 


light nuclei are the reduction in the effects due to statistical factors and the 
simplification in the analysis of experiments upon decreasing the number of 
particles concerned. Helium is the lightest element to emit a charged particle 
on capturing a negative meson, so it was chosen as the subject of our 
investigations; also the He nucleus («-particle) is regarded as a constituent of 
heavier nuclei, which enables us to apply our results to the disintegration of 


carbon (or oxygen) into two (or three) «-particles and one singly charged 


particle. 
The direct capture by nuclei of mesons from a beam has been studied by 


Massey and Corben (1939), Sakata and Tanikawa (1939) and Bruno (1948), 
who found capture times of order 10~* second for mesons of a few volts energy; 


_ their results may be taken as applying to 7-mesons since a strong interaction 


with nucleons was assumed. Fermi and Teller (1947) estimated the time taken 
by a meson of energy 2 kev. to be slowed down and captured into the lowest 


Bohr orbit about a nucleus to be about 10- second, in solids. Thus one can 
, 


assume that meson capture takes place from a Bohr orbit, usually the lowest one; 
this is in agreement with experiments, which show that a meson usually comes to 
rest before creating a star. 

Possible basic capture processes are 


Pigeon, 9 7 Gas (1) 
Piao N Tes me reictere (2) 
Rehr Ny, ere = eas (3) 


where P, N, 7, 7° and y represent a proton, neutron, negative 7-meson, neutral 
7-meson and y-ray respectively. Recent experiments at Berkeley (Aamodt et al. 
1950) on the capture of 7-mesons in hydrogen have shown that in 1H, where 
process (1) is forbidden by conservation of momentum, both (2) and (3) occur, 
while in 7H (2) is not observed and (1) and (3) occur with about equal probabilities. 
In heavier nuclei, where the binding energy is too great to permit (2), (3) is 
expected to be rare and process (1) is here assumed. ‘The validity of this 
assumption can be checked experimentally by considering conservation of 
energy and momentum for a disintegration in which only one neutron is emitted : 
in particular, for the disintegration ‘He + 7 —*°H +n the triton has a fixed energy, 
of about 30 Mev., if no secondary particle (or quantum) is emitted. Recent work 
by Menon, Muirhead and Rochat (1950, referred to as M.M.R.) offers some, 
though inconclusive, evidence on this point from stars which are identified as 
eC +7 —>2 4He+*H +n. 

The discovery by Bjorkland e¢ al. (1950) of a neutral 7-meson which decays 
into two y-rays and the work of Yang (1950 a) showing that such a meson must 
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have spin 0 lead one to expect that the spin of the charged 7-meson is also zero, J 
as does much other indirect evidence. Mesons of spin 0 occur in both the scalar | 
and pseudoscalar theories: the interactions used in this paper, between which | 
it is hoped that experiments might differentiate, are the non-relativistic 
approximations to those given by these two theories, viz. $(r) and o.Vi/(r), 
where ¢(r) is the meson wave function at the position of the nucleon and o is the | 
nucleon spin operator. If instead of estimating actual cross sections one estimates | 
the relative probabilities of one nucleus splitting in various ways after capturing 
a meson and the energy spectra of the emitted particles, many uncertainties due 
to our lack of detailed knowledge of the way in which the meson approaches the 
nucleus, the strength of the interaction and the form of the nuclear wave functions 
are eliminated. For each disintegration we refer to a ‘standard’ energy spectrum 
(defined in §2) which is independent of the meson interaction. This 
spectrum is effectively that due to a constant interaction, but neglects the 
antisymmetry of the wave functions; it is convenient both as a rough 
approximation and for comparing the eflects of different interactions. | 
Another way of differentiating between the two theories, proposed by Yang 
(1950 b), depends on measurements of the relative polarization of the y-rays 
produced when the neutral 7-meson disintegrates. 


§2, GENERAL METHOD 
In this section the method of calculation is described in general terms which 
are later applied to specific problems; where necessary *He is used as a simple 
but typical example. Since it is the relative probabilities of the different modes 
of disintegration of one nucleus which are to be calculated, certain constant 
factors, common to all such modes, are omitted. First order perturbation theory 
gives the transition probability for a particular mode of disintegration as 


[2|HPdpp=[ Wap, writing W=Z|HP,  ...... (4) 


where H is the matrix element for transition to one of the possible final states, 
dp, is the number of momentum states per unit energy of the final system, the 
summation is over all spin states and the integration is over all momentum states. 
The energy spectrum P(e,) of one of the emitted particles is given by 


P(e,) de, x «22 de, | Vip a eee (5) 


where «1? de, dpy’ =dpy. The total kinetic energy E available in the final state 
is given by E=meson mass — difference of binding energies in initial and final 
states; the value of the meson mass is taken as 290 electron masses. 

The matrix element for a disintegration is 


H= | es Se ee, (6) 


where ‘Y’; and ’,, are the wave functions describing the initial and final states 
of the nucleons and G is the meson interaction operator, which may be written as 
Ga GT sr (7) = 

nucleons 


where I’ is the operator which converts a proton into a neutron, and is zero when 
applied to a neutron, and for a scalar meson 


Gy =r), aco agile 


while for a pseudo=scalat meson \G3;=0;-Vj/{x,)0 > | (9) 
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where (r,) is the meson wave function at the position of the 7th nucleon and 
_ ois the nucleon spin operator. Since the ‘correct’ forms of ’,; and Y’, are not 
| known it is convenient to represent all nuclei by gaussian wave functions and all 
_ emitted particles by plane waves. All the wave functions are antisymmetrized 
_ with the help of the quantity «,;..,, where 
+1 if2...sis an even permutation of 1,2...” 
€y...g= 7 —lify...sis an odd permutation of 1,2...n 


0. otherwise. 


_ For example, in the complete disintegration of *He the initial wave function is 


qo exp { — abr? 2 Os BeME Neti «ve vee (10) 
ijk 
and a typical final wave function is 
=(3 ie - 2 Hegre sB ans eae eh at (11) 
where + by, = V exp ene 2 ee 


a, B and €, 7 are the spin and isotopic spin functions, 7,; is the distance between 
the ith and jth nucleons, V is the (large) volume in which the events take place, 
and vk,, ik, and #k, are the momenta of the emitted particles. Following 
Bruno (1948) the value of the constant a used throughout is given by 
a2 —4-0 x 10-8 cm. 

The question of the orthogonality of the wave functions W’; and Vy 
will now be discussed briefly. It will be noted that in the above example 
Y, and ‘’; are orthogonal by virtue of the change in isotopic spin but, in the 
approximation of charge-independent forces, one might expect their spatial parts, 
$= exp {—a Xir,?} and py=Y4 93, to be orogens these may be made so, 
artificially, by introducing ,’ =(1 —dexp {—g dir, “j})%p, With suitable values 
of dand g. This factor will reduce, in the region of small 7, Vip the contribution 
to any matrix element between #, and #%,, as would an increase ink, in this region 
due to the action of nuclear forces. If the spatial parts of the initial and final 
wave functions are orthogonal any constant term in the meson interaction G 
will give no contribution to the matrix element. Using the scalar meson 
interaction Gj it has been found that the matrix element | ‘I’,*G,‘l’; after omitting 
the contribution due to the constant term in G, can be identified with the matrix 
element {‘’,’*G,‘¥; in the limiting case g>a. For d=0-5 we find a/g ~ 0-07 
and a calculation has shown that the proton energy spectra from these two matrix 
elements differ, at any point of non-vanishing probability, by less than about 2° 
and that the difference between the total probabilities of proton emission is also 
about 2%. ‘This provides some justification for the use of non-orthogonal wave 
functions with omission of the constant term. In the case of *He such a constant 
term is automatically eliminated by the SUS PeTy properties of the wave 
functions. 

To obtain W=X|H|? where H=J VeGY, one may first integrate over the 
spin and isotopic spin variables and sum over all final states; this gives W in 


+ The nuclear wave functions ¢, %, and 4, should not be confused with ¢(r), the meson 


wave function. 


iprrs 
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terms of integrals J; involving only the spatial wave functions #, and #, and 
that part O, of the interaction operator G, which is independent of spin. The 
nature of the wave functions is such that each integral J; can be expressed in 


Be ae I * a £0; xp (= pit? — 10,8) de, 12 

= | be Obs = | bp ime iano tetees (12) 
where p; and A, are independent of the operator O and depend on a and k, 
respectively. In this integral, J, the coordinates r are referred to the centre of 
mass of the system and integration over the centre-of-mass coordinate is omitted 
since this gives only a constant factor, after ensuring conservation of momentum. 


Thus | 
W =| byp*,|? x an expression which depends on the interactionG. ...... (13) 


This separation makes it possible to compare the energy spectra obtained from 
different interactions not only with each other but with a standard spectrum 


aoe Woe |Site ihe) he ye ee eee (14) 


To obtain these expressions W one must first determine the meson wave 
function %(r) which occurs in G. The meson is considered to be in the lowest 
Bohr orbit about the nucleus before capture; the wave function for this orbit, 
exp(—yr), which was used in our earlier calculations (Ruddlesden and Clark 
1949) should be modified to take into account the spread of charge over the 
nucleus. ‘The wave function now used behaves like exp(—yr) at large r but 
inside the nucleus it may be written as 


i(r)=A—72+ higher terms; ~ 9) | seas. (15) 


the higher terms make little contribution to the result and are neglected. 


§3. DISINTEGRATION OF °“He-AND He 
The individual cases will now be dealt with in more detail, starting with *He 
which can disintegrate in only two ways, viz. 


SHet+a—>1H+2n and ?He+7—>?H-+n. 


Applying the method and wave functions of the previous section to the first 
of these gives, for final states in which the proton has momentum 7k, and the 
neutrons have momenta #k, and hk,, for interaction G,: W, =J,?+1,° + (1, —1;)* 
where J, is defined as 


L=|p*Oyyr where O,=H(e) (16) 
and, for interaction G,: W,=J,?+J,?+(J.—J;)*, where 

J.=[de* 0.9, with 0,50 ir). 9 (17) 
On integration the standard expression W, becomes 


oy lla — 
Wy= | | bp by P = oe exp (- zh) , writing 6 =6a, 


(277)6 ( a + see 
== exp { — —— since & ~— =E, 
vasege FB 


where M is the mass of a nucleon. 
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To obtain the proton energy spectra these expressions for W must be 


_ integrated over the density of states dp, which, since conservation of momentum 


allows the elimination of one variable (k,), is given by 


2 
ide (=) Ae ee ee, (18) 
We may now write | W dpz= aay | ak, & | W dk, 


_the region of integration being k,?+h,?+k,.k, < ME/h?. Thus the spectrum 
_ of the proton energy «, =/?k,?/2M is given by 


P(e )\deron ede, [ W dk. 


_ The proton spectra obtained from the above expressions W), W, and W, and 


normalized to the same value are shown in Figure 1. 


70 80 90 


4 
! Energy (mel) 


Figure 1. Energy spectrum of protons from *He+7-+1H-+-2n : a, scalar meson; 
b, pseudo-scalar meson ; c, ‘ standard’. 


A typical wave function to describe the alternative process of deuteron 
emission is va; 
Ty oF (3 !) i oP wME 7€ I< ig ke 
ijk 


Qa\3i4 
where t=" (=) exp {—ar,?+41k.(r,+1,)—7k.r,}, 


with spatial part J, =%4.;. Integrating the matrix elements for this process over 
the spin variables and summing over all final states with deuteron momentum 
tk gives W,=31,?/2 and W,=3J,?+J.?+(J,—Jz)?, where J; and J; are again 
defined by equations (16) and (17). On integration, the standard expression W, 
for deuteron emission becomes 


2578 / B\3l2 e : aS 

= P5ag (=) exp (—2 ME/h78), again writing B =6a. 

For two-particle emission dp, dE ={V/(27)?}dk; only one deuteron energy is 
permitted and its probability, { Wdp,, is compared with the probability of proton 
emission. The expressions W,, W, and W, above give the relative probability 
of proton to deuteron emission as 23:1 (scalar interaction), 125: 1 (pseudo-scalar 
interaction), and 60:1 (by integrating W, over the density of states) respectively. 

4He can disintegrate in three ways, viz. 


4He+m7->1H+3n, *He+7->7H+2n, *He+n->2H+n. 


Wo 
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The *He wave function is 
4 
Yin = (4 ieee EXD ls Lar) 2 seg PuPiMiE iN eElEijnls seve (19) 
ajkl 
4 
with spatial part ep RESP (ET) ee ee ee (20) 
‘ 1 
A typical final wave function is 
Yaa (4)? 2 opus PaPing iS eSi€ijkl> Teese (21) 
ajk 
with spatial part Mie Wags Le ee ele ge ee (22) 


where for example 


ang = V3? exp {i(k .r,+k,.r;+k,.r,+k,.7,)} for proton emission, 
ijkl p A=) ats | 3 k PV OLG | p 


2a \3/4 ae é 
dee (= ) exp {— a7, + 37k, .(r,+1;) +72(k,.r,,+k,.r,)} for deuteron 
; emission, 
_ _yik 2av/3\3? exp {— a(n? try? +12) + Hk, (t+, +4,) + ik, .r;} sA) 
hee 7 for triton emission. 
For each mode of disintegration the expressions W,, W, and W, corresponding ~ 


to final states in which the charged particle has momentum /k, are evaluated, — 
giving for proton emission | 


Wy =Z\ (ie 13)" U5 — 1s) Ua alge eee (23) 
W.=2{(Ja—Js)?+Js—Ja)?+(Jo—-Ja?}, eee (24) 
2. Sar? 2ME 
Wy=|[ Yea*bal = pasexp(—See)s ves (25) 
writing « =8a and 
4 F2B2 
E=kinetic energy available in final state = = ; 
1 


where, as before, 


T= | tov* Osher with O,,=$(x,) and J,= | par* Onda with O2,=Vyh(1); 


for deuteron emission 


W223 NN OO ee (26) 
We=(2hit+Jat+Ja)P+20s—-Jay, wees (27) 
a ONS IP 2ME 
Wey (5) (:) fy (- =) 2. a (28) 
for triton emission 
Wyi=2 = yh OO ae ee eee (29) 
W=2(Ji—Jalee ne ote ee (30) 
27\6 2ME 
—/p-1 3/20 lee ef 
W,=V-1(12) (F) exp ( =) Dae (31) 


From these expressions the relative probabilities of proton, deuteron or triton 
emission are found to be 

140:15:1 for scalar interaction G, 

450:20:1 for pseudo-scalar interaction G,, 

320: 30:1 by integrating W, over the density of states, 


Tr 


eX Ato ROR 
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It appears that the probability of triton emission is very small and that proton 
emission will predominate. If the deuteron radius is increased by replacing, 
in its wave function, a by a’ (<a) the only result is an alteration in the probability 
of deuteron emission, from ‘He, by a factor (a’/a)?2(2/3 + a'/3a)-*—putting 
a’ =ha reduces this probability by about 40%. The proton energy spectrum 
P(«,) ce"? | W dk, dk, for each expression W,, W, and W, is shown in Figure 2. 
‘The corresponding neutron spectra are shown in Figure 3; these show a stronger 
dependence on the meson interaction but the neutron spectrum would 
presumably be more difficult to determine experimentally than would the 
proton spectrum. 


Cs 
— 
30~—«40~C*«*SO G0. > ao eC 8 
Energy (Mev) : Energy (Mev) 
Figure 2. Energy spectrum of protons from Figure 3. Energy spectrum of neutrons from 
*He+7->1H+3n: *He+7->1H+3n: 


a, scalar meson ; b, pseudo-scalar meson ; ¢, ‘standard’. 


$47 Din INDTPGRA TION OR 26 AND: 0 
The calculations are now extended to cover the disintegration of carbon 
(or oxygen) into two (or three) «-particles and one singly charged particle by 
ascribing an a-particle structure to these nuclei. As the use of totally 
antisymmetric wave functions for these nuclei would involve too much labour, 
in this treatment the interchange of nucleons between «-particles is neglected ; 
this is equivalent to assuming that the time required for the disintegration to 
take place is short compared with the lifetime of an «-particle in the nucleus. 

The wave function used to describe such a carbon nucleus is 


Y= tLe? exp (-25R,7) ean T gi eby (32) 
1 


where R,; is the distance between the centres of the 7th and jth «-particles, 
5 is a constant and ‘’,’ is the *He wave function of the last section (equation (19)). 
The final wave function to describe the disintegration into two «-particles and 
a singly charged particle is 
Hop = Foe eee exp x1-RitxXe-Rs)}, 0 eee (33) 

where fix, and hy, are the momenta of the «-particles and Y,, is given by 
equation (21). In this case Uir, =4R, and Uik,= x3; conservation of momentum 
gives D?y,=0 and reference to the centre of mass gives X}R,=0. Corresponding 
wave functions and relations hold for nuclei containing more than three 
a-particles. The interactions G, and G, are just as before and each matrix 
element 


Fi a 3 
Vel | ¥arG¥er = Vee | exp {—i(x,-R, + X2- R2)} Por* GV’; exp |- aR, | , 
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omitting the «-particle normalization factors. Similarly, for the disintegration 
of a nucleus composed of  «-particles into n—1 free «-particles and a singly 
charged particle, 


2 nm—T n 
H=V-mvr | exp (1. yR DOR Glee ee (34) 
1 1 


This differs from the matrix element for the corresponding helium disintegration 
only in the space functions; thus integration over the spin variables and 
summation over the final states gives the same expressions for W, and W, as 
previously (equations (23), (24); (26), (27); (29), (30)) where J, and J; are now 


expressed in terms of the wave functions 


mI 
by = exp (- ER?) be, An — exp ( = Xie 2) Por: 
1 


Evaluation of W,=|¢y*,|? for the three modes of disintegration considered 
gives 


Wha Vo * Ofor, proton eNussiOn .) » gee eater (35) 
W, = V-™+) {3 ./(a7)}-*O for deuteron emission, ...... (36) 
Wy= Vat 7 ri/44/3) QO for triton, emission, (jeans (37) 
n42-78 
where Os aoe = aus —yx}, 
no.( na(2ca)" ha 
Pe enrx (38) 
=i, fl ain es oe 
e < gME~ Une 
1—c 2ME 

and a tae a ae et pe ve eo ey eee bee Wetton (39) 


This important factor W, is no longer constant as in the helium case but depends 
on the way the meson rest energy is shared between the kinetic energy of the 
centre of mass of all the «-particles and the kinetic energy relative to its centre of 
mass of the split «-particle. ‘The extent of this dependence is governed by the 
parameter y which is a function of the constant ¢ occurring in the initial wave 
function. O<c<1: the lower limit corresponds to a set of free «-particles; 
at the upper limit the spatial part of the wave function degenerates into 
exp [ —(a/n)X}"r,;7] which cannot be antisymmetrized since all 4m nucleons are 
then equivalent. The wave function is not a sufficiently good approximation to 
allow the value of c to be determined from the nuclear radius. It is easily shown 
that the standard energy spectrum (i.e. the spectrum obtained from W,) of an 
a-particle emitted in one of the above disintegrations is given by 


(1—?) m 
Plc) dee (1 — 1)" #12 dt | ‘ e«(1- a usrl2-4 dy, os... (40) 
J0 yl at) 


where 0 <¢=ne/(n—1)E <1, « is the energy of the «-particle, nis the number of 
a-particles in the initial nucleus And fi) | 2-208 wy/Z according as a proton, 
deuteron or triton is emitted. The value of EF, and hence of y, varies slightly with. 
the particular disintegration considered: in the case of proton emission from 
carbon c=1/2 gives y212. The «spectrum for this disintegration is shown 
in Figure 4 for various values of y. For y=12 the «-spectra of W, and W, have 
also been calculated and are found to be indistinguishable from one another 
numerically: these and the corresponding standard spectrum are shown in 
Figure 5, together with a comparison spectrum of W=constant. Using the same 
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value of y (=12) for carbon, oxygen and neon the standard «-spectra of the 
disintegration of these elements with proton emission are, numerically, almost 
identical. The carbon and oxygen «-spectra for c=1/2 are found to be in fair 
agreement with the experimental results of M.M.R. The relative probabilities 


a 


— ame 
0 5 eS ery 9 0 10 


20a asc 40 
Energy (Mey) Energy (Mev) 
Figure 4. ‘Standard’ energy spectrum Figure 5. Energy spectrum of «-particles 
of «a-particles from 1'C+2-~> from @C+2->2*He+1H+3n with 
2 He-—"H + 33a, y=6; 6, y=12 ; y=12: a, scalar and pseudo-scalar 
c, y= 24. mesons ; b, ‘ standard’ ; ¢, statistical 
factor. 


qQ 30 60 90 120 150 180 
Degrees 


Figure 6. Angular correlation of «-particles from 1#C-++7-—>2 *#He+1H-+3n: a,b,c, scalar meson 
with y ~0, 12, Oo respectively ; d, pseudo-scalar meson and ‘ standard’ (all y). 


of proton, deuteron or triton emission given by W) for yproton=12 are: 
carbon 118:16:1, oxygen 67:13:1. The dependence of these ratios on the 
value of y and on the type of interaction considered is shown in the Table, which 
gives the probability of proton emission relative to that of triton emission from 
carbon (in processes in which two «-particles are also emitted). 


y (0) 6 12 24 (00) 
scalar (54) 78 102 153 (187) 
pseudo-scalar (50) 102 174 238 (343) 
‘ standard ’ (17) 68 118 169 (205) 


When more than one charged particle is emitted in a disintegration the angle 
between two such particles can be measured and the distribution of this angle 
may be found. In the disintegration 


Cee Tee 2 lied Tt oi 2 OP ee (41) 
the angular correlation of the two «-particles is given by 
N(6) do =d0 | W don", ) 
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where dp,” d0 =dp, and @ is the angle between the two «-particles. W,) and W, 
are independent of the distribution of kinetic energy amongst the three 
x-particles and this leads to angular distributions identical with that for a constant 
matrix element. W, depends on the kinetic energy of the centre of mass of the 
split «-particle and yields an angular distribution somewhat different from the 
standard case and depending, slightly, on c. Figure 6 shows the variation of 
N,(6) with c and the standard N,(@) for this disintegration. Preliminary 
experimental results, communicated by Rochat (private communication), are 
as follows: 

Angle @ (deg.) 0-30 30-60 60-90 90-120 120-150 150-180 

No. of events N(@) 5 12 7 6 11 1 


§5. INDIRECT DISINTEGRATION OF ¥C 

The scheme (41) is not the only one which can describe the disintegration 

of 2C on a simple «-particle model; an intermediate state may be involved, 
represented b 
P y 2C 4 7-—>8Be*+4H+3n, | 
®Be*—> 2 *He. j 


An analogous scheme has been suggested by Telegdi and Ziinti (1950) to explain 
the experiments on the photo-disintegration of ?C into three «-particles. They 
obtained agreement with the measured «-particle energy spectrum by assuming 
that the only state of *Be involved was the excited state at 3mev. With the same 
assumption it is possible to calculate the energy spectrum and angular correlation 
of the «particles produced by process (43). The results of the previous section 
show that it is possible to neglect the meson interaction and also the spins of the 
nucleons when calculating the «-spectrum for process (41). This should also 
be true when an intermediate state is involved. 
Taking the wave function of the intermediate state to be 


4 
br=V exp (12kj.t,4ik-R) SOW sevens (44) 


where 7k is the momentum of *Be*, R its centre of gravity and 4(@) its internal 
wave function ( is the distance between the «-particles composing *Be*), leads to 


W, =7®N2V-4(6co4)-8 exp (- 2A een tase i) 32 atts (45) 


ha So unc 
where N= | d(e) exp (—cap?) de. . 


Choosing the value c=1/2, the probability of process (43) relative to that of (41) 
is found to be 0-1 N*a??, The value of NV depends on the spin of *Be*, which 
according to Wheeler (1941) is 0. ‘Thus, taking the *Be* wave function as 


P(e) =(u/7)” exp(—up’), wenn (46) 
the relative probability is of order 20°% for reasonable values of u(~«). More 
recently, scattering experiments (reviewed by Hornyak et al. 1950) indicate 
that SBe* has spin 1 or 2; our.assumptions would then give N=0. ‘This does 
not mean that process (43) is forbidden, for slight changes in the model adopted 
would give non-vanishing values of N. The relative probability should however 
be much less than the 20° found for spin zero. 


The Disintegration of Light Nuclei by Meson Capture 1075 


The distribution of energy « of the «particles is found to be 
PO) =de|e@o de, (47) 


the limits of integration being v=(2M)-¥?|\/(U/2)+/e|, where U is the 
- excitation energy of the *Be*, wv is its velocity and g(v) dv its velocity distribution. 
Figure 7 shows this «spectrum for c=1/2 and c-1; agreement with experiment 
is only obtained for the latter case. The angular correlation of the «-particles 
from this process is shown in Figure 8. 


y 
0 5 10 5 20 : 
Energy (Mev) 0 30 oo, 150180 
Figure 7. ‘Standard’ energy spectrum of Figure 8. ‘ Standard’ angular correlation of 
x-particles from 12C+7-—>*Be* +1H+3n; «-particles from 12C +7 8Be* +-1H+3n; 
SBe*—2 *He: a, y=0; 6, y=123 c, y=24. ®Be*—2 *He: a, y 2/0; 6, y=12. 


§6. DISCUSSION 


In the previous sections the energy spectra of particles emitted from 
the nuclei *=He, *He and #C as a result of z-meson capture have been 
determined for the meson—nucleon interactions %(r) and o.V¢(r). These 
were compared in each case with a standard spectrum, which is independent 
of the meson interaction and depends only on the nuclear wave functions 
assumed and on the statistical (phase-space) factor. The standard spectrum 
neglects the antisymmetry properties of the wave functions, which would cause 
a reduction of intensity at the high-energy end of the spectrum with a 
corresponding increase of intensity at low energies : such a deviation from 
the standard is noticed in all other spectra. 

In the cases of 3He and 4He the standard matrix element is a constant, 
independent of the way the kinetic energy is shared amongst the emitted 
particles, so the standard spectrum is just the statistical one. ‘The dependence 
of the proton spectra of ?He and 4He (Figures J and 2) on the meson interaction 
is probably not sufficient to enable experiment to distinguish between the two 
interactions. Observations on the neutrons emitted from *He, whose spectrum 
(Figure 3) depends considerably on the meson interaction, would probably be 
necessary. 

In the case of #2C, which was treated on an «-particle model, the standard 
spectrum (Figure 4) depends on the parameter y, which represents the degree 
of binding of the «-particles in the nucleus. (The limit y=0, representing such 
strong binding that the «-particles are no longer distinguishable, gives the 
statistical spectrum.) y may also be used to determine the radius of the carbon 
nucleus. For a nucleus composed of n «-particles, by integrating the square of 
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its wave function over the coordinates of all but one of the nucleons, the density 
of nuclear matter is found to be 


4ncar? 


p(r) oe exp (- a) = exp (—277/R,,”), ce ceee (48) 


defining R,?=(n—1+3nc)/2ncn. Thus R,,/R,={(n—1+3nc)/3nc}?. If we 
use the value c=1/2, which gives agreement with the experiments of M.M.R., 
we find R,/Ry,= 1:20, compared with the value 1-44 obtained from the formula 
R,,«<A'% cn, This rather poor agreement must be attributed to the fact that 
gaussian wave functions cannot represent simultaneously both the nuclear radius 
and the average kinetic energy of the «-particles in the nucleus. 

The «-particle energy spectra (Figure 5) from !C (or 1®O) do not distinguish 
between the two meson interactions. This is not very surprising since the 
form of the matrix elements (equation (34)) shows that the interaction only 
affects the momenta of the free nucleons relative to the split «-particle and this 
can be expected to reflect little in the energy spectra of the other «-particles. 
The deviation of these «-spectra from the standard spectrum may be considered 
as due solely to the antisymmetry properties of the wave functions. In the 
centre-of-mass system of the split «-particle the energy of the proton (or neutron) 
relative to the total energy in this system has the same spectrum as in the *He case 
(exactly in the case of the standard and pseudo-scalar spectra and approximately 
in the case of the scalar spectrum). ‘The proton energy spectrum in the usual 
system (centre of mass of all nucleons at rest) requires the inclusion of integrations 
over the momenta of the «-particles and is expected to show even less dependence 
on the meson interaction. 

The pseudo-scalar interaction and the standard matrix element give the same 
distribution of the angle between the two «-particles emitted in the direct 
disintegration of #C (Figure 6) and this is independent of the parameter c. 
The angular distribution due to the scalar interaction depends slightly on ¢ and 
differs somewhat from the standard. 

In §5 the possibility of the indirect disintegration of ”C involving the 
3 Mev. excited state of *Be was considered. ‘The (standard) «-particle energy 
spectrum obtained on this assumption (Figure 7) agrees with the experiments 
of M.M.R. only for c= 1, which is the strong-binding limit of the «-particle 
model used. ‘Thus the indirect process alone is insufficient to account for the 
observed energy spectrum which may, however, be reproduced, within the limits 
of experimental error, by increasing the value of c beyond 1/2 and by assuming a 
corresponding proportion of indirect splitting. The (standard) angular correlation 
of the «a-particles emitted by this indirect process (Figure 8) depends 
considerably on the value of c. The difference between the angular distributions 
of the direct and indirect processes should be sufficient to allow an experiment 
to determine the relative amounts of the two processes and the corresponding 
value of c. Preliminary measurements by Rochat (private communication), 
quoted above, suggest that many of the observed disintegrations follow the 
indirect process. As indicated in §5 the simple «-particle model used here is 
not sufficient to estimate the probability of the indirect process since it gives 
zero unless the spin of *Be* is zero, which is probably not the case. It would be 
desirable to have an independent estimate of the relative probability of the two 
processes: this, used in conjunction with the observed energy spectrum, would 
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determine the value of c required. With this information the angular correlations 
and the proton and neutron energy spectra could be determined uniquely. 
The values of the relative probabilities of emission of the various charged 
particles cannot be relied upon in the same way as can the spectra since, although 
these values are independent of the intensity of the meson interaction, and the 
term exp(—2ME/h?«) which appears in the total probabilities is eliminated, 
the ratio of proton to triton probabilities, for instance, involves a factor 


_« %o Ry.° and the radius of the *He nucleus is not known very exactly. 


It is interesting to note that the use of wave functions for 9H and 4He of the 
type suggested by Irving (1951), which have approximately the correct 
asymptotic form, leads to nearly the same relative probabilities as (31 a). Putting 
n=( and fixing « in each case by the binding energy (m and « are defined in 
Irving’s paper) and neglecting the meson interaction, the relative probabilities 
of triton, deuteron and proton emission from 4He are 1: 20: 320. 

As pointed out by Yamaguchi (1950) it is possible to apply Fermi’s (1950) 
theory of high-energy nuclear events, which Fermi himself describes as the 
opposite extreme to the perturbation theory approach, to the disintegrations 
here considered. This theory supposes that energy, here equal to the mass of the 
m-meson, is suddenly released in the nuclear volume Q and that statistical 
equilibrium is reached. This gives the square of the effective matrix element 
for s emitted particles as (Q/V)*-1, where V is the normalization volume and 
momentum is conserved. “This is independent of any interaction and may be 
compared with the standard expression W, of this paper, which is an 
approximation, independent of interaction, to the square of the matrix elements 
given by first order perturbation theory. The Fermi expression (Q/V)s~1, being 
independent of energy, would give purely statistical energy spectra in all cases. 
In the cases of ?He and 4He therefore, the energy spectra obtained by the two 
methods agree. It has been shown by M.M.R. that the observed «-spectra of 
12C, 14N and 1%O do not agree with the corresponding statistical spectra, whereas 
they may be explained by the perturbation method. ‘The relative probability 
of proton to triton emission given by Fermi’s method depends on (2? oc R® as 
in the perturbation method, but using Q =47R?/3 with R= A" x 1-4 x 10-8 cm. 
gives equal probabilities of proton and triton emission from both *He and #C. 

We conclude that for a very light nucleus (*He) the disintegration is 
practically independent of the details of the nuclear structure but depends on the 
meson interaction to a certain extent. For slightly heavier nuclei (’C) the 
disintegration depends strongly on the nuclear structure and only weakly on the 
meson interaction. The disintegrations of heavy nuclei (1°Ag) can, as is well 
known, be described by the evaporation theory and are independent of both 
meson interaction and nuclear structure. 


POPPE IN DIX 


All the matrix elements occurring in this paper can be expressed in terms 
of integrals of the type 
i nv n 
I, =| Oe exp - Ub r;.r;—-1Up;.t;} ar... Ds PAE (49) 
td ) 
Here we have put D7" 1m,x; =0 and p, =k, —k,, ,,m,/m,, ,,(m, = mass of 7th particle ; 
for He, m,;=nucleon mass M, for C, m;=M or 4M). We have omitted the 
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constant factor due to integration over the centre-of-mass coordinate and 
require X7t'k,=0. The equality 
| exp {—br? —ip.r} dr =(7/b)°” exp (—p?/40) 


can be applied successively to the integrations over r,....f,_; leaving an integral 
of the form : 
| OF (=) exp { Wi nha <i try, :: r,)dr, 


and it is obvious that «, and A,, do not depend on O(r). J,, can be expressed 
simply in terms of the symmetrical determinant B=|),;|, with cofactors B,,. 


- qonl2 me (ae? fO,(r) exp (— Cla in, r)dr 
y= aa? 135 But-Ps | teapot mania fo oe 
where An =) Dae - and 7 = 2B Pil Ban: 


The final bracket of (50) has the values: 1, standard; —3/2«,+A,?/4a,”, scalar; 
7A,,/%,, pSeudo-scalar. ; 
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ie Range-Energy Relation for Protons in Ilford C2 Dry Emulsion 
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_ ABSTRACT. The range-energy relationship for protons of energies up to 12 Mev. in 


- Ilford C2 dry emulsion has been determined with an accuracy better than 0:5°% in energy. 
The particles were produced by bombarding a thin beryllium foil with deuterons having 


an energy of 8 mev. The energy of the incident deuterons was obtained by a method 


which involves the measurement of the two angles of emission at which the two proton 


_ groups of longest range from the reaction °Be(d, p)!°Be have the same energy. 


§1. INTRODUCTION 
HE range-energy relation for protons in Ilford C2 dry emulsion has 
| recently been measured by Rotblat (1950, 1951) and by Catala and Gibson 
(1951) by making use of the elastically scattered deuterons and disintegration 
protons resulting from the deuteron bombardment with various gas targets. The 
energy of the primary deuterons was determined (Rotblat 1950) by comparing, 
under the same conditions, the ranges of alpha-particles from polonium, 
thorium C and thorium C’ with the ranges of alpha-particles scattered from the 
collision of the deuterons with helium at certain angles of scattering. This 
work was recently corrected (Rotblat 1951) in the light of a more accurate 
range-energy relation for particles in air (Bethe 1950) which affects the loss of 
energy of both the incident and emitted particles in passing through the gas 
filling the camera. 

The present measurements were she independently during an investi- 
gation of the disintegration products from a thin beryllium foil exposed 
under vacuum conditions to the deuteron beam from the Liverpool 37 inch 
cyclotron. Although the results are in agreement with those of previous 
workers, the method of procedure is sufficiently different to make this report 
significant. The (d,p) reaction with beryllium is particularly suited for the 
determination of the range—energy relation for two reasons: (a) because of the 
large variation of the proton energy which can be obtained merely by varying 
the angle of observation, (b) because the Q-values for the two proton groups of 
longest range have been accurately determined by Buechner and Strait (1949) 
using a magnetic analyser. 

The energy of the primary deuteron beam was determined by a method which 
involves the measurement of two angles of emission at which these two groups 
of protons have the same range. The range-energy relation has been obtained 
for protons having energies up to 12 Mev. 


§2. EXPERIMENTAL 
The camera used in this experiment has been described by Rotblat, Burrows 
and Powell (1951). The target was a foil of beryllium of thickness 0-325 mg/cm? 
prepared by evaporation. ‘Two exposures were made, one with the foil at an 
angle of 45° with the deuteron beam for measurements at forward angles and 
70-2 
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the other with the foil at 135° for measurements at backward angles. The plates _ 


used were Ilford type C2 with an emulsion 200» thick. The emulsion was dried 


by keeping the plates in vacuum for about 2 hours before the experiment as well | 


as throughout the exposure. The plates were examined with a binocular 
microscope having a 96 x oil immersion objective and 6 x eyepiece. Histograms 
for the range distribution for 23 different angles ranging from 174° to 150° were 
plotted for all accepted tracks of elastically scattered deuterons, disintegration 
tritons and protons. About 2,000 tracks were measured for each angle in a small 
area of width corresponding to a maximum variation of +30’ in the angle. 


5 IRIS SOME Ats) 


A representative spectrum for the disintegration particles at angle 30° is 
shown in Figure 1. An analysis of the histogram shows that groups 1 and 3 
are due to protons emitted in the reaction °Be(d, p)!°Be, in which the 1°Be nucleus 


30 


Number of Tracks 
8 


300 450 600 750 
Range in Emulsion (1 scale division = I-13) 


Figure 1. Histogram of disintegration products of beryllium bombarded by 
7-7 Mev. deuterons at 30°. 


is left in the ground state and first excited state respectively. Groups 2 and 4 
are due to protons from the reaction 1O(d, p)!”O with the residual nucleus 
formed in the ground state and first excited state. Group 5 is due to tritons 
from the reaction *Be(d, t)*Be with the ’Be nucleus left in its ground state. 

The main target contamination was oxygen in the form of an oxide layer. 
Since the peaks in the spectrum due to beryllium change position with angle 
more rapidly than those due to contamination, it was usually possible to separate 
the two components. With the elastically scattered deuteron group and the 
proton group p, it was found advisable to fit their range distribution with gaussian 
curves from which the mean range and the half width were determined. In 
the case of the proton group py, the range was so long that it was far from 
disturbances due to target contamination and the mean range was calculated 
directly as the weighted mean of the group. These measured ranges are the 
projections of the true ranges on the plane of the surface of the emulsion. 


| 
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The true range of a particular group was calculated from the known mean 


__ angle of approach of the particles belonging to that group to the surface of the 
emulsion. 


The energies corresponding to these groups can be calculated from the angle 
of emission if the Q-values and the energy of the primary beam are known. As 
mentioned before, the two proton groups from beryllium have been magnetically 
analysed by Buechner and Strait (1949), who give accurate Q-values for these 


_two groups as Q)=4:576+0-012 Mev. and Q,=1-201+0-007 Mev. The energy 


of the incident deuterons was calculated from the fact that the range of protons 


__ belonging to the ground state of the reaction at large angles is the same as the 


range of the protons from the group due to the first excited state at small angles. 
It was found, for example, that at an angle of 120° the group p, had a range of 


_ 428-81 which was nearly equal to the range of the group p, (429-1) emitted 


at an angle of 17°30’, i.e. the energy of both groups could be taken to be the 
same. From the two known equations in Q) and Q, the proton energy could 
be eliminated by a simple calculation leading to a second-power equation which 
could be solved to find the value of the primary deuteron energy at the effective 
centre of the foil) The problem was more carefully treated by plotting the 
variation between the range and cos@ for the two groups in this small region. 
By selecting the best line which passes through these points one can find the 
cosines of the angles of emission for the two groups of protons of the same range, 
i.e. of same energy. It should be noted that the Be foil is sufficiently thin 
(~0-4 cm. air equivalent) for the mean energy of a group of protons to be 
calculated with sufficient accuracy on the assumption that the nuclear collisions 


_ occur on a plane equidistant from both surfaces of the foil. ‘Thus, no correction 


was necessary for the loss of energy of protons escaping from the foil: this was 
approximately the same for the two groups since one was measured with the 
forward exposure arrangement and the other with the backward arrangement. 
The agreement between the results for different values of the common range 
are shown in Table 1. 


Table 1 


Common range (microns) 429 425 420 415 410 405 400 
Primary deuteron energy (Mev.) 7:696 7:702 7:705 7:700 7:695 7:700 7-699 


From the above results, the mean value of the energy of the deuteron beam 
at the centre of the foil was taken as 7-70 Mev. and used to determine the energy 
of the emitted particles. The accuracy of this method depends mainly on the 
uncertainty in the Q-values of the two proton groups. ‘The energy of the incident 
deuterons was thus determined with an accuracy of within +0:3%. 

For the determination of proton energies for constructing the range—-energy 
curve a correction was made at each angle of observation for the loss of energy 
of a proton in escaping from the foil. The stopping power of the foil material 
was calculated from the shift in the elastic deuteron group at 30° when the foil 
was in the forward and backward arrangements. The value adopted was the mean 
of this value and that given by Livingston and Bethe (1937), namely 0-36 cm. air 
equivalent for 8mev. deuterons. The uncertainty in the stopping power will 
give rise to a very small error since the energy loss of protons in the foil was 
only about 0-02 Mev. 
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The results are shown diagrammatically in Figure 2 (a) and (b) and represen- 
tative data from the mean curve are given in Table 2. The experimental points 
were obtained mainly from the range measurements on the elastically scattered 
deuterons and the two groups of disintegration protons from the reaction 
*Be(d, p)°Be. A few points between 3-72 and 3-87 Mev. were obtained from the 
reaction °Be(d,t)’Be the Q-value of which is 4:67+0-03 Mev. (El-Bedewi, 
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Figure 2 (a). 
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Table 2 
Range Energy Range Energy Range Energy 
(microns) (Mev.) (microns) (Mev.) (microns) (Mev.) 
5 0-49 80 Sole 200 5-495 
10 0-80 90 55355 225 5-905 
15 1-045 100 3°58 250 6°30 
20 127 110 3-80 275 6-665 
25 1-475 120 4-005 300 7:01 
30 1-67 130 4-22 325 HES 
35 1-85 140 4-42 350 7:69 
40 2°01 150 4-605 375 8:01 
45 PEAVY 160 4-795 400 8-325 
50 27325 170 4-975 425 8-635 
60 2°61 180 5-150 450 8:94 
70 2°88 190 beS2 475 LDS 


Energy (Mev) Scale for Curve IV 


Range Energy 
(microns) (Mev.) 
500 9-51 
525 9-78 
550 10-04 
575 10-30 
600 10°55 
625 10-805 
650 11-05 | 
675 11-30 | 
700. mal -545 | 
725 11-785 | 
750. 12-02 
775 12:26 


ft 


i 


= 
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\ unpublished). The point at 0-584 Mev. was taken from measurements on the 


protons from the 14N(n, p)!#C reaction which were recently reported by the 
author (El-Bedewi 1951). The accuracy at any point above a proton energy of 
0-5 Mev. on the range-energy curve is estimated to be better than 0-5%% in the 
proton energy. 
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The Differential Cross Section for High Energy Nucleon- 
Nucleon Collisions and the Mean Square Angle of Scatter 
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School of Theoretical Physics, Dublin Institute for Advanced Studies 


MS. received 25th Fune 1951 


ABSTRACT. In an attempt to account theoretically for the lateral and angular develop- 
ment of the penetrating particles in high energy cosmic-ray showers, we have developed a 
method for deriving the differential from the total cross section, in the laboratory frame of 
reference. An integral equation is set up for the differential cross section, and solved 
assuming that the differential cross section in the centre-of-mass frame of reference is of the 
form R-+S’ cos? 6, where R and S ’are functions of the incident and scattered energies only. 
Using a form of the total cross section derived from the analysis of high-energy cosmic-ray 
experiments, the mean square angle of scatter in a nucleon—nucleon collision has been 
calculated and found to be inversely proportional to the energy of the incident particle. 
The results are compared with experiments. 


§1. INTRODUCTION AND PRELIMINARIES 
OSMIC rays are almost the only reliable source of information on the force 
between two nucleons whose relative velocity approaches closely to that of 
light. To frame a tentative theory in these conditions is as difficult as to 
realize them in the laboratory experiment. _It is easy to see that such information 
—badly needed in the theory of nuclear forces—may well be obtained from investi- 
gating the lateral spread of extended air showers, more especially of their nucleon 
component. A reasonable assumption for the differential cross section that gives 
the observed lateral spread in agreement with observation up to such high relative 
energies will probably be fairly correct. The chief requisite for determining the 


* Both authors now at Department of Mathematical Physics, University of Adelaide, South 
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lateral spread from any assumption for the cross section is the mean square of the 
angle of scattering after one collision. 

The problem of the nucleon cascade differs from the theory of the lateral 
spread of the electron—photon cascade, developed by Euler and Wergeland 
(1940 a, b), Moliére (1942) and others, in two essential points. In the latter case 
there is only electromagnetic interaction, which in the collision between nucleons 
is outweighed by the nuclear force proper, and, for high energy, is even negligible. 
Secondly, with electrons inelastic scattering may be neglected, while a high-energy 
nucleon must lose an appreciable amount of energy by meson production in almost 
every collision. 

This paper is concerned with the two chief preparatory steps, the first being the 
most important: to calculate the differential from an arbitrary total cross section in 
the laboratory frame and then to infer the probable form of the differential cross 
section from previous work on the nucleon cascade which required only the total 
cross section as a function of the energies and accounted satisfactorily for the 
vertical distribution of nucleons in the atmosphere, both in number (Messel 
1951a,b, Messel and Ritson 1950, Janossy and Messel 1951) and in spectral 
distribution (Messel 1951c). The first assumption we have made is that at 
energies greater than 3 x 10? Mev., with which we are mainly concerned, elastic 
collisions may be neglected, in other words that a meson is always produced even 
if it does not escape from the nucleus with which the nucleon collides. The 
second step is merely the computation of the mean square angle of scattering from 
the differential cross section obtained. 

Faced with the problem of finding the differential cross section from a total 
cross section given in the laboratory frame of reference, we have made the further 
assumption that the angular dependence will be simplest in that frame of reference 
in which the mass-centre of the two nucleons is at rest before the collision. More 
especially we assume that, in this frame, the expansion in powers of the cosine of 
the scattering angle 6 may be terminated at cos? @; since it is meaningless to distin- 
guish between the two nucleons after the collision odd powers of cos@ cannot 
occur. ‘The differential cross section must then be of the form R+5S’ cos? 6, 
where R and S’ depend only on the energy of the incident and scattered particle. 
It is remarkable that both these functions are uniquely determined by the total 
cross section and can be expressed in terms of the latter without specializing it. 
The reduction to the laboratory frame is then obvious. 

The most convenient unit of mass for our work is the rest mass of the proton; 
if, then, the unit of time is so chosen that the velocity of light is unity, the rest 
energy of the proton (~10 mey.) will also be unity. The rest mass of an arbitrary 
scattered particle in these units will be denoted by m, thus, if the particle is a 
proton m=1, if it is a 7-meson m=1/6:5. The momentum and total energy of 
the particle in the laboratory frame will be denoted by Pand U respectively, so that 


UAL Pree Pee tel et Ane (1) 


and all quantities referring to the centre-of-mass frame are distinguished by a bar, 


e.g. U and P; since the momentum vector P and the energy U together form a 
covariant four-vector, then 


U=U,U— P,P cos6, Berita C4) 
Pcos0= U,Pcos6— P,U esol) Psin6= Psin 6, AE?) 
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¥ where P and U, are the momentum and energy (which, since m=1, is just the 
/ coefficient of the transformation) of the incident nucleon in the centre-of-mass 
6 system; 6 and 6 are the scattering angles in the laboratory and centre-of-mass 
| frames respectively. The relations inverse to (2) and (3) are 


eeu =U,U+ P Pecos) (5); Pcos@=U,Pcos6+P,U. ...... (6) 


__ The momentum P, and energy U) of the incident nucleon are found in terms of 
P, and U, by substituting U=U), P=P,, and 0=@=0 in (5) and (6), thus: 


Me U, +P, =2U,—1- ...(7); P20 P= 20 (0, tye) 


+ 


The requirements of conservation of energy and momentum impose certain 
restrictions on the energies and momenta of the outgoing particles. The energy 
__ of an outgoing particle whose rest mass does not exceed two proton masses cannot 
, exceed the energy of the incident particle in the centre-of-mass system, that is, 


eae ahd meet Sige Hern (9) 


The proof of (9) is given in the Appendix. 
The inequality (9) imposes an important restriction on the angle of emission of 
the scattered particles in the laboratory system. For it follows with the help of (2) 


_ that UgUi— P, Pcose= Oe £2) cos LU 1) Ps Paw 2a. (10) 


The right-hand side of (10) is very nearly 1— U-1, since Py) U, and P2U. 
_ Hence the angle of scattering cannot exceed the value (2U—!)!” radians, approxi- 
mately, where U is the energy of the scattered particle. Scattered particles are 
_ limited to a cone which becomes narrower and narrower as their energy increases. 
Some cosmic-ray photographs show what appears to be a double cone; it is quite 
plausible that the inner cone be composed of particles from primary collisions, and 
the outer cone of particles from secondary or tertiary collisions, of nucleons within 
the nucleus. The latter will naturally have a much lower energy than the former. 

We shall repeatedly make use of the Lorentz invariance of certain differentials. 
For any particle, the volume element of momentum space P?dP dc d¢ divided by 
the energy U is known to be invariant.* But since UdU=PdP, it follows that 
PdUdcd¢ is Lorentz-invariant; and, since db=d¢, 


PdUide=PdUdeh ae ns ius ee (11) 


Also, if W(U, c,¢)dU dc d¢ is the probability, in a particular experiment, of 
finding the particle in question within the limits U, U+dU, c, c+ dc, ¢, 6+4d¢, 
this probability must be Lorentz-invariant; thus, if W does not depend on 4, 


WU,c)dUde=W(U,c)dUde. saa (12) 
Dividing (12) by (11) we have further 
WU,c)/P=WU,/P. aaa (13) 


Similarly, if W(U',c’; U",c’) dU’ dU" dc' dc” is the probability of finding the 
two nucleons, arising from the same collision, in the specified intervals, 
W(U’,c’; U",c”)/P’P” will be Lorentz-invariant. For our present purpose, we 
have found it unnecessary to determine W(U’,c’; U",c"): on account of the 
symmetry of this function with respect to the variables U’, c’; U",c" we require 
only the function W(U’,c')= Jf W(U',c’; U",c") dU" dc", which is identical with 
that defined above. : 


* ¢ stands for cos @ and ¢c for cos 8. 
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In the following determination of the differential cross section for the inelastic | 
scattering of two nucleons, no distinction will be drawn between neutrons and | 
protons. It is, admittedly, likely that the exclusion principle will make the cross ' 
section for p—p or n—n collisions somewhat different from that for p—n and n—p 
collisions. ‘The reason why we do not need to distinguish between neutrons and 
protons is that the total cross section which we use is averaged over all four types of 
collisions. The differential cross section derived will therefore be automatically 
averaged in a similar way. 


§2. DERIVATION OF DIFFERENTIAL FROM TOTAL CROSS SECTION 
Writing F(U) for our total cross section in the laboratory frame of reference, 

we must have jG) de EU). ee (14) 

Note also that PEC) dU 1h tab 9 asa Pee ee (15) 


It will be convenient to change the integration variable from c to U, keeping U 
constant, and to write: 


U=mcoshu, U=mcosha, U; =cosh®, \U—coshZy 
P=msinhu, P=msinha, P, =sinhv, P,)=sinh 2v ; 
Then, using relations (2), (5), (9) and (13) we find 


-cosh v wo a 7} 
| W (m cosh u, costae Paves ha 2 ooh is egeh *) dymavshid) a0sh aint sinh v F(mcoshu). 
m cosh (v—u) sinh v sinh @ m sinh u 
er ee (17) 
We assume however, (see $1) that 
WO JERU) HS(U)(U Pe ee ee (18) 
so 
— _ coshu—coshvcosh @ ¥ 
W (m cosh i, Ee) = R(m cosh cD) 
+ S(mcoshi)(m coshu—mcoshvcosha)*, ...... (19) 
hence writing 
R(mcosht)=r(u); S(mcoshi)=s(a); F(mcoshu)=f(u), ....(20) 
we get 
cosh a ss 5 oy C(t COStia) Se 
| Heine {r(i) + m*s(a%)(cosh u — cosh v cosh z)?} porary iar f= sinh v f(u). 
HUT: (21) 


The value of u lies between 0 and 2v; anda further integral equation is obtained 
by replacing u by 2v—u in (21), thus: 
~cosh v S ous _,,, d(m cosh #) 
{r(w) + m*s(a)(cosh (2v — u) — cosh v cosh @)?} Sates 
=sinhvf(2v—u). ...... (22) 
Our task is to solve (21) and (22) for r(a@) and s(#). One may proceed as follows: 
subtracting (22) from (21), and writing 


f (u) — f (2v—u) 


| m cosh (v—u) 


sy 4m cosh v'sinh.(u=—®) Ad au) ee ee (25) 
reson : ., Umcoshu 
one finds | eae ma coe (u— a — cosh @} a = 2) 6 eee (24) 
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be One may now substitute this for s(z) in equation (21) and after a straightforward 
' but lengthy calculation find 


H) 7(u) = — sinh* v(d/da){m g'(v — a)sinh a} — 4 sinh of f'(v +a) —f'(v—a)}. 
ae +tcoshvif(v+a)—f(v-a)}. se (26) 


__ Asa check on the above, we substitute the values found for R and S in (18) and 
integrate over both c and U; we find the value of the integral to be equal to 1, as it 
_ should be, since our probability function F(U) is normalized to 1. 
Thus, starting from the general total probability function F(U)dU we have 
without specializing it found the differential probability function W(U,c) dUdc. 
| It is convenient to express W(U,c) in the centre-of-mass frame of reference, 
thus using (23), (25) and (26) 


sinh? vw sinh? v coth “@ 


{f"(e+ a) —f"(o—a)}— 


W(mcosh 4, ¢) = 


(f(o+a)+fi(e—a} 


4 cosh vw 4 cosh v 
ae ; ne sae sinh? v csch? @ *, 
ksinho{f'(v+a)—f’'(v—a)}+ Geren A tts coshv 
1 


x {f(v+a) —f(v—a)} +c? sinh? a sinh? v | - 


4 coshv sinh2a@ 
| 3 cosh @ 


x {f"(o+ a) —f"(v—a)} + neg (+a) tf(e—a)} 


4 cosh v sin 


ps 2 sinh? 4 +3 2 e 
Sear ae eee 


The W(U, c) may be obtained from (27) by using the relations (2), (5), (7), (8) and 
13): 
' We now introduce the special form of f(z) used by Messel (1951 c) 

f(@) =20m(U, — 1)-5(cosh @—1){U)+(m—1)—mcosha). ...... (28) 


This cross section appears to account for the development of the high-energy 
nucleon component of the cosmic radiation throughout the atmosphere both from 
a qualitative and from a quantitative standpoint, and even though it is probably 
not exact, it is unlikely that, when the true total cross section for high-energy 
inelastic nucleon-nucleon collisions becomes known, it will differ radically from 
(28). 

In the ultra-relativistic region where U~ P and U,) = P, the expression given 
by (27) simplifies a great deal. In the centre-of-mass system in this approximation 


W(m cosh a, c) ~ 30m? exp (—5v) exp 3a (expu—mexpu)(1—c?). ...... (29) 
It should be noted that the expression given in (29) is valid either for nucleons or 
mesons. However, since the cross section (28) has not been checked by seeing 
whether it can account for the experimentally observed behaviour of the meson 
component of the cosmic radiation, it should be treated with some caution for 
mesons. 
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§3. THE MEAN SQUARE ANGLE OF SCATTER 
The value of the mean square angle of scatter in a nucleon-nucleon collision is 
given by (@)=[ftan2@W(U,c)dUde, eens (30) 
where we have replaced 6? by tan?@. We prefer to carry out the above integration 
in the centre-of-mass system. Writing 


tan?@=sin26/U,2(cos0+ M)?, sae (31) 
with M=P,U/U,P er (32) 
-cosh v rl 1 os C2 ; " 5 
BN = Barony Re 2 of \72 5) dc. 
(6?) I, ie cosh? a (¢+ My {r(a) + m? sinh? v sinh? @ s(a)c*}d(m cosh a) dc 
vaste (33) 
The c integration may be carried out immediately, yielding 
cosh » r() M-+41 
a \ eet eh 
<6 =|. Sache = {2MIn Moi \ atm cosh i) 
SIAN (ee a eee a - ; M+1 
eel m? sinh is(a){;~8M +2M(2M?—1)In 5 i} 
<d(qmcosh.a)= | "Ss 1) Se eae (34) 


Making use of an ultra-relativistic approximation, and our value of f(z), the integral 
in (34) yields the simple result 


(y= | 
and hence (2) =(4/'U,)(radians)#5° 9) 2) eae (36) 


It may be noted that the integrand in (35) is simply the mean square angle of scatter 
as a function of both the primary and secondary energies. 


cosh v 


320m’ exp 3a exp (—7v)(expv—mexpa)d(dmexpi), ...... (35) 


§4. DISCUSSION OF RESULTS 

It is perhaps surprising that it is possible to derive from the total cross section 
in the laboratory frame of reference, not only the total cross section, but also 
something of the angular distribution in the centre-of-mass frame of reference, and 
thus also indirectly in the laboratory system. It is, however, a fact that isotropy 
in the centre-of-mass system is possible only when certain restrictions are imposed 
upon the total cross section in the laboratory system. This is essentially a 
consequence of the circumstance that the cross section in the laboratory frame is 
defined over a much wider range of energies than in the centre-of-mass frame. 
The differential cross section obtained by our method shows an angular 
distribution with a diffuse maximum at right angles to the motion of the two 
incident nucleons in the centre-of-mass system. This characteristic is not a 
property of the particular form of the total cross section F(U), but is shared by all 
other reasonable cross sections which are homogeneous (with respect to the energies 
of the incident and secondary nucleons) and vanish as the momentum goes to zero 
in the laboratory frame. ‘The differential cross section is negligible for scattering 
in the backward and forward directions; this, however, is only true in the ultra- 
relativistic approximation, and it can be seen from formula (27) that at lower 
energies the contribution from forward and backward scattering increases con- 
tinuously. It is thus not inconsistent with low-energy laboratory experiments. 
It is important to notice also that we have excluded completely the consideration 
of elastic scattering, which is predominant at low energies, and probably contri- 
butes significantly to the scattering near 0° and 180° even at very high energies. 
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The neglect of this scattering is not important for our purposes, as it does not 
directly contribute to the lateral spread of air showers. Our method is applicable 
to any kind of inelastic collision between two elementary particles, and is indepen- 
dent of the kind of particle emitted. We have considered the application of the 
theory to plural meson production in nucleon-nucleon collisions, and with the 
particular total cross section assumed in formula (28) their angular distribution is 
the same as for the scattered nucleons. 

The mean square angle of scatter given by (36) is applicable to nucleon—nucleon 
collisions only; what are usually observed experimentally are nucleon—nucleus 
collisions. Consequently, our present results will only permit us to give a qualita- 
tive picture of what is to be expected. It is almost certain that there occurs a 
cascading effect within the nucleus itself, so that the emitted particles may result 
from either the primary, secondary, or even tertiary, collisions. The latter 
particles will have a mean square scattering angle greater than that which we have 
calculated. We shall give the complete theory of this process shortly. 

In the case of the electron—photon cascade, it is well known that the only 
important source of deflection, giving rise to the lateral and angular spread of the 
shower, is the multiple coulomb scattering of the shower electrons by the nuclei of 
the material traversed. ‘This effect gives, in a path length equal to one cascade 
unit, a mean square angular deflection proportional to 1/E?. On this basis the 
theory of the lateral and angular spread of the electrons and photons in cosmic-ray 
showers has been developed. On the other hand we find that in a single nucleon— 
nucleon collision the mean square angle of scattering is inversely proportional to the 
primary energy. We are thus inevitably led to the conclusion that in considering 
electrons and photons on the one hand, and nucleons an the other, the nucleons 
will have a greater lateral and angular spread. Our work above will thus allow us 
to form the following qualitative picture of air showers as far as high-energy particles 
are concerned. ‘The core of the shower will consist of a mixture of all products of 
high-energy nuclear collisions, but naturally the electrons generated by the cascade 
process will predominate. As one passes from the core to the periphery of the 
shower, the proportion of electrons to penetrating particles will fall. This picture 
appears to be in good agreement with what has acually been found experimentally 
(cf. Eidus et al. 1950, Cocconi et al. 1949, 1950). 
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APPENDIX 


We now give a proof of the fact that the energy of an outgoing particle whose 
rest mass does not exceed two proton masses cannot exceed the energy of the 
incident particle in the centre-of-mass system. The proof given below is quite 
general and holds regardless of the number of particles given off. 

Let us denote by m,, P, and U,,, the mass measured in proton mass units, the 
absolute value of the momentum, and the total energy, respectively, of the 7th 
particle in the centre-of-mass system. 


Since U2=P, 2+P,2+P,2+m2,  vsaee (37) 
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we may regard U, as the length of a vector U, in Euclidean four-space, with the 
three P,,, and m, its Cartesian components. The conservation laws tell us that the 
vector-sum of all the U, lies in the direction of the m axis and that the (ordinary) sum 
of their absolute values is prescribed, i.e. 2U,. If one of the U, is longer than 
half this prescribed total length, it must necessarily contribute more than half to 
the m component &,m,. 
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ABSTRACT. The radiation from a hollow-cathode tube containing rhodium was passed 
four times through an atomic beam of rhodium. Absorption was found in eight lines, 
four of which showed hyperfine structure doublets of separation 0-023 cm~!. The splitting 
could be ascribed to the ground state 4d°5s *F5j.. It was concluded with certainty that 
the spin of the rhodium nucleus 1?3Rh has the value J=}. The relative intensities of the 
components show that the magnetic moment is negative. The observed splitting leads 
to an estimated value of ~=—0-10 nuclear magneton. The values of J and wm indicate 
a configuration (5g9/2)®3p,,. for the last seven protons of the rhodium nucleus in the 
Goeppert-Mayer model. 


§1. INTRODUCTION 
EVERAL investigators (Schtiler and Schmidt 1934, Sibaiya 1937, Williams 
S and Granath 1938) have used hollow-cathode tubes to examine lines of 
the arc spectrum of rhodium, but have been unable to resolve any hyperfine 
structure. Sibaiya (1937) observed doublet structures in two of the ultra- 
violet resonance lines, but Williams and Granath (1938) were able to show that 
self-reversal was the cause of the apparent splitting. 

In the present work, an attempt was made at first to resolve hyperfine structures 
in several rhodium lines with the use of a hollow-cathode tube cooled with liquid 
air and operated at very low currents. The line 5599 a. and the ultra-violet 
resonance lines were examined in this way. When this proved unsuccessful, 
an atomic beam of rhodium was produced and the radiation from the hollow- 
cathode tube was passed through this beam before it was made to enter the 
spectrograph and interferometer. ‘The absorption spectrum showed a sufficiently 
reduced Doppler width to allow doublet structures to be observed in four lines. 
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| The lines investigated with the atomic beam were those arising from a 
/ combination of terms of the configuration 4d*5p with the two lowest levels 
4d°5s4F 5), 72 of the ground term quartet (Sommer 1927). 


§2. EXPERIMENTAL METHOD 
(i) The Hollow-cathode Discharge Tube 


In Schiiler’s original design of the cooled hollow-cathode tube and in various 
' modifications which have since been described (Gurewitch 1949, Fowles 1950, 
_ Arroe and Mack 1950), the cathode forms the outer wall of the tube and protrudes 
through the surface of the cooling liquid. The heat conduction along this metal 
tube leads to a great consumption of liquid air and, unless the tube is made very 
| long, is apt to impair the cooling of the cathode cavity. Also, in this design, the 
_ waxed joint between tube and window has to be kept warm by flowing water. 

The use of metal—glass seals avoids these drawbacks and allows a very simple 
design of discharge tube, several modifications of which have been tried out in 
this laboratory. The type shown in Figure 1 is, apart from a minor modification 
introduced later, that which was used in the present and earlier work. 
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|_ Radiation 
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Cathode Cavity 


Aluminium 
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ae 


Figure 1. ‘The hollow-cathode discharge tube. Figure 2. The atomic beam apparatus. 


The cathode was an aluminium block which fitted tightly into the copper 
tube and which had a cavity of 0:5 cm. diameter and 1:5 cm. depth. A sawcut 
was made down the outside of the block to provide a communicating channel - 
from the small space below the block. Two small screw holes were made in the 
upper face of the cathode block so that this could be removed from the copper 
tube by means of threaded rods screwed into the holes. The tight push-fit 
of the cylindrical cathode used in this work gave adequate thermal contact with 
the copper tube. The removable cathode block was later modified for other work 
in the laboratory by substituting a cone fitting for the cylindrical fitting. The 
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consumption of liquid air when the discharge was switched off was entirely | | 
negligible. This type of discharge tube should be especially suitable for use with | 


liquid-hydrogen cooling. 

The discharge tube was connected to a liquid-air trap containing a few grains 
of charcoal, and to a pumping system and reservoirs of rare gases. No provision 
was made for circulating the carrier gas as the static system remained clean over 
very long periods. Two short pieces (about 30 mg.) of pure rhodium wire were 
inserted in the cathode cavity, and the arc lines were excited in an atmosphere 
of neon at a pressure of 1-2 mm. Hg or of argon at a pressure of 0-6 mm. Hg. 
Currents from 3 ma. to 50 ma. were passed through the discharge. In argon 
at a pressure of about 1 mm. Hg the discharge would usually strike at a voltage 
of about 300 v., and when running at a current of 10 ma. the voltage across the 
discharge tube was about 175 v. 


(ii) The Atomic Beam Apparatus 


Rhodium has a melting point of about 1,950° c., and below this temperature 
its vapour pressure is very low. In order to produce an atomic beam of a density 
sufficient to give absorption in the resonance lines, it was necessary to raise the 
temperature to about 2,200° c. This demanded the use of refractory crucibles 
of very high softening point. 

A number of experiments were conducted with crucibles of beryllia and of 
zirconia, but these did not withstand the high temperature. The only refractory 
material which was found to be satisfactory at the temperatures used was thoria.* 
In the final experiments thoria crucibles of 1-4 cm. depth, 5 mm. internal diameter, 
and 7 mm. external diameter were used; a rim projecting outwards about 
0-5 mm. fitted into the upper turns of the heating filament and supported the 
crucible (see Figure 2). The weight of the crucible contributed to the sagging of 
the filament, but it was found difficult to support the crucible independently. 
Attempts were made to provide a stand on which the bottom of the crucible 
could rest; this took the form of a vertical tungsten rod with a flat end, or a wire 
bent over at the top, but in all cases the crucible was punctured: obviously the 
thoria became too soft at high temperatures to withstand mechanical pressure. 
Also the high electrical conductivity of thoria at these high temperatures led 
frequently to short circuits. 

The heating coil was a helix of tungsten wire of 18 turns, inside which the 
crucible fitted. ‘The diameter of the wire was 0-83 mm. and that of the helix 
was about 8 mm. ‘The ends of the coil were spot-welded on to two upright 
tungsten rods. Windings of thin nickel wire assisted in the spot-welds at the 
extremities of the filament wire; in the hotter region near the helix, platinum 
was used. In addition, the filament was tied to the uprights at the right-angle 
bends by means of tungsten wire of 0-2 mm. thickness. In this way the filament 
was prevented from leaning over when the tungsten began to soften at high 
temperatures. 

However, a sagging of the helix under its own weight and that of the crucible 
did take place when the tungsten softened. ‘This gave rise to a separation of the 
turns at the upper part, with a consequent loss of heating efficiency and a short- 
circuiting of the turns at the lower part of the helix. In order that the separation 


* The authors are indebted to Mr. Powell, of Messrs. Johnson, Matthey, for his advice in 
this matter. 
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q of the turns of the helix should be correct at the highest temperatures used 
| (i.e. after sagging), the helix was initially made with uneven spacings of the 
| turns. 

The crucible assembly was surrounded by two coaxial cylindrical radiation 
| shields made of molybdenum foil. ‘The tungsten rod carrying these shields and 
_ the two tungsten rods serving as supports and current leads for the heating coil 
/ were sealed through a Pyrex ground joint. This was fitted by means of Apiezon 
_ wax into the lower end of the Pyrex tube, the details of which can be seen in 
| Figure 2. The Pyrex tube was surrounded by a water jacket. 


At a power input of 2,000 watts (40 amp., 50 v.) the rhodium could be 


: evaporated at a rate sufficient to produce absorption. A charge of 0-5 gm. of 


rhodium lasted for about 30 minutes of evaporation. The pressure in the tube 


as read on the Penning gauge was maintained below about 2 x 10-4 mm. Hg during 


evaporation. At this pressure the ‘shadow’ of the collimating aperture was 
clearly defined by the rhodium deposit. 
The collimation of the atomic beam was determined by the aperture of the 
crucible (5 mm.), the diameter of the glass neck of the outer tube, which was 
2 cm., and the vertical distance, about 10 cm., between the two. According to 
the formula given by Minkowski and Bruck (1935), the half-value width of the 
- lines should be about one-fifteenth of that in an ordinary gas at the temperature 
of the crucible. This gives a value of 0-007 cm“! for the line width to be expected, 
though it must be remembered that the shape of an absorption line due to an 
atomic beam is different from a normal Doppler distribution. The value 
of the collimation was chosen in order to match the Doppler width with the 
instrumental! resolving limit (see the next paragraph). 


(ui) The Optical System 
In the final experiments the radiation from the hollow-cathode tube was 
made to pass four times through the atomic beam: this was necessary when 
absorption in the weaker lines was studied. Figure 3 shows schematically the 
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Figure 3. The optical system. 


optical arrangement. The light emerging from the window W of the discharge 


- tube was deflected into a horizontal direction by the mirror M,. The lens 


L, formed an image of the cathode cavity near the window M, which was coated 
with aluminium. A concave mirror M, of radius of curvature 37-5 cm. was placed 
so that the image of the hollow cathode was at its centre of curvature. After 
the second and third passages through the beam, a second image was formed 


by M, at Mk, slightly displaced laterally from the first. ‘The radiation was then 
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reflected back through the beam to the mirror M, and then to the lens L, and te 
the spectrograph slit. ; 

The numerical aperture of the beam of light passing through the absorption 
chamber was 1/15. This value was chosen in such a way that the full numerical 
aperture at the hollow cathode, namely 1/5-5, was used, and the spectrograph 
(a numerical aperture of 1/30) was filled with light; at the same time, the image 
of the hollow cathode formed on the slit of the spectrograph was magnified: 
five times. 

The spectrograph was a Hilger instrument with quartz—rocksalt—quartz 
achromats of 150 cm. focal length and 6 cm. diameter. With a quartz Cornu: 
prism, the linear dispersion at 3400 a. was about 20 a/mm. 

The high resolution instrument was a Fabry-Perot etalon placed between} 
the collimator and the prism of the spectrograph. The aluminium-coated 
surfaces of the etalon had a reflectivity of about 81° in the near ultra-violet 
which corresponds to a theoretical resolving limit of about 0-008 cm™! with 
a 5 cm. spacer, which was used in most experiments. This value is about the: 
same as the half-value width of the absorption lines caused by the atomic beam. 
A pure isotope !*Hg lamp, supplied by the General Electric Company, proved 
most valuable in adjusting the etalon for parallelism of the plates. When used 
with a rare earth filter, it provides a source which is both more highly 
monochromatic and more intense than any ordinary source, as for example: 
a cadmium lamp. 

Ilford Special Rapid plates were used to photograph the resonance lines j 
of rhodium in the near ultra-violet. 


§3. RESULTS 


in emission, the lines 3434 a., 3692 a. and 3396 a. arising from transitions 
to the ground term 4d°5s‘Fy). of rhodium were found to be slightly broader: 
than the lines 3658 a. and 3528 a. which arise from transitions to the slightly ' 
higher term 4d*5s*F,),, and which showed the expected Doppler width. But 
even at the lowest current used, 3:5 ma., no structure could be resolved. 

With the atomic beam, absorption was observed in eight lines of the are: 
spectrum of rhodium. In four of these a doublet splitting was clearly resolved, 
and in the other four no structure, but only a single absorption line was observed. 
The structures are indicated in the Table. 


‘Transition Line Structure Separation 
4d®5s 4F 9.—4d®5p *Gizyje 3434 a. Doublet —0-022,+0-001; cm-? 
*F gi0 3396 A. Doublet —0-023, cm7 
Das 3692 a. Doublet —0-023,+ 0-002 cm! 
4 oa 3502 a. Doublet Not measured 
4d85s *F,/2—4d®5p *Goje 3700 a. Unresolved 
{Fo H9 3528 A Unresolved 


*D5)2 3658 A. Unresolved 
xGa 3323 A. Unresolved 


The negative sign in the values of the separation indicates that the component 
of lower intensity lay towards longer wavelength. A photometer tracing of the 
atomic beam absorption in the line 3434 a. is shown in Figure 4, with a tracing of | 
the background superimposed. The two components are marked by arrows. 
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No attempt was made to achieve great accuracy in the values of the separation, 


| The doublet splitting of the line 3434 a. was measured with a travelling eieconcene 
on four plates, that of the line 3692 a. on two plates, and that of the line 3396 a. on 
only one plate. 

The emission line was excited at a current just below that at which self- 
reversal took place in the hollow cathode. The operating current was usually 
~ about 30 ma. in neon at a pressure of 2:0 mm. Hg, or about 17 ma. in argon at a 
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Figure 4. 


pressure of 0-3 mm. Hg. At these currents the line was broadened and flattened 
_ by self-absorption, but not reversed. The exposure time ranged from about 
5 minutes for the strongest line 3434., with only two passages through the atomic 
beam, to about 15 minutes for the weaker lines, with four passages through the 
beam. The emission line alone was always photographed on the same plate 
under the same conditions of the discharge as in the exposure with the atomic 
beam. 
§4. DISCUSSION OF THE RESULTS 

Rhodium consists of only one stable isotope of atomic number 45 and mass 
number 103. Hence the observed doublet splittings must be ascribed to hyperfine 
structure caused by the nuclear spin. In transitions between terms whose J values 
are all greater than 5/2, a doublet splitting can only be caused by a nuclear spin 
of J=}. Any other possible value of J would have caused a structure of four or 
more components. 

Since there are no penetrating electrons in the configuration 4d°5p, the 
splitting of the terms of this configuration can be expected to be small compared 
with that of the term 4d°5s4F,,.._ The validity of this assumption is borne out 
by the fact that all four lines arising from the combination of the ground term with 
terms of the 4d°5p configuration show the same splitting, within the accuracy 
of the measurements. To this degree of accuracy, the-splitting can therefore be 
ascribed entirely to the ground term 4d85s4F,),, and is taken to be about 0-023 cm. 


73-2 
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The fact that the component of longer wavelength is the weaker one shows ~ 
that the hyperfine structure of the term 4d°5s4F oj. is inverted. In this term the | 
spin of the s electron points in the same direction as the total angular momentum; 
the inverted order of the hyperfine structure terms therefore shows that the 
magnetic moment of the nucleus is negative. The calculated value of the intensity 
ratio of the two components is 11/9 for all the doublets. This is consistent with 
the appearance of the photometer tracings, but no actual measurements of 
intensities were made. 

To a first approximation, pure Russell—-Saunders coupling can be assumed 
and the interaction of the nucleus with the d electrons can be neglected. The 
ratio of the doublet separations of the two terms 4Fy,. and ‘F,). arising from the. 
configuration 4d85s is found to be 7:4. This explains the fact that no structure 
was observed in the four lines connected with the term *F 7). 

On the same assumptions, the value of the nuclear magnetic moment can be 
calculated by the usual approximation methods (see, for example, Kopfermann 
1940). The value found was »=-—0-10; nuclear magneton.. An attempt to . 
estimate the effect of the d electrons on’ this value of » was made by assuming 
that the relation between the interval factor A(*F,).) and the splitting factors 
a, and aq, given by Crawford (1935) for the 4d#5s4F 5). term of lanthanum, applied 
also to the 4d°5s4F,,. term of rhodium. The result of including the influcess 
of the d electrons was a reduction of the numerical value of ju by about 10°%. This 
is of the same order of magnitude as the other uncertainties in this approxttta 
calculation of ». The final value » ~ —0-10 n.m. may be taken to be accurate 
within about 25%. 

The comparative ease with which absorption was observed in lines arising 
from transitions from the state 4d°5s4F,). was due to the very high temperature 
of the atomic beam oven. At a temperature of 2,500° k. the number of atoms 
thermally excited to the state *F,,. which is raised by 1,530 cm above the ground 
state *Fy). is one-third of the number in the ground state. ‘This is an interesting 
case of some advantage being derived from the use of a high oven temperature, 
which from other points of view (Doppler width, ease of operation) is regarded 
as a drawback. 

From the mason of nuclear shell structure, Goeppert-Mayer (1950) predicts 
a spin of either } or 9/2 for the rhodium nucleus. Since the experimental value _ 
is } the configuration (529/2)®3py)2 for the last seven protons is suggested. This — 
is toned by the yalue of the nuclear magnetic moment — 0-10 n.m. which lies 
vey clese to the line J=/—4 in the Schmidt diagram (Schmidt 1937). Since 
I=14, I must have the value 1, and the ground state of the nucleus is a 2P,, 1/2 State. 
This assignment is in disagreement with that quoted by Feenberg and Hammack 
(1949) as being most preferred, namely 


(1s)*(2p)*(3d)i(4f)41 (Sg)2(2s)* 
re seal 
in which the orbits (5g)!® in square brackets refer to the ground state and 
(5g'*(2s)! to a metastable state. If the order of these terms is inverted so that 
(5g)!*(2s)! describes the ground state, the spin } is correctly given, but the small: 
value of . cannot be reconciled with the interpretation of the ground staté as an 
Sj), term. It must therefore be concluded that at Z=45 the:3p orbit is’ not 
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j pushed upwards as much as is predicted by Feenberg and Hammack, and that the 
4 scheme of Goeppert-Mayer gives the better representation of the rhodium 
nucleus. 
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ABSTRACT. Consideration is made of the predictions of the Bloch spin wave model 
within its region of approximate validity, for the magnetic scattering of very slow neutrons. 
by a ferromagnetic. It is found that observable peaks of intensity in the angular distribution 
of scattered neutrons occur due to spin wave absorption. Neutrons are also scattered 
with phonon absorption and by an analogous calculation it can be shown that similar 
peaks occur due to this vibrational scattering. A possible method is suggested of 
differentiating between the two kinds of scattering in order to provide a check of the Bloch 


model. 


Sst NMR OD Ci lON 

LASTIC magnetic scattering of slow neutrons by ferromagnetics has been: 
the object of much theoretical study (see, for example, Halpern and 
Johnson 1939). In fact, in most of the work that has been done it is. 
implicit that the elementary magnetic moments of the terromagnetic are not 
changed in direction by the scattering. This is broadly true at low temperatures. 
and for the slow neutrons used, as the neutrons will not be able to supply the 
exchange energy that the ferromagnetic would gain if an atomic magnetic 

moment turned out of the direction of magnetization of its domain. 

But at higher temperatures, say, approaching half the Curie temperature, 
measured in absolute units, there will be in each domain a fair number of 
disoriented magnetic moments. ‘Thus one would expect an effect in which the 
scatterer gave energy to the neutrons by means of a turn round of atomic magnetic. 
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moment into the direction of domain magnetization. This would result in a loss | 
of exchange energy by the scatterer. The cross section for this effect should 
increase with temperature and also with the wavelength of the incoming neutrons. 
If then the incoming neutrons had a wavelength greater than the Bragg reflection 
wavelength there could be no elastic magnetic scattering and energetic 
considerations would make an energy transfer from the neutron to the scatterer 
extremely unlikely. Thus, any magnetic scattering observed would ‘be from the 
inelastic scattering in which the neutrons gain energy. 

The main object of the work presented here is to investigate the order of 
magnitude and distinguishing features of this effect. ‘The difficulty arises that 
the problem of ferromagnetism is by no means solved, and thus.some doubt is 
cast on the validity of the wave functions that are used for the initial and final 
states of the scatterer. 

At present the only wave functions that can be used for scattering calculations 
and that take into account the co-operative nature of ferromagnetism are those of the 
Bloch (1930) spin wave model. An objection to the spin wave picture is that while 
the wave functions give a good approximation to the solution of the Schrodinger 
equation at low temperatures, the problem of finding a wave function to give 
a good solution near the Curie point has not been overcome. However, for 
temperatures up to about a third of the Curie point the wave functions of the 
spin wave picture should be reasonably good. 

Thus though the Bloch spin wave model is imperfect, it seems to be the 
only means, at present, of approaching the problem of inelastic magnetic 
scattering. Comparison of the results with experiment to investigate the 
predictions of the Bloch model are a useful further check on its validity. 

Perturbation theory, and a simplified model of one electron per atom, are used. 


§2. CROSS SECTION FOR MAGNETIC SCATTERING 


After Heisenberg had demonstrated that the exchange interaction between 
atomic electrons provides a suitable mechanism for ferromagnetism, Bloch (1930) 
was able to give a simple model. Consider a crystal lattice containing one 
valence s-electron in each atom and with these electrons responsible for the 
ferromagnetism. ‘Then an approximate eigenfunction of the total Hamiltonian, 
including the exchange interaction between nearest neighbours, for just one 


disoriented electron spin in a domain containing N ferromagnetic electrons 
is (Sommerfeld and Bethe 1933) 


P= NM(N!)P PX + Q)X exp (AR) F4(74)Sr/2(S1) 
Q 1 


. wg EP) SRG) oe eae) Sunn) ee (2.1) 


k is a vector of the form k=(27/a)(m/Ny, m2/N2, n3/N3) where a is the side 
of the unit cell, assumed to be cubic, Ny, N,, N, are large integers with 
N,a, N,a, Na the sides of the domain of the crystal considered and nj, m3, ng is 
an integral vector. (+Q) is the operator of permutation of both the electron 
space and spin coordinates (r, s), according as the permutation is even or odd. 
R, is the position vector of the /th atom and runs over all the atoms, the spin of 
the electron on the /th atom being the one that is reversed. ‘fH’; is the space wave 
function of the electron on the /th atom and S the spin wave function, 
3/2 denoting an electron spin oriented in one direction and S_,,. an orientation 
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the opposite direction. The electron wave functions WV), ¥, ..., y are 
qual in shape, but are centred round different origins, namely, the nuclei of the 
arious atoms of the crystal. N-"(N !)-1? is a normalizing factor. 

The normalized wave function for no disoriented spins in the domain is 


Q'= (MIP ECE QE) Suialse) oe PN)Sip(s) . . » Py Sipe(sy)- 


Taking s-electrons instead of d-electrons as carriers of the ferromagnetism 
is not significant. For the orbital current of the d-electrons in a crystalline 
7 ferromagnetic is quenched and the ferromagnetism is wholly due to the spin. 
i) Scattering by a single crystal is investigated where the quantum mechanical 
: scattering unit is the magnetic domain. As the scattering is inelastic the 
contribution from the different scattering units will superpose incoherently. 
_ In a magnetized crystal the scattering unit is the whole crystal. 

& Suppose now that a neutron beam, taken unpolarized for simplicity and with 
_ the neutrons having wavelength greater than the Bragg reflection wavelength, 
_ is inelastically scattered, the scatterer going from state ® to state ®’. Thus one 
| spin wave of wave vector k is annihilated. The scatterer is the ferromagnetic 
_ domain and perturbation theory is used. Then the scattered wave of the neutron 
| and the domain for the neutron at a large distance 7 from the scatterer is (Halpern 
-and Johnson 1939), with unit density of incident neutrons, 


i 2m 
e T= art exp (ip’r)x,0'(’, x.,p’|V|®, x. Pp). ---+-- (2.3) 


_ Primes indicate final states. p is the wave number of the neutron, y, the spin 
| _ state of the neutron, Mis the neutron mass, (h?/2 M)(p? — p’*) = E’ — E, Ebeing the 
energy of the scatterer. For a body-centred cubic lattice, for instance 
- (Sommerfeld and Bethe 1933, eqn. 6.20) 


é E=E,+E(k), E(k)=8J(1— cos 4k,a cos 4k, cos 3k,a), ...... (2.4) 


J being the positive exchange integral, and it is found for any crystal lattice, 
when |k|a is small, that : 
I ATS Vee NI paella Ee Ra EEE (2.5) 


| A 
EEG UPN OE 
so that a (P iD?) A Sele oe ae (2.6) 


V is the energy of interaction between the neutron andthe ferromagnetic electrons 
(Halpern and Johnson 1939) 


r is the position vector of the neutron, r, that of the /th electron, a, is the Dirac 
matrix vector for the /th electron and p is the magnetic moment of the neutron 


CHM C) VS = ee ke ov (2.8) 


y being the ratio of the neutron magnetic moment to the nuclear magneton. 


(P' xe1 0 Box ®) = [Pel Sie ¥ eS COM ep dtd 


[r—P 


(xe lau] x,) ex fip—p')e EE exp (R(k—q)YINW® (2.9) 
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The origin of the vectors r and r, is the centre of atom / (the choice of atom /is, 
of course, quite arbitrary and immaterial); q=p’—p. Equation (2.9) can be: 
evaluated by the method given by Halpern and Johnson (1939). It is then | 
found that the differential cross section per atom for scattering of the neutron | 
into the angle dw,, with absorption of the spin wave k is 


Qey\2 1 1p’ F 
(=) Wale p oe S selielerets (2.10) 
where F=[¥*(7,)¥(%) exp (—iq . r,)4?7, as (2:15 
S=1+¢,2, e=(p—p'/lp—p'|=-allq|  -.-+-. (2.12) 


and e,, is the component of e in the direction of magnetization of the domain. 
dw, is the differential scattering (solid) angle. 


D=|Bexp (iR(k~4)} 


DS x exp {7R,(k —q) } |? 2 {exp 7k,(k—q)}[?. ...... (2.13) 


In the first factor on the right-hand side of (2.13) summation is carried out over . 
all the unit cells in the crystal and in the second over all atoms in any one 
elementary cell; therefore 

_ sin? (k—q),3Nja sin? (k— q),3N,a sin? (k— q),3.N34 : - 
Dank (k 4) io sin Eada (kaye, 41 ee 


Up to now the scatterer has been taken to have only one disoriented electronic 
spin, that is only one spin wave. On the assumption that there are few spin waves 
and that they do not interact, the differential cross section per atom for 
annihilation of a spin wave k, if there are initially ZL, spin waves with wave 


vector k is 24,\ 2 f 
(= y ) eee diay DF*S. 


N? 16 
This corresponds to equation (2.10). The approximation of non-interacting 
spin waves is probably adequate when the proportion of disoriented spins is 
small. And for temperatures up to about one-third of the Curie temperature 
the number of wrong spins is small. 

In this approximation equations (2.4) and (2.6) still hold. 


mc2 


SS CAINE RUNG BYa TAT SING I CRiyEs uPA 
For variations of k or p’ large compared with 27a/N,, we may substitute 
for D the expression 
(27/a)?N,N,N,N?X8(k + p— p’ —27n/a). eee (3.1) 


N, is the number of atoms in any elementary cell. Thus N,N,N,N,=N. The 
summation over n is over all integral vectors including zero for which 
x,exp (277R,n/a) 40; that is 27n/a is in the reciprocal lattice. We shall consider 
the wavelength of the neutron as greater than the Bragg reflection wavelength. 
It can be shown that for this case scattering with excitation of a spin wave, and 
the neutron losing energy, does not take place into appreciable angles, for 
p <p, and k (given by (2.4) with a change of sign on the right-hand side) cannot 
be large enough for k+ p—p’—2zn/a to be zero forn0. The case of elastic 


scattering gives rise to a function &,5(p—p’—2z7n/a) and there is no elastic 
angular scattering for long wavelengths. 


Slow Neutron Scattering by Ferromagnetic Crystals I1OI 


The total cross section, for the initial neutron wavelength greater than the 
Bragg reflection wavelength, is, ignoring the small angle contribution from n =0, 


a Wel ee Caen, (3.2) 
n+0 
| 2ey\? rp! ae 
where o,= (=) | dk [Feo yFX(p ~p') 


d(k+p—p’—2nmija)a(lte?)Ly eee (3.3) 
This follows from equation (2.15) by summing over k and. replacing the 
summation by an integration: 
: ‘N,N2N3a° 1, 
espa d kk. 
L,, is the number of spin waves in the scatterer having wave vector k. 
The spin waves obey Bose statistics, so the cross section for scattering at an 
absolute temperature T is given by substituting in (3.3) 
=n (exp ye) Pel tg ec, (3.4) 
where k is Boltzmann’s constant, and E(k) is defined in (2.4). Suppose that 
ka is small enough for (2.5) to hold. Then from (2.6) and (3.3) the following two 
relations are fulfilled 
pr =p tak? oo. (3.5) PD 0s cr ce (3.6) 
where «=2J Ma?/h?>1, p,=p—2an/a. Then (3.5) and (3.6) give, by 


elimination of k, 


(a—1)p*—Zp pacos@+(a+1)p?—Opr=0 ....5. (327) 
where ¢ is the angle between p’ and p, and where 
2n\2 4 da\? , 4 
229=(=) ee dP (=) n? npr, nog oat (3.8) 
2 2 9 we as 1/2 
whence p =p, eee NET EOS 7 a Se tetas (3.9) 
CC 


Both the positive and negative signs in the numerator give equally valid 
solutions. For, take p, as polar axis and take the usual coordinates p’, 0, ¢ for 
the end point of the vector p’. To find the possible values of p’, satisfying 
(3.5) and (3.6), that lie in any plane ¢ =constant (which, of course, contains the 
vectors p’, p, and k), first take a z axis perpendicular to the plane, as shown in the 
Figure. Then construct the surfaces of revolution z =p” and z=p"+a(p’—p,). 
Project the line of intersection of the two surfaces perpendicularly on to the 
plane ¢=constant. The end point of p’ then moves on this projection in the 
plane ¢=constant. As every plane $=constant is equivalent the end point of 
p’ moves on the surface of revolution given by rotating the projection round the 
polar axis. This shows how to each value of 0, ¢ correspond two values of p’. 
From (3.5) and (3.9) 

kh? = p?(2a cos? 6—2%+2+(%-1)0+2cosé {a? cos? 6—« 
HEE (OTOP aaa) waa (3.10) 
For the sake of definiteness, take now the scatterer as being an unmagnetized 
crystal of a body-centred cubic lattice. The wave number corresponding to the 


Bragg reflection wavelength is 27/a,/2. Thus if p=27«/a, «<1/\/2. Taking 
the vectors of the reciprocal lattice as axes, consider as examples of the two cases 


p=(2z7/a) (x, 0,0), p=(27/a) (k/+/3, k/4/3, k/+/3). 
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In Table 1 are given the values of p,? and of Q for the lowest values of n. As 
will be seen later the higher values are not so important. In the body-centred 
lattice only those values of n occur, the sum of whose components is even (see 
equation (2.14)). Numbers are given for the particular case x = 4. 

It is seen that Q<1. Though it is possible for Q to be zero when « has its 
limiting value (1/./2 in the body-centred case), for smaller values of « it is seen 
from Table 1 that Q~1. Unless otherwise stated the values of Q in what 
follows are such that Ow~1. . 

The limits of @ are, from (3.10), 


Ojamf1 (o, SLO uo2e sin? 0. ae (3.11) 

also Bip2=Oje Odie | Seen (3.12) 
LI MG ds kO 3 

Now RRO (W/2Ma)’ OOO OIC (3.13) 


where © is the Curie temperature. For a body-centred cubic J/R© =1/6-1 
(Sommerfeld and Bethe 1933). Thus, taking the constants and body-centred 
lattice structure of iron, a=2°85 a., © =1,043 0°K., «=62. From (3.12) 
| ka? =(O/a)p,2a? =(40QO/a)(ap,/27)?. eae (3.14) 
Thus the expansion of E(k) (equations (2.4), (2.5) ) is certainly valid for 
(ap,/27)? <1, for a body-centred cubic. 


0 9 P p’ 


a, paraboloid z=p” ; 6b, paraboloid z=p?-+-«(p’—p,)? ; c, intersection of A and B ; d, projection 
ae = 
of C ; e, cone for phonon absorption. (P,K)?=p? OP,=p, OP’=p’. 


To evaluate the cross section (3.3) the values of L, are needed. For ka small 
enough it is found, using in turn (2.5), (3.12), (3.13) 


EXE)/RT = (apij/27)-O(110/2)s,. soe See eee (3.15) 

‘It is thus evident that at room temperatures or higher and for (ap,/27)?<1 a 

good approximation for L, is given by 
T2MRT Ree 1 

«"FepeQ ~ T0O@p ae > 

In Table 2 values of k°a? and o, are given, worked out on the assumption that 

(2.3) and\(3.16) are valid. As the crystal is assumed to be unmagnetized e,,2 is 

replaced by its average, }. 

This assumption is made only to simplify the formulae. In fact, from the 

point of view of theoretical prediction it is rather better to havea magnetized 


Slow Neutron Scattering by Ferromagnetic Crystals L103 


crystal, for then question arises, such as for instance, the smallness of domain 
size affecting the theory. For a magnetized crystal the cross sections o, given 
_ below must be multiplied by 3(1+7,,2), where 7,, is the component of n in the 
direction of magnetization. 


Using (3.3), (3.9), (3.11), (3.16) and (3.17) it follows that 
Zery\eres pi 27 cos 8 i 1 

=) i =|) ae = Fe Oenla Kk); 

= ae) pe es aes vn :)(s5) Oana ee 


4 The factor 2 in front comes from the fact that to each value of cos @ in the range 
| (3.11) correspond two values of p’. 


m 


|Z 
> Table 1 
eB p=(27/a)(x, 0, 0) p=(27/a\(k, x, «)/ V3 
i 2 2 
_ n Pi =P, O 
} (4) Q (2) ss 
i K=} K=} Keak K=3 
\ 2 : : ; 
(0, 0, 0) K 7 0) 0 K2 4 0 0) 
(4, 1, 0) 2—deee 5 — : 2—2-3e-+u? 11 —— 0-77 
i 5 | 
2—2k a » 
o- —1, 0) 2—2«+k? yi Sie 5 2K 2°25 Dane 0°88 
13 24-20 12 2 
Bei ere 2H. G 2+k? 225 Se 0-88 
13 Pg) 12 242:3 
eS | oe gutter Gt | coe 0 
Z ar site “AE 
L 5 ee 4 92293 
(1, 0, 1) DOE ee ee ee eH 
ne. ae calie . : a ONS FL 
’ 5 pees) 4 2 
Pee tat le 2+? 2-25 Ss 0-88 
fF i} 
13 +2 12 2 
(—1 0 1) 2- 2e+K? ws ny, ayaa 3 2-2 2:25 ae 0:88 
13 242 iW) Sane) 
Be OD 2b 2eret be ee Teele bese = 834 a 0-92 
ie Stee OC etek 
9 2 358 
0, 1, 1) Dol js mae : 2—-2:3«+K? = 14 a 0:77 
pee 
9 7) 8 5 
m0, 1, —1) 2+K , ae : +h 2-25 = 0-88 
ss ee 
9 2 
i, —1, 1) 2+K2 ; i 2 - : 2+ 2-25 oa 0-88 
sigh S51K 
9 2 8 5 2+2:3k 
— — 2 as = 3+ Ke? 2. : : 
«0, —1, —1) 2+k rm Te 9 24+2:3K+K 3-4 523ene 0-92 


Equation (3.17) is not an exact formula for o,, but within the limitations 

_ given is a good approximation. It can be seen from the column C in Table 2 how 
good the approximation for L, is, and from column A how good the 
approximation Jka? for the energy change is. It is seen, for example, that for 
=(27/a)(4, 0,0) the o, given is a very good approximation for n=(1, 1, 0) 
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order of magnitude. 
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—1),while for the other values of n tabulated it holds in 


From column A it is evident that F?(2mn/a—k) can be replaced by 


F?(2mn/a) since |k|<|2zn/a|. 


Various estimates of F have been made by 


various authors in connection with elastic magnetic scattering of neutrons 


by iron. 


If then we take, consistently with the previous assumptions for the 


lattice constant and the Curie temperature, the form factor to be the same as 
that for the 3d electrons of iron, we can make use of the previous estimates. 
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Hamermesh (1942) calculated F using Hartree wave functions for free iron 


atoms. 


Steinberger and Wick (1949) took into account the Fock exchange 
integrals of a 3d electron with all other shells and with its own. 
F?(2nn/a) are given in Table 3. 
n’=7°. 


Values of 


N(n?) is the number of values of n such that 


= 
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It is evident from Table 3 that %,,:..0, is small. For the increase in N(n?) is 


more than compensated by the decrease in F?, while the value of any particular 


6, decreases rapidly with increasing n due to ahe factor 1/exp{E(k)/RT}—1 in 
the integrand. ‘Taking the higher, and probably better, value of Steinberger 
and Wick for F?, it is seen from Table 2 and equation (3.17) that the cross section 
for scattering into a solid angle of 7Q/« radians around the direction, approxi- 
mately, of (p—27n/a) for p=(27/a)(3, 0, 0) and n=(1, 1, 0), and T=400°K. is 
about 1/50. A value of 62 has been used for «. Since 7Q/« is about 1/25, this 
corresponds to a differential cross section of approximately 0-5 barn per steradian. 
Taking, as usual, the structure for iron it is seen that the mean free path for this 
n-type scattering (with n=(1, 1,0), p=(27/a)(4, 0,0), T=400°K.) is about 
580 cm. 


Table 3 
n? y 4 6 8 10 iD 
N(n®) 12 6 24 12 24 6 
F?(27n/a) H. 0:23 0:078 0:036 0-014 0-0064 0:0028 
S. & W. 0-36 0:2 0-1 0:048 0:025 0:014 


H,.=Hamermesh ; S. & W.=Steinberger and Wick. 


It is in the high-temperature range that the assumption of non-interaction 
of spin waves begins to part with reality, so little confidence can be placed on the 
predictions of (3.17) for T above half the Curie temperature. 

The calculations of this section have been carried out for a crystal of ‘ ideal 
iron’, which has the lattice structure and Curie temperature of iron, but only 
one ferromagnetic electron per atom. ‘The construction of a model with two 
spins per atom, which approximates closely to real iron (2.2 spins per atom), 


~ is more complicated, though it presents no fundamental difficulties. In such a 


theory the general features of the present model would be preserved. The 
peculiar angular distribution of the scattering depends on the energy and 
momentum equations—equations (3.5) and (3.6). In a two spin per atom model 
these equations would still hold, but theré would be two possible equations 
corresponding to (3.5). Consider the spin of the two electrons, which wil! be 
coupled (Mott and Jones 1936) as quantized in the direction of magnetization. 
Then ‘one equation (3.5), with a constant «,, will correspond to electron spin 
transitions —1->+0. The other, with constant «,, will correspond to the electron 
=Spp transition 0-1. Both «, and «, will be somewhat smaller than «, because 
ey Corgespond to one energy change while the Curie temperature corresponds 
to two. Thus with the two spin model there will be the two kinds of scattering 
(that is, with «, and a, respectively) superimposed. Decrease in « leads to an 
increase in the angular range for each part o(m) of the scattering (see equation 
(3.18) ), while leaving the‘differential cross section unaltered. ‘Thus with the two 
spin per atom model an increase in the angular spread round the directions 
p —2mn/a, and an increase in the total cross section, is to be expected. 

For neutrons of wavelength greater than the Bragg reflection wavelength 
elastic nuclear scattering will be small and isotropically distributed. It is due to 
presence of isotopes and nuclei with spin and is less than 0-25 barn per steradian 
(Cassels and Latham 1948). ‘This small isotropic background of elastic scattering 
provides no obstacle to the observation of spin wave peaks. 

However, the neutron beam undergoes thermo-vibrational scattering 
(Weinstock 1944), both through magnetic and through nuclear forces. In fact, 
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the very slow neutrons considered will be scattered with absorption of phonons | 


in the same way as they are scattered with absorption of spin waves. 

In order to gain some insight into this effect, the Debye approximation— _ 
that the energy of a phonon of wave vector k is fic, |k| where c, is the constant 
average sound velocity—-is used. In this approximation the difference between 
transverse and longitudinal sound velocities is neglected for simplicity. The 
approximation of a sound velocity. independent of wavelength holds for 
he, |k| <}©) where ©, is the Debye temperature (or |k|<7/ainiron). The effect 
of phonon absorption can be found by a calculation similar to that just performed, 

with (2.Mc,/h) |k| put for «Rk? on the right-hand side of equation (3.5). For 
example, the only change in finding the scattering angle is that the paraboloid 
b in the Figure is replaced by the cone e. ‘Thus the maxima that occur due to 
phonon absorption are centred in the same directions as these due to spin wave 
absorption. 

For the maximum corresponding to p=(27/a)(4, 0, 0), n=(1, 1, 0) phonon 
absorption gives a differential-cross section of 0-5 barn per steradian (with nuclear 
vibrational scattering dominating magnetic vibrational scattering). The angular 
range of the peak is of the same order as the angular range of the corresponding 
spin wave absorption peak on the one spin per atom model. The Debye 
approximation holds for the low frequency phonons absorbed to make this 
peak. But for p=(27/a)(4, 0, 0) n=(—1, 1, 0) OP, (see Figure) is greater than 
in the preceding example, and thus the cone e (phonon absorption) cuts wu very 
much outside the intersection of b witha. ‘This shows, if the Debye approximation 
holds, that the phonon peak has a very much greater angular spread than the 
corresponding spin wave peak. In fact, the Debye approximation is probably 
beginning to break down for the high-frequency phonons absorbed to make this 
peak. But it can be qualitatively seen that the peak is much broader and less 
intense than the corresponding spin wave peak (Table 2). 

An increase in the differential cross section of 0-5 barn per steradian, such as 
is predicted for some spin wave peaks, should be well observable with the neutron 

fluxes available from piles. Another possibility of observing the effect comes 
from the fact that (equations (3.12), (3.5)) the energy of the scattered neutrons 
is approximately (h?/2M)(p,?+p*). For p=(27/a)(%, 0,0) n=(1, 1, 0) this is 
about five times the energy of the incoming neutrons. The different variation 
with n of spin wave and phonon peaks (as outlined in the last paragraph) might 
well serve to distinguish the two. The magnetic properties of spin wave peaks 
also offer a means of distinction. 

For a magnetized crystal and unpolarized neutrons the cross section of 
Table 2 should be multiplied by 3(1+7,7). This factor varies between 
* and 3, according to the relative orientation of n and the direction of 
magnetization. Ifthe Bloch spin wave model were correct, the variation of cross 
section with the direction of magnetization of the crystal could presumably be 
observed. ‘This variation with magnetization is characteristic of spin wave 
absorption, and thus may provide a test of the Bloch model. (A, different, 
variation occurs also for the magnetic part of the scattering with phonon 
absorption, but calculation shows this scattering to be almost certainly negligible.) 

The above analysis is for monochromatic neutrons, but if the direction of the 
incoming neutrons, that is the direction of p, is the same as the direction n, 
a peak will be observed in the direction p—2zn/a even if the neutrons are not 
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monochromatic. For such an experiment, of course, the neutrons would all 
have wavelength greater than the Bragg reflection wavelength, both in order that 
this inelastic magnetic scattering should be of considerable magnitude, and that 
elastic scattering should be small and isotropic. 


§4. CONCLUSION 

It has been shown that, assuming the Bloch spin wave model of a ferro- 
magnetic, considerable inelastic magnetic scattering of very slow neutrons in 
which the neutrons gain energy takes place. The essential reason for this is that 
for such slow neutrons the temperature of the scatterer is considerably greater 
than that of the neutrons. Though the total cross section for scattering of this 
kind is not large, the scattering has a peculiarly limited angular distribution 
(provided the scatterer is a single crystal and the neutrons are monochromatic), 
due to the energy momentum relations, and in particular the small ‘ effective 
mass’ of a spin wave. This limitation results in the differential cross section 
being large for those angles into which scattering takes place. Considerable 
nuclear scattering with phonon absorption occurs, and gives peaks centred in 
the same directions as, though not generally coincident with, the spin wave 
peaks. In a magnetized crystal the variation of spin wave absorption scattering 
with direction of magnetization might be used to distinguish it from vibrational 
scattering. 

For two reasons the effect of spin wave and phonon absorption is most 
important for neutrons of wavelength somewhat greater than the Bragg reflection 
wavelength. Firstly for neutrons of wavelength less than the Bragg reflection 
wavelength it can be seen from the mathematical methods of §3, and it is indeed 
obvious from energy considerations, that this type of scattering is less important. 
Secondly for wavelengths greater than the Bragg reflection wavelength elastic 
nuclear scattering is small and isotropic. ‘Thus the peaks of intensity discussed 
here appear more clearly. 

The calculations have been carried out here for iron, but similar, though not 
so marked, peaks are to be expected for nickel. 
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Antiferromagnetism by the Bethe Method 
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ABSTRACT. 'The Bethe method for order—disorder is applied to a lattice of magnetic 
dipoles with antiferromagnetic nearest-neighbour interactions. Magnetization and entropy 
are calculated for all values of temperature and magnetic field, and an anti-parallel ordered 
phase is found below a definite transition curve. The computed curves agree qualitatively 
with those observed experimentally with certain paramagnetic salts at low temperatures. 


OME recent experiments (Garrett 1951, Kurti 1951, Gorter 1951) on 
paramagnetic salts at low temperatures have yielded results which may 
be interpreted according to a theory proposed by Sauer and Temperley 

(1940). The anomalies in the entropy and magnetic susceptibility are considered 
to be due to anti-parallel ordering of the lattice of atomic spins, due to magnetic 
dipole-dipole interactions. ‘The crystal behaves as if it were antiferromagnetic, 
the anti-parallel state existing when the temperature and external magnetic field 
have values lying below a certain curve in the H,T plane. 

The theory of Sauer and Temperley uses the Bragg—Williams (1934) 
approximation for the energy. It is well known that better approximations 
are available to solve order—disorder problems, in particular the method of Bethe 
(1935), which takes rather more account of local order. Unfortunately, the 
Bethe method assumes that the ordering of the lattice is due to the interaction of 
any atom with one or more shells of neighbouring atoms. In the case of 
dipole-dipole forces this assumption is rather artificial, since the interaction is 
highly anisotropic, and falls off only as the inverse cube of the distance between 
the spins, so that a very large number of shells, which would interact among 
themselves in a complicated way, would be needed. Nevertheless, it is worth 
considering a very simple model, to see how far it can give results in agreement 
with experiment. 

Let two lattices interpenetrate in such a way that each atom of one lattice 
has z nearest neighbours on the other. At each lattice point let there be an 
Ising dipole, i.e. a spin variable S which can only take the values +1 and —1. 
Ifz andj refer to neighbouring atoms, there is to be an interaction energy Ju?S,S,,, 
where J is positive, S; is to be the magnetic moment of the 7th atom, measured 
along a direction in which there is an external magnetic field Hy. The con- 
figuration of lowest energy for such a system is one in which all spins on one 
lattice point to the right and all on the other lattice point to the left. 

This problem has been solved formally by the Bethe method by Firgau (1942). 
We follow her procedure, but in the more elegant formulation suggested by Weiss 
(1948). Suppose, as in Figure J (a), that wS, is the moment of a central atom 
ona site where the dipole points to the right in the lowest-energy state (right-hand 
site), that is, Sy takes the value +1 when this dipole points to the right. wS, 
is the total moment of the shell of z atoms surrounding S,; S, takes the values 
z,2—2,2—4,.... —z, and its sign is to be positive when all the dipoles of the 
shell point to the left, as in the diagram, i.e. S, = —XS;, where the sum is over 
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sites surrounding Sp. Again, Hy is the field acting on S,; as usual in the Bethe 
_ method, this must be the external field, since the internal field acting on the 
| central atom is taken account of by considering explicitly the interaction between 
| Sy, and the shell S,. However, H, is a hypothetical internal field acting on S, 
alone, and represents the combined effects of the external field and of the 
| neighbours of the shell (apart from the central atom). The positive senses of 
_ these fields are indicated in Figure 1 (a). 
The partition function for this cluster may be written 
3 g! 
geen |e CEMA CEA 
— pi (— BH So — pH, S, — WI SoS,)/RT} 


ee ey) 
=X Pa («+ =) oe (on) ete (1) 
where x) =exp(—pH)/RT), x =exp(—wH,/RT), y=exp(—p2J/RT. 


3 <—-—S <— I, Ss, —> ss, —— 
b. 4—>  &—— <—s <— 
z (a) (2) 
We Figure 1 (a). Right-hand site. (6) Left-hand site. 
\pee 
Z For atoms on, and surrounding, a left-hand site we define the fields and 
. _ spins according to the scheme of Figure 1 (6), and obtain similarly a partition 
function : , gy by ae 
fy =f) (24+5 + 
; ee (5 a) + z( ey 
_ where ¢)=exp(—pH)'/RT, € =exp(—puH,'/RT. 
: The average moment of each type of atom may now be determined; thus mp, 
_ the average value of So, refers to an atom on a right-hand site, and is given by 


a 0 
iy = — 1% =~ [log (P-£)] 


£42) —2(ay+ AV] /[n(2+2) +4 (9+ 2) 
lS) sle* a) IG) (ora) | 


and similarly m, =(wS,/z),,= — — oes — ? Tog (p,f2) | ete: 


But the two lattices are ues and are subject to the same external field. 
mother words, H,-= — H,; te. £&,=1/xp. 

Moreover, since any two right-hand sites or any two left-hand sites are 
equivalent, whether considered as occupied by a central atom or yr by an atom in 
a shell, they must possess the same average moment. That is ,My = my! and my! =m. 

We now have sufficient equations to eliminate the unknown fields H,, H,’. The 
partition function for the whole crystal can be written (P.F.) = [(p.£.)(p.£.)NPe+ = 
since there are 4 sites of each type, but each atom is counted z+1 times, being 
either a central atom of asite or in the first shell of z other sites. By differentiating 
(P.F.) with respect to T and H, all the thermodynamic and magnetic properties 
may be deduced. 


Writing @= (= + 2) (y+ =). b= (e+ =) (é if *) and using 


6 and ¢ as parameters, it is possible to find a quartic equation for x)”, which may 
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be solved numerically. In this way a set of simultaneous values of Ho, T, the 
magnetic moment M, and the entropy S is obtained, and by choosing various 
values of and ¢ any number of such sets may be computed. Isomomentals, 
isentropics, etc. may then be drawn by numerical interpolation, in an obvious 
manner. The results of these calculations (for s=8) are plotted in Figures 2 
and 3. The following points should be observed. 


N, 
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Figure 2. 
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Figure 3. 


(a) Outside a line labelled ‘phase boundary’ the equations have only one 
solution, which is obtained by putting =1/x. ‘This corresponds to absence of 
long-range order, so that right-hand and left-hand sites have identical properties ; 
this is the solution quoted by Firgau. Inside the boundary, however, there is 
an additional solution which has a lower free energy, and corresponds to long- 
range anti-parallel order. ‘Ihe two solutions join on the phase boundary, but 
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not smoothly, so that there should be no abrupt change in magnetization or 
entropy, but jump discontinuities in susceptibility and specific heat should appear 
| as the boundary is crossed. The first-order transition postulated by Sauer and 
a Temperley below T,/3 does not appear. 

t (b) In zero external field there is a Curie point at the temperature T', where 
_ J p?/RT = 5 log|z/(z—2)]. At very low temperatures the zero-field susceptibility 
) behaves as i, 
4 Xo = aE exp (— 22d u?/RT) 


Bend the entropy vanishes as Nk exp(—2Jyu2/RT (but these formulae do not 
hold up to the Curie point itself). Above the Curie point the Firgau formula 


is is valid, i.e. iy RT = Ce 


(c) At zero temperature, when the field reaches the value zuJ, the magnetic 
moment jumps from zero to its maximum possible value »N. This corresponds 
_ to the turning over of all the anti-parallel spins into the direction of the field, 
_ which is now stronger than any possible internal field. It is to be noted that the 
_ anti-parallel phase extends for a certain distance above the level of this critival 

field at some temperatures. Although definitely predicted by the equations, this 
_ phenomenon does not seem to be physically significant, but is a result of the 
~ somewhat artificial model-employed. It has been assumed that long-range 

order, and the definite identification of right-hand and left-hand sites, can be 

called into action whenever they are required. But the solutions of the equations 

near the phase boundary correspond to states in which dipoles on both right-hand 

and left-hand sites have average moments in the same sense, although with 
_ different magnitudes. It seems very doubtful whether a state of this kind could 
subsist in a large crystal, with the quantum-mechanical exchange forces always 
tending to give mixtures of many states with the same energy total. In any 
case, aS we approach the transition region the Bethe model becomes more and 
more approximate, and takes no account of the intermediate state, in which 
long-range order has broken down while there is still a large degree of order 
over a range of a few lattice distances from any given site. A line has been drawn 
in the figures showing where the average moment of an atom on a left-hand site 
vanishes and changes sign. This may well be taken as the limit to the region in 
which the low-temperature solution is trustworthy. Outside the phase boundary 
the high-temperature results are also probably correct. 

(d) These curves differ in other respects from those of Sauer and 'Temperley. 
It is not possible, of course, to discuss the effect of the macroscopic shape of the 
specimen, since only nearest-neighbour forces have been included. But it 
should be noted that the lines of constant magnetization all tend towards the 
point H = zyJ, T=0, whereas Sauer and Temperley seem to show the isomomentals 
sloping more sharply away from the temperature axis, and cutting the transition 
curve again at a finite temperature. 

(e) The experimental data are, as yet, somewhat meagre. The only case 
in which a substantial area of the phase diagram has been obtained is Garrett’s 
study of a single crystal of cobalt ammonium sulphate. His results have been 
displayed on a diagram equivalent to, and directly comparable with, Figure 2. 
Such comparison shows very good general agreement for the isomomentals, but 
the entropy, on our model, seems too dependent on the magnetic field. ‘This is 
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to be expected, since we have neglected all entropy save that due to the disordering ~ 
of the magnetic dipole lattice. Our method is also defective in not taking full | 
account of the intermediate-range order, which contributes some negative — 
entropy above the Curie point even in zero field. 

Garrett finds no discontinuities of susceptibility or specific heat. We may 
assume, therefore, that the isomomentals should not meet at an angle, as in our 
figure, but that they should be replaced by smooth curves in the transition region 
(i.e. between the phase boundary and the ‘boundary of the anti-parallel field’). 
The boundary of the antiferromagnetic region, according to Garrett’s definition, 
should then lie on the minima of these smoothed isomomentals, that is, somewhere 
between the two lines enclosing the transition region. 

(f) The number RT,/uH, which Garrett finds to be equal to 1-2 for cobalt 
ammonium sulphate, on our model turns out to be 227! log[z/(z—2)], which is 
slightly less than unity for real lattices. Although this figure may seem to agree 
a little better with experiment that Sauer and Temperley’s value of 2, it is difficult 
to assign any great significance to this fact. Garrett’s salt actually contains two » 
differently oriented dipoles per unit cell, which seriously complicates the analysis. 
Because of the simplicity of the model, it is not very profitable to discuss precise 
values of Curie temperatures, zero-field susceptibility, atomic constants, etc. 
It has already been pointed out by Kurti (1951) that the observed Curie points 
for substances showing these anomalies indicate that the forces are indeed 
magnetic dipole interactions (or are at least of that order of magnitude), and our 
analysis can do no more than support that assertion. We have obtained a rather 
general phase diagram for antiferromagnetics, to which any substance should 
conform qualitatively; more exact treatment, taking into account the nature 
of the forces and the exact crystal structure, must inevitably use more refined 
mathematical techniques, which go beyond the scope of this paper. 


Note added in proof. Since this paper was written, C. G. B. Garrett (J. Chem. 
Phys., 1951, 19, 1154) has applied the internal field method to a model almost 
identical with our own, and obtained a theoretical phase similar to Figure 2. 
The chief difference is in the shape of the phase boundary, which is shown 
‘bulging’ much further beyond the critical field. Since the Bethe method 
is an approximation of higher order, it is likely that the exact solution will not 
show this anomaly at all. 
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ABSTRACT. Silver crystal surfaces which had been prepared by either electrolytic 


etching or electrolytic polishing were studied by means of the electron diffraction reflection 
technique. It was found that fine structure effects, in the form of split spots or streaks 
through spots, gave information concerning the topographies of the various surfaces of the 
crystals. ‘The fine structure effects showed certain asymmetries which do not immediately 
follow from the Laue diffraction theory. The possible causes of these asymmetries are 
discussed, and it is concluded that they were probably due to distortion of the silver crystals, 
but that some might have also been due to refraction effects. 


§1. INTRODUCTION 


LECTRON diffraction reflection patterns obtained from single crystal 
E surfaces are considerably influenced by the topography of the crystal surface. 
Atomically flat surfaces are characterized by the effects due to refraction 
of the electrons upon entering and leaving the surface (Thomson and Cochrane 
1939, Chapter VIII). With other surfaces, diffraction patterns arise when 
electrons penetrate through sub-microscopic projections covering the surfaces 
(Thomson 1931), and the shape and size of such projections can introduce 
important effects. Firstly, the electrons may be refracted when entering or 
leaving the projections by certain facets; Miyake (1938) has reported an example 
of this effect in which the diffraction spots are split up into groups of close 
subsidiary spots. Similar refraction effects have also been found during high 
resolution examination of certain oxide smokes (see, for example, Cowley and 
Rees 1947). Secondly, Laue (1936) has shown that around each reciprocal lattice 
point due to a crystallite with well-developed facets, there is an intensity region 
which can be approximated to a number of spikes or ‘stachel’ perpendicular to 
each of the well-developed facets. When the usual Ewald sphere construction 
is applied to such a reciprocal lattice, the spikes cause some of the diffraction 
spots to be split up into groups of subsidiary spots. If, therefore, a crystal surface 
is covered by many similarly shaped projections, then such groups of subsidiary 
spots should occur on electron diffraction patterns from that surface. ‘The effect 
has been observed several times (Kirchner and Lassen 1935, Briick 1936, 
Cochrane 1936). 

In the present paper, an account is given of the application of electron 
diffraction to the study of various surfaces of single crystals of silver, with the 
object of obtaining information concerning their topography, as a preliminary 
to the study of the growth of chemical compounds on them. Fine structure 
effects are observed, and these can be partially explained in terms of reciprocal 
lattice spikes. Some of the results suggest, however, the existence of spikes 
which are not symmetrical about the centre of the appropriate reciprocal lattice 
point Tbe Laue diffraction theory, which predicts the spikes, also predicts 
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that they should be symmetrical about the centre of the reciprocal lattice point. 
One of the principal aims of this paper is to report and discuss these asymmetrical _ 
effects. 

§2. THE PREPARATION OF THE SILVER CRYSTALS 

Silver crystals were grown by slowly cooling beads of molten silver, of a few 
millimetres diameter, in a continuously evacuated quartz tube. About half of 
the specimens consisted of more than one crystal grain, and were therefore 
rejected. Flat surfaces were ground on those specimens consisting of only one 
crystal grain, by careful grinding on grade 0 emery paper, and finishing on 
grade 0000 paper. The surfaces were etched, and the orientation of the crystal 
axes determined by means of electron diffraction. Flat surfaces could then be 
prepared parallel to any chosen lattice plane, using a specially designed jig for 
clamping the crystals during the grinding. Crystals with flat surfaces parallel 
to the (100), (111), (110) and (311) type lattice planes were prepared, these being 
the four most closely packed planes. 

The final preparation of the surfaces was carried out in one of three ways: 
(i) 35% nitric acid etching, (ii) electrolytic polishing, or (ii) electrolytic etching. 
Nitric acid etching was found to be not very controllable, and to leave the surface 
very uneven, so that only preliminary experiments were carried out with this 
technique. 

The electrolytic polishing was performed according to the method of 
Shuttleworth, King and Chalmers (1947), and very bright even surfaces could 
be obtained after about twenty minutes polishing, with a potential of 0-7v. 
across the cell. The brown film which formed during the polishing could be 
easily brushed off during the rinsing of the specimens in distilled water. 

Electrolytic etching was carried out in the same bath, using a potential of 
0-05 v.; the cell current density was then about 1ma/cm*. The etching was 
usually continued for two hours on each specimen, immediately after the specimen 
had been prepared on grade 0000 emery paper. 

Before examination in the electron diffraction camera, the specimens were 
always washed in distilled water, alcohol and benzene before being dried. 


§3. THE ELECTRON DIFFRACTION APPARATUS 
Two different electron diffraction cameras were used, one having a cold cathode 
electron gun and the other a hot filament electron gun; both were fitted with 
magnetic focusing coils. ‘The electron accelerating potentials were in the range 
30-50kv. and the camera lengths were 40 and 35cm. respectively. The 
reflection technique was used exclusively. 


§4. THE GEOMETRICAL RELATION BETWEEN RECIPROCAL 
LATTICE SPIKES AND THE POSITIONS OF SUBSIDIARY 
DIFFRACTION SPOTS 

During the analysis of some of the diffraction patterns it became necessary 
to calculate the positions on the photographic plate of subsidiary spots due to 
certain reciprocal lattice spikes. A suitable approximate relation has been 
derived. 

Let the crystal surface be cut parallel to a given lattice plane, and let the 
electron beam direction be inclined at an angle 7 to this surface (Figure 1 (a)). 
Let the beam also be set off from being perpendicular to a certain reciprocal 
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plane (yz plane) by the small azimuth angle ¢. If a point P in that reciprocal 
lattice plane gives rise to the diffraction spot P’ (Figure 2), with coordinates 
Y and Z with respect to the undeviated spot as origin, then it can be shown that 
the coordinates of a spot P,’ due to a spike through P are S and U, where 


Ds sii =i — sh tan Bee Sn (1) 
F= (Y?+Z?)—2L(Zsini— Y cosisin ¢) 
(Z tan B—Y sin y)—(sin y cos sin é +siniztan B—cos¢cos y)L" 


B and y define the direction of the spike, as given in Figure 1(b). The xy plane 
is parallel to the crystal surface, and the yz plane is parallel to the above-mentioned 
reciprocal lattice plane containing P; L is the camera length. 


and 


Crystal Surface parallel to 
xy plane 


S —-s p’ 
Ewald Sphere i : 
a 
(ses as y p’ 
y : 
| 
I 
| 
rf \Z 
Definition of 
B and y 
| 
Electron | 
Beam { 
2 of ee x @0 
@ (a) as Central Spot 
Figure 1. Reciprocal lattice and Ewald sphere Figure 2. Definition of S and U. 
construction. 


These formulae may be used to calculate the pattern of subsidiary spots 
due to a given system of spikes, with a given crystal setting relative to the electron 
beam direction. 


§5. THE DIFFRACTION PATTERNS FROM THE ETCHED SILVER 
CRYSTAL SURFACES 


35% nitric acid was first used as etching agent, and this produced matt 
surfaces which gave rise to cross-grating patterns of spots which were sometimes 
split into groups of subsidiary spots, as previously reported (Pashley 1949). 
The etching was not very controllable, however, and the surfaces tended to 
become visually uneven. The electrolytic method of etching was therefore 
finally adopted, and this was found to produce surfaces which gave rise to 
cross-grating patterns of spots often having very marked fine structure. Diffuse 
Kikuchi line patterns were also occasionally obtained, as well as some weak diffuse 
streaks connecting the cross-grating spots along the [111] directions. 


(i) (100) Surface 
Figure 3 (Plate *) shows a reflection pattern from an electrolytically etched 


(100) surface of silver, with the beam set along the [011] azimuth. Some of the 
diffraction spots are resolved into groups of five subsidiary spots. When the 


* For Plates see end of issue. 
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electron beam was set along the [001] azimuth, some of the cross-grating spots . 


were split up into groups of four subsidiary spots. The two groups of spots are 


illustrated in Figure 4. The subsidiary spots can be accounted for if it is assumed 


that through each reciprocal lattice point there are spikes along the following 
crystal directions: [010], [001], [011], [011], [120], [120], [102], and [102]. 
Explanation in terms of a refraction effect such as observed by Miyake (1938) 
does not seem possible in this case. 

The angle of incidence of the electron beam with the (100) crystal plane was 
calculated from formulae (1) and (2), assuming the subsidiary spots C and D 
(Figure 5), for the (600) spot of Figure 3, to be due to [010] and [001] spikes. 
o@ was assumed zero because of the symmetry of the pattern. The values of 
B and y for all of the spikes are listed in Table 1. 


o @ e e @ e 
® 
(a) (b) 

[oil] Azimuth [001] Azimuth 


Figure 4. Groups of subsidiary spots observed from etched (100) surface. 
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Figure 5. Calculated pattern of Figure 3. 


Table 1. Calculated Values of 8 and y 


Spike [010] [001] [011] [017] [120] [120] [102] [102] 
B (deg.) 0 ra 0 26-29 +26:3, 0206-3 063 
y (deg.) +45 45. 0) 4-90. 90445 4a ea 


The expected pattern of spots for the proposed system of spikes, calculated 
according to formulae (1) and (2), is plotted in Figure 5. A comparison with 
the observed pattern of Figure 3 is given in Table 2. The spots C and D are due 
to the [010] and [001] spikes, and the spots A, B, E and F to the [120], [120], 
[102] and [102] spikes. The [011] spike gives rise to the spot G. 

Although the calculated and observed positions of subsidiary spots agree 
well, it is noticed that no subsidiary spots with U positive (spots A and B) are 
observed. It was at first thought that this could be explained by the absence 
of the [120] and [102] spikes; these spikes would be expected to appear, however, 
since the (102), (102), (120), and (120) planes are equivalent and are placed in 
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symmetrical positions with respect to the surface. A photograph (Figure 6 
(Plate)) was therefore taken with 6 =4°, by changing the azimuthal setting of the 
crystal. Here subsidiary spots due to the Ewald sphere passing both in front 
and behind reciprocal lattice points are observed, as illustrated in Figure 7. 
Splitting due to the [120], [120], [102] and [102] spikes is thus recorded on one 


Table 2. Comparison between Observed and Calculated Positions of the 
Subsidiary Spots 


Sub- Sub- 


c diary S (mm.) ~ U (mm.) Bidiaky SS (mm.) U (mm.) 
spot Obs. Caleta Obs o Cale spot Obs Calc Obsas WCale 
(600)A — 1-02 ee Ore) V(622)a: — +1:58 ae 1-10, 
(600)B — ={-02 —==—+0'72, (622)8 ee Sher == 6460-99 
{600)C +1:0 +0-97 0 0 (622)C +14 +1:49 0 0 
(600)D —1-:0 —0-97 0 0 (622)D —1-2 —1-36 0 0 
(600)E +1:0 +0:92 —06 -—0-65 (622)E +14 +1-41 —0-7 —0-99 
e {600)F  —1-:0 —0-92 —06 —0-65 (622)F —1:2 -—1:30 —0-7 —0-91 
(511)A — — +0-68 —-—~+0-48 (133)A — +0:-95 a= 20-67 
(511)B = —(65 = +046 (133)B —  —0-83 =—- 40-56 
{511)C +0:7 +0-66 ) 0 (133)C 0-9 +0-96 0 0 
11D) 0-6. = 0:63 (cet ae) (133)D —0-:8 —0-83 0 0 
fe (511)E  +0-7 +063 —0-4 —0-44 (133)E 0-9 +0:97  —0-5 —0:68 
_ (511)F —06 —060 —0-4 —0-42 (133)F —0:8 —0-84 —0-5 —0-56 
(711)A 1-79 =e 1-25) (333)A — +1:-08 —=) 20-76 
(711)B = — 1-70 — e119 (333)B —  —1:00 = = 0-7t 
(AIC +16 +1-68 0 0 (333)C 1:1 =+1-07 0 0 
(711)D  —1:3 —1-60 0 0 (333)D —0-9 —0-93 0 0 
Mid yes i-6" £157 —1-0: —1-10' <1 (333)B hl -1-41:05.. +06 0-73 
ape. 1:3 - —1-50. —1:0/ —1-05 (333)F ..—0-9. —0-91 —06 0:64 
(422)A ——+0:70  — .+0-49 — (533)A — +1:73 — +1:21 
(422)B == 0-63 —- +044 (533)B ae 1-49 ==) See eod 
(422)C +0-7 +0-69 0 0 (533)\C_ +1:6 +1-66 0 
(422)D —0-6 —0-62 0 0 (533)D° —1:3 =1-44 0 0 
(422)E +0-7 +066 —0-4 -—0-46 (533)E- +1:°6 +1:559 —1:0 —1-11 
(422)F —06 —060 —0-4 —0-42 (533)F —1:3 —1:38 —0-9 —0-77 


° Intersection 
of Sphere 
e e e / with Plane 
Sphere passes in front of Plane 
e 


Figure 7. Intersection of Ewald sphere with (011) reciprocal lattice plane, for pattern of Figure 6. 


photograph, but still no subsidiary spots with U positive are observed. 
Occasionally there were signs of subsidiary spots with U positive, for cases where 
the Ewald sphere passes very close to the reciprocal lattice point; this was more 
noticeable at the [001] azimuthal setting, because the [102] type spikes were 
then more nearly perpendicular to the Ewald sphere. All the effects can be 
explained, therefore, if it is assumed that the [102], [102], [120] and [120] spikes 
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are much longer on the side of the reciprocal lattice point nearer to the origin | 
than they are on the opposite side. It is estimated that the ratio of the lengths | 
of the spikes on the two sides is approximately five to one. i | 

The [010], [001], [011] and [011] spikes showed no such asymmetry. This 
does not necessarily mean that the spikes are symmetrical, since, for the particular 
case of facets perpendicular to the surface, the spike can be due to the superposition 
of two asymmetrical spikes, because two parallel facets facing in opposite directions 
are then possible. 

Apart from the asymmetry effects, the patterns from the etched (100) surfaces 
suggest that the (010), (001), (011), (011), (102), (102), (120) and (120) facets 
are well developed during the etching process. A later paper will show how the 
presence of some of these facets is also indicated by the way in which silver halides 
orient on the same surfaces. 

(ii) (100) Surface 

When the (110) silver surface was electrolytically etched, the diffraction spots 
were not split into groups of subsidiary spots but had various shapes. At the 
[110] azimuthal setting the spots were often diamond shaped, being 
predominantly streaked down towards the shadow edge. In Figure 8 (Plate) 
the spots are star shaped, with three sets of streaks through them. At other 
azimuthal settings other streaked spots were observed. ‘The patterns obtained 
varied in detail from experiment to experiment, but the spots always had streaks 
through them, and these streaks were always perpendicular to certain crystal 
lattice planes. At various times streaks perpendicular to the following lattice 
planes were observed : (110), (110), (100), (010), (001), (021), (021), (201), (201), 
(210), (210), (120), (120), (113) and (113). 

The streaks were often asymmetrical about the centres of the Laue spots, 
being stronger on the side nearer to the central spot. The more nearly the streak 
was perpendicular to the shadow edge, the more was it asymmetrical. Streaks 
parallel to the shadow edge were symmetrical. 

(ui) (111) Surface 

When the (111) surface was electrolytically etched, two distinct types of 
surfaces were obtained, according to conditions. If the crystal was cut with its 
surface accurately parallel (to better than one degree) to the (111) lattice plane, 
the electrolytic etching gave a bright flat surface very similar to a cleavage face. 
If the crystal was not so accurately cut, or if the surface was given a previous 
etch in nitric acid, then a matt surface was obtained. 

The bright flat (111) surface gave rise to diffraction patterns such as shown in 
Figure 9 (Plate). This is similar to a pattern from a cleavage face, and indicates. 
a high degree of surface flatness. Rotation photographs taken from these surfaces 
showed that the diffraction spots were displaced down by refraction during the 
entrance and exit of the electron beam from the surface. An inner potential for 
silver of approximately 20 Vv. was indicated. This is consistent with the value 
of 23-2 v. recently obtained by Tull (1951). It is thus concluded that this type of 
(111) surface contains areas which are atomically smooth, in contrast to other 
etched surfaces which are usually rough. It would appear that there is a 
preference for silver to etch with extensive (111) planes as its boundary, since 
nitric acid etching of the silver crystals caused the appearance of bright flat. 
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patches of about one miilimetre in dimension parallel to the cee lattice planes 
(Pashley 1949). 

The matt (111) surface gave cross-grating patterns of spots which were 
heavily streaked asymmetrically, as shown in Figure 10 (Plate) for the [211] 
azimuthal setting; the streaks are in [120] and [102] directions. Occasionally 
_ there were signs of the appearance of subsidiary spots, but there is insufficient 
evidence on this point to allow any definite conclusions to be drawn. At other 


| azimuthal settings, streaks perpendicular to the following lattice planes were 


observed : (012), (102), (210), (201), (110), (101), (011), (113), (131) and (311); 
all showed asymmetry. 
(iv) (311) Surface 
When the (311) surface was electrolytically etched, patterns showing fine 
_ structure effects were also obtained. Similar subsidiary spots to those observed 
from the (100) surface were sometimes found, and asymmetrical streaks 
perpendicular to (110), (311) and (210) type facets were also observed. 


§6. THE DIFFRACTION PATTERNS FROM THE 
ELECTFROLYTICALLY POLISHED SILVER CRYSTALS 

Electrolytically polished silver crystal surfaces gave rise to diffraction 
patterns consisting of the following components: (i) cross-grating patterns of 
spots which were usually diffuse and strongly streaked out towards the shadow 
edge, these spots were due to the main silver crystal; (ii) weaker spots due to 
silver which was twinned on the (111) planes of the main silver crystal; (iii) weak 
sharp ring patterns; (iv) diffuse Kikuchi line patterns. Figure 11 (Plate) shows 
a typical example, with prominent twin spots present. 

The cross-grating patterns were usually quite extensive, and often two Laue 
zones of spots were considerably overlapped; they must therefore have been due 
to electron penetration through very small projections on the surface. The 
breadth of the diffraction spots indicated that the projections were probably 
no larger than a few atomic dimensions across. ‘The surfaces were certainly 
not atomically flat, and even when an accurately cut (111) surface was 
electrolytically polished, the surface was also covered by small projections, 
despite the fact that electrolytic etching of such surfaces produced atomically 
flat areas (see §5). It is therefore concluded that the mechanism of electrolytic 
polishing does not lead to the formation of an atomically flat surface for the 
particular technique used. The spots were always elongated downwards; 
this effect has been observed previously (Kranert, Leise and Raether 1944, 
Raether 1946, Halfawy 1948) and Kranert et al. have explained it as due to a smooth 
crystal surface which is wavy on an atomic scale. 

The electrolytically polished silver crystal surfaces thus gave rise to 
diffraction patterns which also indicated some asymmetrical etfects. 

The presence of twinned silver suggests that some silver is re-deposited 
during the electrolytic polishing; this is supported by the fact that electrolytic 
deposition of silver on to the crystals is also found to lead to twinned growth. 
The presence of randomly oriented silver on the electrolytically polished surfaces 
might also be due to re-deposition. The twinned and randomly oriented silver 
could be removed by less than five minutes’ electrolytic etching, and was therefore 
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confined to a very thin layer on the surface. ‘The re-deposition might be 


connected with the brown film which formed on the crystal surface during | 


electrolytic polishing and which disappeared when the surface was brushed with ' 
water after the crystal had been taken out of the electrolytic cell. 

The diffuse Kikuchi line patterns varied in brightness from specimen to 
specimen. 


§7. DISCUSSION 


Whereas electrolytic polishing caused the silver crystal surfaces to become 
smooth and flat, the electrolytic etching, with one remarkable exception, produced 
rough surfaces which were matt in appearance and which were bounded by 
certain well-developed facets. It is interesting to note that the facets were always 
of the (100), (110), (210) and (311) types, and that no (111) oblique facets were 
ever observed, although the octahedral planes are the most closely packed, and 
might have been expected to be the most stable. Cochrane (1936) observed that 
prominent (111) type facets were produced on some electrolytically deposited 
cobalt surfaces, and Kirchner and Lassen (1935) observed well-developed 
(111) facets on condensed layers of silver. In the present experiments (111) 
crystal boundaries were only produced when they were parallel to the surface; 
they then formed large flat patches on that surface. It is also interesting to note 
that if the (hk/) facet occurs on a given surface, all of the crystallographically 
equivalent facets do not necessarily occur, but only those which make the same 
angle with the surface as does the (hR/) facet. 

With the electrolytically etched surfaces there is the difficulty of explaining 
the observed asymmetrical effects. ‘These effects were of two distinct types: 
the subsidiary spots, shown in Figures 3 and 6, which were due to asymmetrical 
reciprocal lattice spikes, and the asymmetrical streaks as shown on Figures 8 
and 10. 

Asymmetrical intensity regions around a reciprocal lattice point of silver are 
inconsistent with the Laue diffraction theory for an undistorted lattice. The 
intensity of a point (h,, fz, 3) in the reciprocal lattice is proportional to 
Dinan, = LL*SS* where L,,,,,,= % exp 2nt(lh, +mh, +nhg) 


l,m,n 


and Span a =f,(9) exp 27t(hyg, + hogs + h;23). 
g 


(/, m, n) are the coordinates of the origins of the unit cells in the crystal lattice, and 
(21, £o) 3) the coordinates of the individual centres within the unit cell, referred 
to the cell origin. /f,(@) is the atomic scattering factor of those centres. 
LL* for h,=h,' +4, h,=h, +8, hz=hs'+y is equal to LL* for h,=h, oe 
ho =h,' —B, hg=h, —y provided hy’, h,', h,’ are integers. LL* is therefore 
symmetrical about each reciprocal lattice point. Although SS* can, in the 
general case, be asymmetrical about certain reciprocal lattice points, this is 
not always so, and for the case of the face-centred cubic structures (e.g. silver) 
SiS* is symmetrical about all allowed reciprocal lattice points. In the present 
experiments, there should thus be complete symmetry of the Laue spikes about 
each allowed reciprocal lattice point. The following explanations have therefore 
to be considered: (i) refraction effects; (i1) a change in crystal structure near the 
surface, such that the cell structure aio S can cause ge B NSIS (iii) crystal 
imperfections and distortions. 
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Although asymmetrical reciprocal lattice spikes do not appear to have been 
previously observed, one-sided streaks through diffraction spots have been, and 
Yamaguti (1935) considers the streaks to be due to refraction occurring as electrons 
enter and leave surface projections via certain facets. He does not explain, 
however, why refraction gives rise to streaks instead of subsidiary spots, as 
would be expected. Figure 12 shows the subsidiary spots which would be 
expected to appear on Figure 10 due to refraction at the (120) and (102) facets. 
The refraction displacements are very sensitive to the angle of incidence of the 
electrons with the facet which causes them, so that small variations in the angle 
of incidence would tend to cause the subsidiary spots to streak along [120] and 
[102], and perhaps coalesce to give a continuous streak. The variation in angle 
of incidence could perhaps be caused by facets with wavy surfaces, or might also 
be caused by a mosaic crystal which gives rise to facets which are rotated through 
small angles with respect to one another. The diffraction patterns indicate that 
the angular range due to this latter cause is less than one degree. It must also be 
remembered that a slightly convergent electron beam is used. It might therefore 
be possible to overcome the objection that separate subsidiary spots should be 
observed if refraction effects are appreciable. 
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Figure 12. Subsidiary spots due to facet refraction, for pattern of Figure 10. 


In the present experiments there are three further objections to the 
explanation that the observed asymmetrical streaks are due to refraction: 
(i) Although predominantly one-sided, the streaks pass through the diffraction 
spots; the shorter upward part of the streak cannot be explained by refraction. 
(ii) Except in special cases, most of the electron beam should enter and leave 
the crystal via facets which make large angles with the beam direction. Only 
a small part of the beam would thus be expected to suffer refraction, and streaks 
as strong as those observed would not be expected. Miyake (1938) only obtained 
strong subsidiary spots due to facet refraction when he set his beam several 
degrees off azimuth, so that much of the electron beam passed through the 
appropriate facets. (ili) A combined effect due to refraction at two different 
surfaces, one during entrance and one during exit, would be expected to occur. 
No such effects are observed. 

Refraction cannot, therefore, entirely explain the observed asymmetrical 
streaks, but it might contribute strongly; thus the asymmetrical streaks which 
pass through the Laue spots (see Figures 8 and 10) might be due to the 
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superposition of a streak arising from a Laue spike and a streak due to refraction. 
That the streaks on Figure 10 can be entirely due to asymmetrical spikes, which 
are tangential to the Ewald sphere, does not seem very likely, since such long 
spikes would be expected to give rise to prominent subsidiary spots at certain 
azimuthal settings, and no such spots were observed. 

The cell structure factor S, as previously mentioned, can sometimes introduce 
asymmetry into the reciprocal lattice, but not with a face-centred cubic structure. 
If, however, the lattice structure near the surface were changed (e.g. to a 
hexagonal lattice) then asymmetry might result. The change would have to be 
such as to introduce no extra reciprocal lattice points, and for this reason an 
explanation on these lines is considered unlikely. 

The final effect to be considered is that of crystal distortion. If a crystal 
lattice is distorted, then the symmetry of LL* is no longer maintained. A simple 
example to consider is that of a crystal (Bravais lattice) with one set of plane 
spacings changed so that the distance between the mth and (n+1)th plane from 
the surface is given by af(n), where a is the normal spacing of the set of parallel 
planes, and f(z) is some arbitrary function. Consider a rectangular block of 
crystal with the crystal axes parallel to the edges of the block, and let the 
(001) plane spacings be modified. Then the Laue factor may be written as 

L= % exp2mi(lh, + mh.) X exp 21 f (n)heg ; 

l,m Nn 

f(n) is arranged so that the spacings of the (001) planes near the surface are larger 
than normal, and other parallel planes gradually approach the normal value as 
n increases. It can be shown (unpublished work by Mr. Skwirzynski) that the 
spike due to the (001) facets will be asymmetrical, being much longer on the 
side nearer the origin of the reciprocal lattice. A qualitative picture of the result 
is given by considering that the gradual increase of the crystal spacings causes a 
gradual decrease of the corresponding reciprocal lattice distances, so that the 
diffraction spots are streaked on the side nearer to the central spot. The change 
in plane spacing may be regarded as a non-linear stretch of the crystal lattice 
along a direction perpendicular to the particular crystal planes. 

An increase of the spacings of the planes parallel to and near to a surface | 
facet would therefore cause the spike due to that facet to become asymmetrical. 
An important result which follows from the theory is that the reciprocal lattice 
should contain a plane across which the asymmetry should reverse in direction. 
This reciprocal lattice plane is perpendicular to the direction along which the 
crystal lattice is stretched, and passes through the origin. The spikes on either 
side of this reversal plane will be longer on the side pointing towards the plane. 
The reciprocal lattice will still be symmetrical about its origin. 

In order to apply these calculations to the explanation of the observed case 
of the asymmetrical (210) type spikes, it is necessary to choose a unit cell with no 
structure factor (Bravais lattice) so that two of the principal crystal axes are in 
the (210) plane and the third perpendicular to it. All of these conditions cannot 
be fulfilled simultaneously, and the third crystal axis has to be chosen (for 
example) as the [110] axis (referred to the cubic frame). The other two axes 
are the [121] and [001], both lying in the (210) plane, and the unit cell only 
contains one atom (Figure 13). This means that the crystal is stretched along 
the [110] axis, which is not perpendicular to the (210) planes, so that the 
(210) planes are also displaced parallel to one another during the stretching. 
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_ This type of distortion should lead to the (110) reciprocal lattice plane being a 
plane of reversal for the asymmetrical (210) spikes. Since this plane of reversal 
is not observed on the diffraction photographs, non-linear stretching along the 
[110] crystal axis does not supply a sufficient explanation of the asymmetry. 
_ Other stretching directions which can be considered by this simple treatment 
also do not fit as regards reversal of the asymmetry. The treatment of a 
non-linear stretch in any general direction is more difficult to deal with, since the 
stretching would have to be resolved into its three components along the principal 
crystal axes, and a triple summation evaluated. This would not lead to a simple 
rule for the plane of reversal. 


ft) 


Figure 13. 


The important result of these considerations is that lattice distortion can 
cause the Laue spikes to become asymmetrical about the centre of the reciprocal 
lattice point, and it seems reasonable to suppose that a certain distortion, not 
necessarily a simple stretch, would probably lead to the observed asymmetry. 
It is not easy to decide which distortion has to be chosen, because each choice, 
which is to a large extent a matter of guesswork, would have to be checked for 
such details as the occurrence or otherwise of planes of symmetry reversal in the 
reciprocal lattice, and for the non-appearance of extra reciprocal lattice points. 

Summarizing the conclusions which follow from the above discussion, we 
‘can offer the following explanation of the asymmetrical effects: (i) the 
asymmetrical spikes are due to facets on the surface, the asymmetry arising from 
lattice distortion; (ii) the asymmetrical streaks observed on many of the 
diffraction photographs are due to the superposition of effects due to the refraction 
of electrons upon entering and leaving the crystal surface via certain facets, and 
effects due to Laue spikes arising from the same facets. It is possible that the 
asymmetry effects observed with the electrolytically polished surfaces are also 
due to crystal distortions near the surface. 

No explanation of the asymmetry in terms of dynamical interaction or other 
effects such as absorption has so far proved to be satisfactory. 

If a change in crystal spacings near the surface is sufficient to account for 
the observed asymmetry, then this could possibly arise if, during the electrolytic 
treatment of the crystals, some foreign atoms should become adsorbed at the 
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surface, and pack between the first few atomic layers. It is well known, for 
example, that silver readily dissolves oxygen. Any such adsorbed atoms, — 
particularly if they are of low atomic weight, would otherwise be difficult to ') 
detect by electron diffraction. 

It is interesting to note that asymmetrical effects in the form of streaks 
through spots have been observed during x-ray studies in aluminium-—copper 
mixtures (Preston 1938a, b, Guinier 1938). Preston (1938a) has suggested 
that this asymmetry arises from a varying lattice parameter accompanied by a 
varying structure amplitude. 

Laue (1937) has stated that the appearance or otherwise of subsidiary spots 
(due to reciprocal lattice spikes) around a normally forbidden spot should give 
some information concerning the fine structure of the facets. With the present 
experiments, the spikes were so short that it is considered that any spike 
surrounding a forbidden spot would be too weak to observe. No conclusions 
can therefore be drawn about the finer details of the facets. 
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ABSTRACT. A very general calculation has been made by Fréhlich giving the dielectric 
constant of a material in terms of the spontaneous fluctuations of the electric moment of the 
material which would occur in a field-free region. This calculation is now extended one step 
further so as to include the first term in the field dependence of the dielectric constant. By 
way of illustration the general theory is applied to Kirkwood’s model for a dipolar liquid. 


(In order to preserve continuity of argument it has been found necessary to reproduce much 
of Fréhlich’s work.) 


§1. INTRODUCTION 

TTEMPTS to set up a theory of the field dependence of the dielectric 

constant of a material have nearly always related to water, since the most 

obvious practical application of this saturation effect is that of finding the 

field strength in the neighbourhood of ions in an electrolyte. Debye (1928) 

attempted to estimate this effect by applying Mosotti’s hypothesis, but the 

experiments of Malsch (1928) showed a very much smaller effect for water than 

that anticipated by Debye’s estimate. More recently Booth (1951) has given 

a treatment of the saturation effect using Kirkwood’s analysis, and has applied 

the results to water. This application, however, contains a mistake, as has been 
pointed out elsewhere. 

In the following paper Frohlich’s (1949) calculations for the dielectric constant 
are extended so as to include the first term of the field dependence. A general 
method is developed for dealing with this correction term, which is then applied 
to Kirkwood’s (1939) model of a dipolar liquid. 


§2. CALCULATIONS 
2.1. Preliminaries 


Following Frohlich we consider an infinite, homogeneous, isotropic dielectric 
which is not permanently polarized. This is split up into equal regions, each of 
which is just large enough to have the same dielectric properties as a macroscopic 
specimen. A spherical region containing a large number of these smaller regions 
is then considered, and the dielectric constant calculated in terms of the 
spontaneous fluctuations of the large sphere embedded in its own medium. 
The particles of the sub-regions inside the sphere will be treated according to 
classical statistical mechanics, but the outside will be considered as a continuous 
dielectric of macroscopic dielectric constant 


Gate th a ee eee nee (1) 
where «¢, is the field independent static dielectric constant. 


* Based on E.R.A. Report L/'T 258. 
+ Now at National Standards Laboratory, Sydney, Australia. 
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Suppose that the large spherical region consists of a number of charges — 
e,, each of which can be described in terms of its displacement r, from its position _ 
in the ground state of the whole system. If we put 

X= (Ly; fae ese op Ee eee) eal a es 
then the average value of the projection of the moment of the sphere in the 
direction of the macroscopic field will be 
Mu _ JM(X) cos 6 exp {— U(X, E)/RT} dX (2) 
BE expt =UGG ERT {2X |) See 

where U(X, E) is the potential energy of the system in the configuration X in 
the presence of a field EF, M(X) is the moment of the sphere in terms of the charge 
displacements, and 0 is the angle between M(X) and E. 


2.2. The Cavity Field 
The potential energy is given by Frohlich as 


UX, E)=UX)— 5 ats FM(X)Ecos8 sees (3) 


where U(X) is the energy in the absence ues a macroscopic field, and the second 
term represents the interaction of the cavity field with charges of the spherical 
region. This second term will not be correct for our case as the dielectric constant 
of the medium is now a variable function depending on the field strength. If 
it is assumed that equation (1) can be written as 
e= e,1 OL) oer elk pen ree (4) 

where 6 is a small constant, then the first order approximation to the cavity 
field can be obtained by expansion of the potential in spherical harmonics. The 
cavity field will, of course, consist of a uniform and a non-uniform part, and it 
can be shown that the latter is of a smaller order than the correction in the former. 
In fact to a first order the cavity field is 

G= BeqiZeorel JES = aa ek ee Eee (5) 
where E,, is the macroscopic field at a great distance from the cavity and 


Se 2 
c{148 (5-5 Ea) \. 
We may then give the potential energy as 


U(X, E)=U(X)— 5 = 


{ME COS:0. 2 sien Laas (6): 


2.3. The Fundamental Equation 
In view of (6) we make the expansion 
exp {— U(X, E)/kT} 
=exp {— U(X)/kT} E aC: oe a cos 6 
C? M*(X)E? CG Sage 
D Repe cos" A+ Ra Tea cost 0+... |, 
where C=34)(2e4+1) . a eee (8) 


varies very slowly with E compared with the term of which it is the coefficient. 
Substituting (7) into (2) and noting that terms of the form 


JM"(X) cos” 6 exp {— U(X)/RT} dX 
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f q will vanish if m is odd (they are proportional to an odd power of the moment in 


the E direction in the absence of a macroscopic field) we find 
(C/kT){M?(X)E cos? 6 exp {— U(X)/kT} dX 
+ (C3/6k?T*) [ M4(X)E° cost 6 exp {— U(X)/RT} dX 
~ Jexp {— U(X) kT} dX + (C?/2k?T?) | M?(X)E? cos? 6 exp {— U(X)/kT} dX 


retaining only terms of order up to E*. Expanding the denominator by the 
binomial theorem 


a eal M?(X) cos? 6 exp {— U(X)/kT} dX 
zn — i M*(X) cos! 6 exp {— U(X)/kT} dX 
: me — | [ 122(X) cos? 0 exp {— U(X)/R Tax | scat (9) 
where Peete = UX RT). oe Phe (10) 


Since the field E may have any direction relative to M(X), we must average over 
all values of 6. This gives 


ey 3 ae l C3 ES ange M?2)2 
M,= CME 50 C aa 8M ees) Gt Ra Meer et (11) 
where M? and M# are the mean square and fourth power respectively of the 
spontaneous polarization of a sphere of the dielectric embedded in a large 
specimen of the same material. 

Equation (11) is the basic result from which we calculate the dielectric 
constant. Before proceeding to special cases, it is necessary to develop 
equation (11) further and show that, in the general case, it does lead to a 
polarization which is directly proportional to the volume. 


2.4. Discussion of Ve 


Suppose now that the large sphere is composed of N units (molecules or other 
groups of atoms) in such a way that each unit has on the average the same 
polarization in an external field. Let us suppose also that each unit contains 
the same number, k, of elementary charges which in each case can be arranged 
: a similar way relative to each other and relative to their surroundings. Let 

=(f)1, To, . .. - %j,), denote the whole set of displacements in the jth unit, and 
ACs) the electric dipole moment of this unit. Then 


N N 
M(X)= = m(x,;) and M*(X)= % m(x)M(X). —...... (12) 
j=1 j=1 
nes N . 
Hence by definition M?= & J | m(x,)M(X) exp {— U(X)/RT} dX. ...... (13) 
j=1 


In the jth term of this summation the integration is carried out firstly over the 
the whole spherical region except the jth unit, and subsequently over this unit. 
This gives Frohlich’s result 

| m(x;) m*(x,) p(x,) dx;, ee ray 


4= 1 
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M(X) exp {— U(X)/kT} dX; 
where we have put m*(x;) = ee Ena ee Fagen en (15) 
pfu) = LexP{= UXYRT} aX, 
? fexp {— U(X)/kT} dX 
It will be seen that m*(x,) represents the average moment of the whole sphere» 
for a fixed set x, of displacements of the jth unit leading to a moment m(x;), and 
p(x;) is the probability of finding the jth unit with this particular set of © 
displacements. 

Practically all terms in the summation (14) will be equal. The quantity 
m*(x,) represents the average moment of the large spherical region when the 
jth unit is fixed with a moment m(x,). Now by a simple theorem in electrostatics 
the same moment m*(x,) is contained in any other sphere surrounding the 
jth unit, provided only that this sphere is sufficiently large to be treated 
macroscopically. Hence m*(x,) differs from m(x,) because of short range 
interaction, and is independent of the position of the jth unit provided that it is 
not too near the surface of the sphere. The number of units for which this will © 
not hold can be made very small compared with the total, so we may write 

MN ma, nee” ae era (17) 

where m.m*= fm(x,;)m*(x;) p(x) dx, se ee (18) 
where j may refer to any of the units except those near the surface of the sphere. 
In view of the above interpretation of m* equation (18) may be written somewhat 
differently. Using (15) and (16), together with the fact that short range - 
interaction determines m*, we may write 
Jma(x;)M(Xj) exp {= U(AG)/RT} dX, | 

jexp d= UC ede ne 
where j again may refer to any of the units except those near the surface and 
X= (py Xpgay s 2+ + Xjys 00+ Xm—1> Hm) denotes the set of displacements of any 
set of units so chosen as to include the jth and all those units in short range 
interaction with the jth, but subject to no other restrictions. 


m.m*= 


2.5. Discussion of {3 M4—5( M?)?} 
NEN, 
We may write M*= x & {m(x,).M(X)}{m(«,).M(X)}, 
k=11=1 
in virtue of equation (12) so that 


3Mi-3 3 =’ (m(w,) -M(X)) (me Ma) 


+35 1 __, fae) MOG) MOQ, os. (20) 


the averages being over all configurational coordinates. The first double 
summation in (20) accounts for the case in which / refers to a unit which has at 
least one unit in common short range interaction with the kth (there are assumed 
to be 7 such units excluding the hth itself), while the second double summation | 
covers all other cases. | 
To perform the average of the second term of (20), we divide the large sphere 
into two regions 1 and 2, such that region 1 includes all units from 1 to s including | 
the Ath and all units in short range interaction with it, and region 2 includes all 
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| units from s+ 1 to N including the /th and all units in short range interaction with 


Mepit-” Then let X,=(x,,....%3,....«,) and X,=(%,1,....%).... #y) denote 
___ the set of configurational coordinates as shown, and we may write for the potential 


energy in the absence of a field 
OX) =A Xa Xs) = U(X, )AUK) tates (21) 


___ if a small term (surface interaction energy between regions 1 and 2) be ignored. 


I Also MEMO emer) 8 ar (hh i rien (22) 


Using (22) in the second term of (20) yields four terms of which two vanish on 
averaging to give 
N ye ee eS nt 
35S fm(e,). MX} m=). MO} 


7. 


3% “2 [im@). MOSH ame) MED) 


k=1 l=1,k+n 
+ {m(x,) .M(X2)} {m(x,) . M(X,)} | 
SaaIN Nomar IC ete (23) 
since on using (21) in the averaging we find 
{m(x;,) . M(X})} (ma(«,). M(X,)}=(m . m*)? 
and {m(«,) - M(X2)} {m(x,) . M(X})} 
={m(x,) . m(x,, x,4,)} {ma(x,) - M(x, Kp4r)} 
= 3{m(x,) . m(x,, X,4,)} {ma(#,) . M(x), %,4,)} 


=+(m.m*)? 


where m(x,, «;,,) denotes the total moment due to the jth unit and all those in 
short range interaction with it. 
To perform the average of the first term of (20) we write 
M(X) = m(«,, %,5,) +M(X=%)4, X441) sed A, 2 

where (%,,%,,) denotes the set of coordinates (%,, %,41,-+-+%n4,) and 
(X—x,, %,,,) denotes the set of all other coordinates. Using (24) in the first 
term of (20), and performing the summation over /, we obtain four terms of which 
two vanish on averaging to give 


3 2 fale) MEX) - MX} 


= Aa {M(x p,) «WAH py X45) }MO(X py Wen) 


NP Sol og Ee ee ee ce 
+3 ae {m(x,) « M(X — X45 Xp 4n)} {M(H us Xn) M(X — 2%}, Xpan)} 


The second term of (25) can be simplified by performing the Pilgwene menos 


=Nm.m*(NV—=y—1)m.m*  —— avaeas (26) 
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Then using (23), (25) and (26) in (20) we have finally 
3 M4—5( M2)? 
N Picola 2 sae Nes 2 Ltt 
a8) oe {rn (20 5.) «TIC py 5,45) }O(% py B+ n) 


—5 Mat Damm) Se (27) 
and the expression (27) varies linearly with the volume of the macroscopic sphere 
considered. Further analysis of the first term of (27) requires the help of an 
explicit model. 

2.6. The Dielectric Constant 
Returning to equation (8), we write 
eps Gane OE 
ODE eer ee te 
introducing an error which can be shown to be at least an order of magnitude 
smaller than the correction being calculated. Since 


4a Mz 
c= 1 = E ~V. se eeee (29) 
by definition and Frohlich gives ee 
4a 3e, M? 
Sia. ie 3V 2e,+ 1 RT’ so0oKed (30) 
we have on inserting (28), (29) and (30) in (11) 
val 4a 1 3€5 3 ‘Ee 774 7 72\2 
€=e,+ 750 (a4) pReM — 5(M?)*}, cece (31) 


where by virtue of (27) the correction term in (31) is independent of the volume 
of the macroscopic sphere chosen. 


§3. SPECIAL CASES ' 
Equation (31) will be applied to two cases in which (27) can be calculated with 
the help of a detailed model. 
3.1. Assembly of Non-Interacting Dipoles 


In this case it is easily seen that (27) becomes 


3M4—5(M2)?=—-2Nu ae (27 a) 
where there are N dipoles each of moment p, and (31) gives 
Sai Se, Nee 
=e Seeey santa ee 
€=€3— Ze (5 i) Bp Not afte (31 a) 


where there are N, dipoles per unit volume. (A simple direct calculation verifies 
this result.) 
3.2. A Dipolar Liquid (Kirkwood’s Formula) 


An attempt to extend Kirkwood’s formula for polar liquids can be made as 


follows. ‘The liquid will be considered as a continuous medium of dielectric 


constant m* containing dipoles of moment wu. (different from the vacuum moment 
p,, but connected with it by 
ty = Spl te? 2) ee ese (32) 


viz. 
Hig - : . 3e,  47Nop 


\v 


kaa 


iw far WER ance PRR 
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_ for the case of spherical molecules). This amounts to assuming that the 


molecular polarizability is due to electronic displacement. If, in addition, we 
assume that the main part of the dielectric saturation is due to dipole orientation 
and not to electronic displacement then at will become 
471 ( 3e 3s 

$ A 2)\2 
750 tt) we aM = 5( 472} Lame te (31) 
Here «, is calculated from the specialized Kirkwood formula given by Frohlich, 


S S66 


Clee Cosy) aenn 6 aes (33) 


Es 


Wr degeoe SRT. 
where z is the average number of nearest neighbours and cos y is the average 
of the cosine of the angle between neighbouring dipoles. It is assumed in the 
derivation of (33) that only nearest neighbours make any contribution to m*. 
When calculating the average of the first term of (27) on the Kirkwood model 
we may, to a good approximation, take y to refer to nearest neighbours only. Then 


{m(x,) - M(X;,, % pn) }MV(X py XH 4n) 
=p4(1+2 cos y)(g+1+2z2 cos y+42(z—1) cosy’) ...... (34) 
where cos y’ is the average of the cosine of the angle between dipolar directions 
of two nearest neighbours of a given molecule. For the second term of (27) we 


have at once oe — 
(this ea (les Cosiy)tewl vi 8 mci tic (35) 


so using (34) and (35) in (27) and applying the result to (31b) we find for the 
extension of the specialized form of Kirkwood’s formula 


4rN, 3€, 


+ (727 +2) cos y — 62? cos? y +52*(z+1)(cos y)? 


ercat b)cosy (lis cosy)|/ °° a ls 


§4. APPLICATION TO EXPERIMENTAL WORK 


Since there is some experimental work available on the subject it seems 
desirable to apply the results of (31c) to a calculation of the dielectric constant of 
water. Using the model of Bernal and Fowler (1933) in which each water molecule 
has four nearest neighbours at the vertices of a regular tetrahedron, and the 
bond between neighbouring molecules is directed along the O-H bond of one 
molecule towards the oxygen ion of the other, we find by elementary 


3, cos? y=4, cos y’ =x, which gives 


trigonometry cos y= 
2 2 
_14q f ye heap (=H) PENG eee aes (31d) 


3 3 kT 3 Rei? 


where (32) has been used, and the approximation 3¢,/(2e,+m°?)=3/2 made on 


account of the high dielectric constant of water. 
Substituting numerical values: Avogadro’s number =6 x 108, molecular 


weight = 18, k=1-4 x 10-1, »,=1-9 x 107** and n? =1-77 we find 


2 
ie (31.6) 
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or at room temperature (T'=300° x). 
€=63-43x10-VE2 nn ae we (317) 


where £ is now in volts per cm. 
If we define the incremental dielectric constant as 0D/dE, then for water 
at room temperature 


dDiob=63-129%10 2 Eee (36) 


On the basis of Debye (1928) theory, Malsch (1927) calculated that he should 
obtain about 14°% decrease in the incremental dielectric constant of water at 
room temperature in a field of 100,000 volts per cm., but in his initial experiments 
he failed to detect any effect. In later experiments (Malsch 1928) he found that 
the dielectric constant did, in fact, saturate according to a law of the form of 
(36), and he measured a 0-7% decrease in the incremental dielectric constant of 
water at room temperature, in a field of 250,000 volts per cm. The present 
theory predicts a decrease of 1:2°% in the incremental dielectric constant at this 
field strength. ‘This agreement must be regarded as satisfactory in view of the 
approximate nature of the theory and the scarcity of relevant experimental work. 
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LETTERS TO THE EDITOR 


Approximate Solution of a Transient Heat Flow Problem 


The theoretical solution of a geophysical heat conduction problem has been obtained by 
approximate operational methods; it is thought that the results and the procedure used may 
be helpful in other types of ‘ diffusion’ problems with analogous boundary conditions. 

The problem is that of measuring the thermal constants of an infinite medium of low 
conductivity by the following method. A thermal ‘probe’ is inserted into the medium 
containing a source of heat and some form of thermometer, and a record obtained of 
temperature within the ‘ probe’ versus elapsed time. The thermal constants are to be 
evaluated from this record. 

A similar method has been used before (Hooper and Lepper 1950), but the simplifying 
assumptions made in the theory are too extreme for the present problem. The practical 
restrictions on the latter are: (i) the probe cross section must be circular and its outside 
diameter of the order of 4 cm. or greater, (ii) thermal contact resistance at the boundary 
between probe and external medium is likely to be high, and must be taken into account. 
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By using both ‘large time’ and ‘small time’ approximate solutions for the temperature 
in the probe it becomes possible to measure the contact resistance under the conditions of 
the experiment and to eliminate its effect. 

The problem, basically, is that of solving the heat equation for an infinite region 
bounded internally by a (hollow) circular cylinder with the standard contact 


_ resistance boundary condition at the probe surface (Carslaw and Jaeger 1947 a, pp. 17-18), 


— K,06,/@p=H(0,—6,) and zero initial temperature. 
The differential equation and boundary conditions have been subjected to the Laplace 


. _ transformation and large and small time approximate solutions (Carslaw and Jaeger 1947 b, 


Goldstein 1932) obtained for a number of probe arrangements. 

Examination of the possibilities shows that a hollow metal probe with constant heat flux 
across the inner surface is both experimentally practicable and convenient for calculation of 
results. 

With the simplifying assumptions that the probe wall is perfectly conducting and that 
heat flow is entirely radial, the following results are obtained : 


Large times (in the limit a+b, H +, this solution reduces to one given by Carslaw and 
Jaeger (1947 a, p. 283)): (T=h,7t/b? > 1) 


aQ 2 __ ahs 2K 1 
6,(t) = | i n4T—y+ + yp4ln ar y+1 im, (Ina y4 a) + +0 (x3) |. 


Small times: (T=h,?t/b? < 1) 


ye Pipest6r li y= 
~~ = a 2 1 3 
6,(t)= t =H vm) * | Océ?) |, 


where 6,(¢) is the probe temperature, t is the time, Ky, h,? the thermal conductivity and 
diffusivity of the external medium, H the ‘ outer conductivity’ at probe surface (a measure 
of the contact conductance), Q the heat flux (constant) across inner surface of probe, a, b the 
internal and external probe radii, ~a=(b?—a?)8,c,/b, where 5,, c, are the probe density and 
specific heat, y is Euler’s constant=0:57722... 

Correction for finite probe conductivity is not difficult and may be necessary unless the 
probe walls can be made thin. 

An upper limit can be placed on the error introduced by finite probe length. Restric- 
tions are not particularly severe if temperature measurement is made at the centre; in the 
writers’ work, for example, a length/diameter ratio of about 20 is expected to reduce the 
error from this cause well below the experimental error. 

In order to determine thermal constants, two experiments will be made. From the 
results of a ‘large time’ experiment a plot of @, against log t will enable the values of K, 
and (log h.2+2K,/bH) to be determined to a first approximation. From the results of a 
“small time’ experiment, H can be determined roughly by rectification of the initial curve. 
Then a first approximation to /,? can be obtained and, if necessary, the original value of H 
improved using further terms of the small-time series. 

If 2K,/bH is sufficiently small, then it is obvious that this term can be neglected and the 
small-time solution will not be required. 

Experimental work using these results is in its early stages, and a full report on both 
theoretical and practical aspects will be made later. 


Department of Physics, J. H. BLackwELt. 
University of Western Ontario, A. D. MISENER. 
London, Ontario, Canada. 
25th September 1951. 


CarsLaw, H. S., and Jazcer, J. C., 1947 a, Conduction of Heat in Solids (Oxford: University 
Press) ; 1047 b, tes Methods in Applied Mathematics, 2nd edition (Oxford : Univer- 
sity Press), Ch. XII 

GOLDSTEIN, S., 1932, Proc. fond. Math. Soc., 34, 51- 

Hooper, F. C., and Leppsr, F. R., 1950, Amer. Soc. Heat. Vent. Eng. F., August, 129. 
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Observations on the Interactions of y.-Mesons Underground 


During the course of experiments on penetrating bursts carried out in Holborn 
Underground Station, at a depth equivalent to 60 metres of water, the results of which are 
described elsewhere (George and Trent 1951), we have made some observations on the 
behaviour of penetrating particles traversing large thicknesses of lead. 

The apparatus, shown in the Figure, consisted of a number of counter trays separated 
by lead screens of various thicknesses. The events recorded were coincidences between 
one of the groups of tray T (a, b, c, d, e in figure), the centre counter of tray M, and one 
or more of the counters of tray B. In order to record which of the groups of counters of T 
and which of the counters of B were discharged in any coincidence, these counters and 
counter groups were connected to a 25-channel hodoscope of orthodox design. 


A 


Our results may be separated into two groups, firstly those in which a single counter 
or counter group was discharged in each bank, and secondly, events in which the discharge 
of a single counter group in T and the centre counter in M was accompanied by the 
discharge of two or more counters of B. 

In the first case the large majority of coincidences could be explained in terms of the 
passage’ of a single, undeviated particle through the apparatus. A small proportion of 
cases, however, were such that explanation in terms of a single, undeflected particle was not 
possible. These events could be explained either as a result of the production by a 
penetrating particle of one or more knock-on electrons in the centre portion of the apparatus, 
or else by the penetrating particle suffering large angle scatter. 

From considerations of geometry, and from some subsidiary experiments on the 
production of knock-on electrons, we conclude that a large proportion of these events may 
be due to scatter, and further, the angular distribution of this scatter is such that it could 
not be due to coulomb interaction alone. We therefore conclude that nuclear scatter is 
occurring, and assuming that the particles observed are u-mesons, the cross section deduced 
is 2x 10-*8 cm” per nucleon. On the other hand, the results may be explained by the 
presence of a small proportion (~1°%) of particles with nuclear geometric collision 
cross section, or by deviations of the w-meson spectrum from the form assumed for the 
purpose of calculating the contribution due to coulomb scatter, viz. N(E)dE=dE/(E+E,)', 
where Ey equals 1:2 x 101° ey. 

This result is not strictly comparable with the results obtained at sea level by various 
other workers, and in any case, with the exception of the work of Amaldi and Fidecaro 
(1951), most other results on the scattering of mesons appear to have been carried out by 
means of cloud chambers, with consequent very different conditions of geometry, etc. 
Our results do not compare well with those of Amaldi and Fidecaro, who obtained a 
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considerably smaller cross section in their experiments with an apparatus similar to ours. 
Nevertheless, the figure of 2x 10-®8 cm? per nucleon appears reasonable in view of the 
values obtained by cloud-chamber workers (Sahiar, private communication). 

With regard to the events in which two or more counters were discharged in bank B, 
‘simultaneously with a single counter in banks T and M, our conclusions remain 
substantially as suggested previously (George and Trent 1949), i.e. that groups, primarily 
pairs, of penetrating particles are produced with a cross section of 5 x 10-9 cm? per nucleon. 
We have, however, been able to ascertain more accurately the contribution due to knock-on 
electrons produced in the lead A, by introducing lead slats between the counters of bank B 
and noting the reduction in the number of coincidences with two or more counters 
discharged in B. From these later results we conclude that these events may not all be 
explained in terms of knock-on electrons, but that pairs of penetrating particles are 
produced, with a cross section having a lower limit of about 4 x 10-®° cm? per nucleon. 

Further, as a result of the absence of a transition effect when an additional 30 cm. of 
lead was inserted at A in the telescope, the tray T being moved to position 'T’ for this 
purpose, we conclude that at least one of the penetrating particles is either s‘rongly absorbed, 
as would be a 7z-meson, or else is generally of low energy (about 200 to 300 mev.). 

The results of further work on these phenomena are in process of publication. 


Birkbeck College, FE. P. GEORGE. 


London. PEs ERENT. 
2nd October 1951. 


AMaALDI, E., and Frpecaro, G., 1951, Phys. Rev., 81, 339. 
GerorcE, E. P., and TRENT, P. T., 1949, Nature, Lond., 164, 838; 1951, Proc. Phys. Soc. A, 64, 733. 


Nuclear Spin of 3{}Chromium 


Previous observations (Bleaney 1950) of the paramagnetic resonance spectrum of the 
chrome alums have failed to reveal any hyperfine structure. This is due partly to the 
low abundance (9:4°%) of the single odd isotope of chromium, mass 53, and partly to its 
presumably small nuclear moment. By using a crystal of potasstum chromium selenate 
alum, diluted with the corresponding aluminium alum, and grown from heavy water, 
we have managed to reduce the line width sufficiently so that a small satellite line on either 
side of the main line due to the even isotopes could be observed. If this is ascribed to 


Hyperfine structure of 33Cr in the electronic transition M=}<—> —} of 


dilute potassium chrome selenate alum. 


the expected hyperfine structure of *?Cr, the overall splitting (= 50 gauss) is of the expected 
order of magnitude, but the nuclear spin could not be determined owing to intermediate 
lines being obscured by the strong line of the even isotopes. Relative intensity measurements 
suggested that the spin is 3/2, but the accuracy was insufficient to be certain. 

These measurements have been repeated on a sample of chromium enriched in the 
isotope 53. Four lines are clearly visible in the photograph of the spectrum shown in the 
figure, confirming the value 3/2 for the nuclear spin. Fitting the spectrum to the usual 


Hamiltonian H =28H.S+ D{S2—S(S+1)}+A S.1 
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with S=3/2, the values of the constants are found to be 
g=1-976+ 0-002, |2D| =0-21 cm-1, |A|=0-0018;+ 0-0001 cm-. 


In a crystal containing ordinary water instead of heavy water the value of |2D]| was 
0:19 cm~—!, while g was unaltered. 2D is the electronic splitting in zero magnetic field. 

An estimate of the nuclear magnetic moment of ®*Cr can be obtained by comparison 
with the spectrum of the isoelectronic ion V++ (Bleaney, Ingram and Scovil 1951). In 
each case the configuration is 3s?3p°3d?, 4F, and Abragam and Pryce (1951) have shown > 
that hyperfine splitting should be zero except through admixture of the state 3s3p®3d*4s 
by configurational interaction. Quantitative calculation of this admixture is not possible 
from first principles, and one writes A = —2y{fx (1/r*)k, where « is an empirical coefficient. 
It turns out that the values of the product {(1/r?)x} vary by not more than 20° for the 
four ions V++, Mn++, Cott and Cu*+. Physically this means that the contribution to 
the magnetic field at the nucleus due to configurational interaction is constant to this 
accuracy. If this holds also for Cr++*, one obtains the value 0:-45+0-1n.m. from the 
ratio of the hyperfine structures of Cr+++ and V**, using the known magnetic moment 
of *1V (Knight and Cohen 1949). The empirical data which justify this method are 
somewhat scanty, and the result should be accepted with some reserve. The sign is not 
determined, but is expected to be negative (Schawlow and Townes (1951) predict 
—0-:67 n.m.). 

The magnitude of the hyperfine structure in zero magnetic field is of interest in work 
on adiab-tic demagnetization below 0:01° k. Since the crystalline electric field splits the 
electronic levels (+ 3/2) away from the levels (+ 1/2) by about 0:3° K., they must be treated 
separately. The hyperfine structure splits the former into four equally spaced doublets 
with an overall separation of 0:011,°«. The splitting of the (+1/2) levels is more 
complicated. There are two singlets at —0-0060 and +0-0046, and three doublets at 
—0:0055, +0-0041 and +0-0020, these values being in degrees Kelvin relative to the 
position of the unsplit line. 

The authors wish to record their gratitude to Dr. P. F. D. Shaw, who converted the 
chromium to the desired chemical form, and to the Director and the Electromagnetic 
Separator Group of the Atomic Energy Research Establishment, Harwell, who supplied 
the enhanced chromium 53. Our measurements show that the sample provided contained 
at least 9794, of the isotope 53. 


The Clarendon Laboratory, B. BLEANEY. 
Oxford. K. D. Bowers. 
6th October 1951. 


ApracaM, A., and Pryce, M. H. L., 1951, Proc. Roy. Soc. A, 205, 135. 

BLEANEY, B., 1950, Proc. Roy. Soc. A, 204, 203. 

BLEANEY, B., INGRAM, D. J. E., and Scovit, H. E. D., 1951, Proc. Phys Soc. A, 64, Gor. 
Knicut, W. D., and CoHEN, V. W., 1949, Phys. Rev., 76, 1421. 

ScHAwLow, A. L., and Townss, C. H., 1951, Phys. Rev., 82 268. 


Interference in the Directional Correlation of 
Conversion Electrons 


Two previous communications (Gardner 1949, 1951, denoted hereinafter by I and [I) 
have investigated the directional correlation between the two conversion electrons 
successively ejected from an atom whose nucleus undergoes transitions A>+B—-C in 
cascade. In I both nuclear transitions were assumed to be pure electric multipole in 
character; in IT pure electric-pure magnetic and pure magnetic-pure magnetic cascades 
were considered. It was also remarked in II that, to the approximation imposed by the 
initial assumptions, cascades with either or both transitions a mixture of electric and 
magnetic multipoles could be treated by taking just a weighted mixture of the appropriate 
pure cascades, there being no interference terms. The absence of interference terms was 
also noticed by Fierz (1949) in treating conversion-y correlation in the approximation 
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of I and II. The present note shows that this non-appearance of interference terms in 
c-c and c~y correlations derives essentially from two assumptions, namely, that the 
conversion electrons are ejected from an initial s-state, and that a non-relativistic treatment 
of internal conversion is justified. If either (or both) of these assumptions be disallowed 
the correlation formula will, in general, contain interference terms. 

In I and II the final-state total wave function %, is expressed as a linear combination 
of the final-st te nuclear wave function ¢, (subsequently dropped as a constant) and the 
final-state electron wave functions %,, #2, the coefficients of the various terms in this 
expansion being the familiar products of Wigner (Clebsch—-Gordan) coefficients for the 
vector addition of angular momenta. Now if the transition A+B is a mixture of electric 
and magnetic multipoles 4, breaks down into ‘ electric’ and ‘ magnetic’ parts, #,2 and 
a"; similarly with }, if B+C is a mixture. Hence y%, must be written as a mixture, say 
arbs®+6s,™, so that in using the formula I(6)~ X [*%,dr for the relative angular 

My 
distribution of the two electrons (see I and II) we have the following product to consider : 


(oct ypeP + BAS) (cabs B+ Bb.) =| ayfy® |?+| Bib .™ [2 a* By TPM + oP * yp Byy™. 

The first two terms on the right-hand side are the pure multipole terms, fully discussed 
in I and I]; it is the second two—the cross terms—which would give rise to interference 
effects, but which vanish in the above approximation. To see this one can use the 
expansion of #,¥ as a linear combination of the %,", #8, the expansion of %,™ as a linear 
combination of the %™, 4", and the corresponding complex conjugate expansions for 

<* and 4". On substituting these expansions in f*?,™ and Y,"%*™, and multiplying out, 
the general term in the resulting expression is found to contain a factor of the type : 

LX Chay Cove’ 
tobe 168 

where, in the notation of I and II, the C’s are Wigner vector addition coefficients, a’, b’, c, c’ 
are angular momenta and a’, f’, y, y’ their magnetic components. In particular the 
c, c’ are the orbital angular momenta of an ejected electron associated with a magnetic 
2!-pole and an electric 2!+1-pole respectively. (This is the only type of mixture of practical 
interest.) If c, c’ are different the above sum of products of Wigner coefficients vanishes 
by their familiar orthogonality property (see, for example, Racah 1942); that is to say 
PrPb.™ and W,°s%™" vanish since every term in their expansions vanishes. Now for 
electrons initially in an s-state the non-relativistic theory of internal conversion predicts 
c=l—1 (from the magnetic 2!/-pole) and c’=/-+1 (from the electric 2!+1-pole). In the 
relativistic treatment, however, there are in general two alternative values of orbital angular 
momentum available to an electron associated with a particular (electric or magnetic) 
multipole, according to whether or not ‘ spin-flip’’ occurs. Therefore one has to consider 
two values for c and two values for c’, and as these four values are not all different the 
interference terms do not all vanish. Again, suppose we remove the restriction to 
s-electrons, and consider conversion electrons ejected from an initial p-state, say. Each 
electron then has an initial orbital angular momentum which can in general combine in 
three possible ways with that which it acquires from the nucleus, so, even in the 
non-relativistic treatment, there are three alternative values for c and three for c’; and since 
these are not all different, not all the interference terms vanish. When considering specific 
experimental examples one must, as always, bear in mind the operation of any special 
selection rules which may prohibit a common value of ¢c and c’ and so prevent interference. 

The inadequacy of the non-relativistic theory for calculating magnetic internal 
conversion coefficients, even at relatively low energies, has been pointed out by Drell (1949) 
who gives a semi-relativistic treatment (neglecting higher powers than the first of 
hyv/mc?). In conversion correlation the error is probably not as serious since one has here 
to deal with a ratio of two matrix elements which are affected more or less equally by the 
first-order relativistic correction. Fully relativistic calculations for conversion correlation 
in pure multipole cascades are being carried out by M. E. Rose and coll.borators at the 
Oak Ridge National Laboratory, Tennessee; interference effects are being investigated 
by S. P. Lloyd at the Princeton Institute. I am indebted to Professor Drell, Dr. Rose and 
Dr. Lloyd for correspondence; also to Professor Racah who kindly sent me the manuscript 
of his forthcoming paper (Racah 1951) dealing with the general angular correlation problem 
(all kinds of particles, with or without interference). 
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In conclusion, this opportunity is taken to correct an earlier statement concerning a 
paper by Berestetzky (1948). The treatment of internal conversion in that paper is 
completely non-relativistic, and not, as erroneously reported in II, semi-relativistic in the 
sense of Drell’s calculations. 


Division of Atomic Energy, J. W. GarDNER*. 
National Research Council of Canada, 
Chalk River, Ontario. 
30th August 1951. 
* N.R.L. Postdoctoral Fellow. 
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Theoretical Discussion of the Beta Decay of *P 


The apparent contradiction between the empirically second forbidden ft value of **P 
and its allowed shape beta spectrum has not yet been explained. The beta spectrum has 
been measured by Siegba n (1946), Langer and Price (1949), Agnew (1950) and Warshaw, 
Chen and Appleton (1950), and there can be little doubt that the spectrum has the allowed 
shape for energies greater than about 250 kev. Below this energy, Agnew and Warshaw 
et al. have reported a marked excess of electrons over the number predicted for an allowed 
spectrum; this, however, may be due to a second low energy spectrum with an end point 
at about 250 kev. and an intensity about 1 or 2% of the main transition. 

None of the second forbidden correction factors are capable of explaining the shape 
of the **P spectrum. The only correction factors which do so are those for the allowed 
transitions, for first forbidden transitions with unit change of spin on the tensor and polar 
vector interactions, and for first forbidden transitions with no spin change on the axial 
vector and pseudoscalar interactions. Since no gamma-rays are observed, the transition 
must be to the ground state of *?S, and since this has spin zero and the spin change in the 
transition is not more than 1, the spin of **P must be 0 or 1. Since the transition must be 
either allowed or first forbidden, some mechanism such as an additional selection rule is 
necessary to explain the long lifetime. 

According to the shell model, neutron and proton levels are filled in the order 
S1/2) Ps/2, P1/2> Gs/2) Si/2, ds/2, - . . (this ordering of the s,/. and ds. levels is necessary to 
explain the observed spins of odd nuclei). The shell model interpretation of the decay 
of **P is therefore that a ds;. neutron is transformed into a s,/, proton; this involves no 
parity change, and consequently the transition cannot be first forbidden and must be 
allowed. 

It has been shown (Pursey 1951) that for single particle transitions involving an orbital 
angular momentum change of two units, and a total angular momentum change of one unit 
the matrix elements fo and {fo (which are usually classed as allowed) are of order (v/c)?, 
where v is of the order of the nucleon velocities. Since according to the shell model the 
decay of *?P is of this type, this may account for the long half-life of the transition. The 
correction factor for the transition, assuming the tensor interaction, is the same as Cp 
given by Konopinski and Uhlenbeck (1941) if [a is replaced by [Bo and foxr by Jfaxr. 
If the magnitudes of the matrix elements [fo and {fa xr are assumed to be (v/c)? and up/c 
respectively, where p=(1/274)A}/* is the nuclear radius in units of the Compton 
wavelength of the electron, then the theoretical lifetime depends only on the sign of the 
cross term in the correction factor and on one parameter x=v/cp which may be chosen to 
fit the observed lifetime. Good fits are obtained by taking x=12 (v/c=0-14) if the cross 
term is negative and x=5-5 (v/c=0-064) if the cross term is positive. These values of the 
ratio of u/c to p are rather larger than those usually taken : the writer’s work on the ratio 
of fa to Jr (Pursey 1951) suggests, however, that they are not excessive. The theoretical 
spectra for both values of x used are experimentally indistinguishable from allowed 
spectra. A similar explanation of the transition is, of course, possible using the axial vector 
interaction. 
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Sherwin (1951) has measured the angular correlation between the electron and neutrino 
emitted from *P, and found it to be given approximately by (1+ cos @). From the 
theoretical work of Greuling and Meeks (1951) it is clear that this angular correlation is 
consistent only with allowed transitions with no spin change on the polar vector interaction, 
and with first forbidden transitions with no spin change on the axial vector interaction; 
it is not consistent with the explanation of the spectrum shape and lifetime given above. 
The spectrum shape, lifetime and angular correlation can all be explained by assuming 
‘that the transition is allowed, and that the interaction is mainly tensor or axial vector with 
an admixture of about 1°% of polar vector; this explanation, however, seems very artificial. 

A direct measurement of the spin of °*P would be of great value. If the spin is greater 
than 1, then it is impossible to explain the shape of the spectrum on the present theory of 
beta decay. If it is 1, then the lifetime and the spectrum shape can both be explained 
simply and naturally on the basis of the shell model, but Sherwin’s angular correlation 
measurements are in contradiction with the theory given by Greuling and Meeks. If the 
spin is 0, the angular correlation can be explained, but the shell model level scheme breaks 
down, and the problem of reconciling the lifetime and spectrum shape remains unsolved 
unless the interaction is assumed to contain an admixture of the polar vector interaction 
so small as to be unobservable in any other experiment. 

The writer wishes to express his gratitude to the Department of Scientific and 
Industrial Research for financial assistance during the course of this work. 


Department of Natural Philosophy, D. L. Pursey.* 
University of Glasgow. 
2nd October 1950. 


* Now at Physics Department, King’s College, London. 
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The G.L.E.E.P. as an Absolute Standard of Thermal 
Neutron Density 


The Graphite Low Energy Experimental Pile at Harwell has been calibrated for use as 
a standard of thermal neutron density. It can be used within the range 10 neutrons/cm® 
to 2x 104 neutrons/cm?, the density being known throughout the range to +5%. We 
have considered as therma! neutrons those with energies below the cadmium ‘cut-off’, 
i.e. below about 0-4 ev. 

We have made an absolute measurement of the density of thermal neutrons at the 
centre of the pile in terms of the current from an ion chamber permanently situated 
1 metre from the centre. The pile is normally operated only at low powers, with small 
temperature variations, and with a constant amount of absorbing material in the pile. Thus 
the neutron density distribution remains constant with time, and the chamber current 
provides a measure of density at the centre. 

The ion chamber is of the boron-coated type (Rossi and Staub 1949), and the ion current 
from the chamber flows through a wire-wound 100 MQ resistor, kept in a temperature- 
controlled cabinet. Resistances of down to 5 MQ can be selected from the resistor. The 
voltage produced across the resistor is balanced using a potentiometer, any out of balance 
being observed by means of a directly coupled amplifier. The potentiometer reads from 
1 to 100 volts with an accuracy better than 0:01 volt. Chamber currents between 0:01 and 
20 microamperes can therefore be measured, corresponding to the density range quoted 
previously. 

Now if o is the absorption cross section of an element for neutrons of velocity v, and 
these neutrons have a density p, then the total number of thermal neutrons absorbed per 
second by the element is [{pvo dv, where the integral extends over the whole thermal 
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velocity range. If the element follows a 1/v cross-section law in the thermal region, then 
o=Gpty/v, where oo is the cross section at some standard velocity vp. ‘Thus the total rate 
of absorption of thermal neutrons by the element is o9v% | p dv. 

It is apparent, therefore, that an absolute determination of the total thermal neutron 
density in the pile can be made by a measurement of the neutron absorption cross section 
of an element and a determination of its rate of absorption of thermal neutrons, provided 
the element follows a 1/v law in the thermal region. 

The element chosen for the absolute density determination was sodium. This has 
only one stable isotope, and absorption of thermal neutrons leads always‘to the 14-9 hr. 
activity of #4Na. In addition, sodium follows a 1/v law in the thermal region, and a 
measurement of the absolute f-activity acquired by a sample of sodium was used as a 
measure of the rate of absorption of thermal neutrons. ‘The sodium was used in the form 
of carbonate, neutron absorption by carbon and oxygen being negligible. A measurement 
of the activity acquired by a similar sample irradiated inside a cadmium box enabled us 
to allow for the effect of neutrons with greater than thermal velocities, both in the 
determination of the rate of absorption of thermal neutrons and in the measurement of 
cross section. 

The determination of the absolute f-activity of the sodium was made using a liquid 
counter which had previously been calibrated for 74Na. The calibration was made using 
a sodium sample whose absolute f-activity had been determined by both a defined solid 
angle and a coincidence method (Putman 1950). The two methods gave calibrations which 
agreed to better than 1%. 

The absorption cross section of sodiurn at 2,200 m/sec. was obtained by comparison 
with that of boron using a swing method (Colmer and Littler 1950), a value of 710 barns 
at 2,200 m/sec. being assumed for boron. 

Another measurement of absolute thermal neutron density has been made by a 
determination of the absolute B-activity of a sample of manganese irradiated in the B.E.P.O. 
(Price 1951), combined with a measurement in the G.L.E.E.P. of the neutron absorption 
cross section of manganese at 2 200 m/sec. The manganese used was in the form of 
sulphate, and the f-activity was determined by a defined solid angle method. 

The conditions under which the B.E.P.O. is used are such that the neutron density 
distribution varies appreciably with time. Instead of expressing the density in terms 
of ion chamber current, therefore, a manganese foil was calibrated in terms of absolute 
thermal neutron density. This foil was irradiated in the B.E.P.O. under conditions known 
with respect to those used for the irradiation of the mangan:se sulphate, and the activity 
of the foil was measured on a Geiger—Miiller counter, using a standard uranium foil for 
checking the counter sensitivity. Other manganese foils and Geiger—-Miuller counters were 
then standardized for measurement of absolute thermal neutron density. 

A comparison of the sodium and manganese density measurements was made by 
irradiating some of the standardized manganese foils in the G.L.E.E.P. The absolute 
density of thermal neutrons determined by the two methods agreed to better than 2%. 

The neutron density is checked periodically by measuring the radioactivity of other 
manganese foils irradiated in the centre of the G.L.E.E.P. against that of a radium standard. 
No significant change in the density has been observed during the past eighteen months. 

In addition to making the density measurement, we have used the method described 
by Anderson et al. (1947) to determine the mean velocity of the thermal neutrons in the 
G.L.E.E.P. This was found to be about 2,690 m/sec. at a pile temperature of 300° x. 

We would like to express our thanks to Miss A. M. Wildblood, who carried out the 
B-counting of ?4Na for us. 


Atomic Energy Research Establishment, D. J. cItrrER: 
Harwell, Didcot, Berks. Bab pwockErn 
12th September 1951. Bae sERICE: 
fo, H. L., Fermi, E., WATTENBERG, A., WEIL, G. L., and Zinn, W. H., 1947, Phys. Rev., 
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Obituary Notices 


LESLIE JOHN COMRIE 


“THE death of Dr. L. J. Comrie, F.R.S., in London on 11th December 1950 at the age of 
fifty-seven has removed an original and independent figure from the ranks of British 
Scientists. 

Leslie John Comrie was born on 15th August 1893 at Pukekohe, Auckland, New 
Zealand. He was educated at Auckland Grammar School and Auckland University College, 
where, after a promising career, he graduated M.A. in 1916, with First Class Honours in 
Chemistry. He served in the First World War with the New Zealand Expeditionary Force, 
incurring severe wounds, with consequent amputation of a leg. It was doubtless this 
incapacitation, together with the handicap of increasing deafness, that caused him to recon- 
‘sider his career, for after the war he remained in Britain as a New Zealand Expeditionary 
Force Scholar, studying astronomy and computation at St. John’s College, Cambridge, where 
he was an Isaac Newton Student from 1920 to 1923. In 1923 a thesis, devoted to the history 
cand. detailed examination of methods of predicting the occultation of stars by planets, 
gained him the degree of Ph.D. 

Two years in America followed, the first as Assistant Professor of Mathematics and 
Astronomy at Swarthmore College, Philadelphia, the second as Assistant Professor of 
Astronomy at Northwestern University, Evanston, [llinois. In 1925 he was appointed 
Deputy Superintendent of H.M. Nautical Almanac Office, Greenwich, and in 1926 entered 
upon his duties there. He was no stranger to Greenwich, having already, in 1920, spent 
two months at the Royal Observatory as part of the conditions of his Isaac Newton student- 
ship. In 1930 he became Superintendent, a post which he held until 1936. Without 
-delay he used the opportunity to embody in the next issue of the Nautical Almanac, for 1931, 
certain changes which he had without doubt considered for some years, notably the use of a 
‘standard equinox. The Almanac of 1931, breaking away from the tradition of almost a 
hundred years, incorporated, among other radical changes, the new ephemerides. 

Comrie’s imagination, essentially direct and pragmatic, had been captured, from his 
first contact with calculating machines, by their unexplored possibilities. He saw that with 
the rapid improvement in calculating machines many accepted procedures, legacies from the 
-pre-machine era, were wasteful in time and effort and should be superseded or re-drafted in 
the newer idiom. He reconsidered such processes as interpolation, numerical integration and 
-subtabulation from first principles, and devised simple and powerful repetitive techniques, 
requiring the minimum of copying and capable of being mastered and understood by 
operators of only moderate mathematical training. His ingenious mind ever busied itself 
with inventing methods for turning the latest devices, such as transfer levers, split keyboards 
or twin machines, to practical use. It was characteristic of him, for example, to harness 
together National Accounting Machines in order to build up and print tabular values from 
differences of high order, or to use punched cards and the Hollerith machine, as he did in 
1928 with E. W. Brown’s Tables of the Moon, to sum large numbers of harmonic terms. His 
methods, severely simple and simply described are explicit in the prefaces or introductions 
to the various tables published under his editorship, notably in Interpolation and Allied 
Tables, reprinted from the Nautical Almanac for 1937. 

In 1936 he resigned from the Superintendentship and began a new and independent 
venture, the founding on business lines of a computational firm, the now celebrated Scientific 
Computing Service Limited. This, at the time, was a bold step, though soon justified by 
4ncreasing success, by requests for computational help on a large scale by the universities, 
by industry and by Government departments. With the outbreak of the Second World 
War the demand, official and unofficial, for such services was extraordinarily accelerated. 
‘The office in Bedford Square served as training ground for many computers, who entered 
with moderate qualifications and passed out as experts. By 1944 the staff numbered thirty; 
-and after the war, with the business now firmly established and enlarged, Comrie found the 
need of an assistant manager, of similar interest and authority in all matters concerning 
computation. This was Dr. J. C. P. Miller, who now succeeds him. 

Comrie had now for many years been recognized as an authority of international 
reputation on mathematical tables. He had been the first Director of the Computing 
Section of the British Astronomical Association, Secretary to the Mathematical Tables 
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Committee of the British Association, 1930-36, member of the Committee on Mathematical 
Tables of the National Research Council, Washington; and from 1948 to 1950 he was a 
member of the Mathematical Tables Committee of the Royal Society. He was elected a 
Fellow of the Royal Astronomical Society as early as 1919, and served on its Council 
from 1929 to 1932. He was also a Fellow of the Royal Statistical Society, and of the Physical 
Society. He was an active member of the International Astronomical Union, attending all its. 
conferences from 1928 onwards, and was elected President, in 1932, of its Commission on 
Ephemerides. His native country recognized him by honorary fellowships, of the Astro- 
nomical Society of New Zealand in 1933 and of the Royal Society of New Zealand in 1945. 
In March 1950 he was elected a Fellow of the Royal Society, but lived only ten months to- 
enjoy this culminating honour. 

His papers, though they number about a hundred, represent only a small part of his 
work. A proportion of these, mostly appearing in the Monthly Notices of the Royal 
Astronomical Society, are on astronomical matters, but the larger proportion (often appearing 
in such journals as the Office Machinery Users’ Association Transactions) deal with table- 
making and the application of various types of machine to problems of computation. Of 
great importance are his reviews and notes on errata in existing mathematical tables, a 
characteristic example being the article in Mathematical Tables and Aids to Computation, 1 
(1945), pp. 391-399, on Jahnke and Emde’s Tables of Functions. His erudition in this field 
was immense. 

The tables for which he was personally responsible bear on every page the marks of his 
exacting standards, imposed not least upon himself. First, by rigorous systematic 
checking and counterchecking, the values must be impeccably accurate; they must also be 
interpolable with the maximum of convenience, readable with the maximum of ease, and 
typographically elegant. No one can doubt, perusing, for example, the third and fourth 
enlarged editions of Barlow’s Tables, issued in 1930 and 1941, or, perhaps most con- 
spicuously, the new Six-Figure Mathematical Tables in two volumes, logarithmic and natural,. 
published by Chambers in 1948, that such aims have been achieved to the full. These last 
tables call for special mention; they resulted from the decision of the publishing firm, on the 
centenary of the old tables in 1944, to take the advice of a modern expert in regard to their 
re-issue. ‘They are rightly described in the Editor’s preface as “‘ one of the greatest of all 
British table-making enterprises ’’. 

In character Comrie was intolerant of imprecision, direct and outspoken to a degree 
sometimes startling, though this was really a concomitant of his essential simplicity and 
singlemindedness. The writer of the present notice was never a fellow committee-man with 
him; but it is no secret that his abruptness of manner in delivering an opinion, his refusal to: 
use any of the ordinary arts of compromise, did not always meet with a favourable response. 
This was, however, only one side of the man, reflecting his inveterate honesty. 

Perhaps one who was bound to Comrie by ties of friendship, of birth in New Zealand at 
about the same period, and of similar service and similar though less severe experience im 
the First World War, may conclude on a specially personal note. Comrie as a young man, 
before departing for France and suffering so severely, was physically vigorous and energetic, 
a good rifle shot and tennis player, pursuits in which his qualities of precision and decisive- 
ness found a ready outlet. It is more than probable that the disabilities that permanently 
precluded such active outlets produced a certain frustration, which in its turn was overcome 
by the unremitting energy which Comrie poured into everything that concerned astronomy, 
machines, tables and computation; for, as has already been stated, the amount of work he 
did in his lifetime cannot by any means be summed up in terms of publications. Not least 
was the help he unobtrusively and willingly gave to others in the field, of which a fitting 
memorial exisis in the Index of Mathematical Tables of Fletcher, Miller and Rosenhead, 1945, 
already classical. New Zealanders visiting London, and astronomers and computers from 
America and Europe can testify to the warm hospitality he gave them, ably seconded by his. 
wife; but very few knew of his many acts of secret and private generosity, some of which have 
only come to light since his death. The picture that lingers most in the mind is of Comrie as 
host, cut off by deafness, it is true, from the cross-currents of conversation, but himself 
sitting quietly happy at the evident enjoyment of his guests. It is a question whether, in 
giving full rein to his proclivities, he did not drive himself too hard along a single track; but 
his work is of great, permanent and visible value, and he has set standards, in all that concerns 
table-making and computation, that can hardly be surpassed. 


A. C. AITKEN. 
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LEONARD GEORGE GRIMMETT 


I WORKED as assistant to L. G. Grimmett for 18 months, 1934-35, in the Physics: 
Department of the Radium Beam Therapy Research at The Radium Institute, London. 
I have had as transient a contact with many others, yet they have passed out of my 
recollection and interest. In contrast, Grimmett made a lasting impression. I have always 
felt a surge of pleasure on meeting him again, or on hearing of his exploits, and was very 
grieved to learn of his sudden death from heart failure on 27th May 1951 in Houston, 
Texas. He had a unique and admirable personality. 

The two pillars of Grimmett’s career were his love of science and of music; he was a 
pianist of the first rank. On leaving the Holloway County Secondary School he paid the 
expenses of a physics degree course at King’s College from his earnings in restaurant and 
theatre bands. He made a voyage as ship’s pianist to Buenos Aires. Undoubtedly this 
double life was a heavy strain. To me he deplored that, at a time when he was engaged 
in this struggle, his Professor was completely discouraging, believing that a degree could 
not be taken in this eccentric fashion. Nevertheless it was taken (in 1926), and so well that 
he was encouraged to do three years’ research. In 1930 he became Assistant Physicist 
at the Westminster Hospital Radium Annexe, working under Professor H. T. Flint, and 
together they developed one of the earliest hospital physics departments. Grimmett was 
a pioneer in the development of radium beam units for the treatment of malignant disease; a 
two-gramme unit designed by him in the early 1930’s is still in use. Methods of auxiliary 
physical measurements, such as those of dose distribution round these units, which are 
now routine, also claimed his attention. 

He collaborated with Dr. F. G. Spear (1932) in an early quantitative study of the 
influence of gamma-ray intensity on the inhibition of cell division in tissue cultures (Spear 
and Grimmett 1933). This work caught the eye of Professor McLennan, Chairman of the 
Executive Committee of the Radium Beam Therapy Research, who persuaded him to 
join that body. He served it, and the Radiotherapeutic Research Unit of the Medical 
Research Council which assimilated it, until 1944. 

Grimmett excelled in the design of apparatus; he was an artist in temperament. The 
common view that if a scientific apparatus works it is satisfactory was foreign to him.. 
It must not only work but be well designed, elegant, and properly finished. Furthermore, 
each stage in its development must reach the same standard. It was not sufficient that 
the drawing should be adequate for the mechanic’s use; the penmanship must be excellent 
also. He saw everything in a big way. The idea of improving the collimation of a radium 
beam by fitting a solid gold nozzle to the teleradium unit was typical, and he made it a 
reality (Flint, Grimmett, Rock Carling and Cade 1934). He spent time and money 
without inhibition or fear to get the results he wanted. He supplied most of the ingenious. 
ideas for the 5-gramme radium unit (later 10 grammes), with pneumatic transference of 
the radium (Grimmett 1937), replicas of which have been installed in hospitals all over the: 
world. He was one of the first to realize the importance of the total absorption of radiation 
by a patient (now known as ‘ volume dose’), as distinct from the absorption in the 
malignant growth itself. He tackled its measurement practically by the construction of a 
full-size human shape of parallel celluloid sheets with narrow air gaps which formed. 
ionization chambers (Grimmett 1942). 

He believed in the need for higher voltages in x-ray therapy, and in 1939 began the: 
design of what ultimately became the first high pressure Van de Graaff generator to operate 
in Britain. Harassed by the war, the work was moved from the bombed Radium Institute 
to the Imperial College of Science and then to the Hammersmith Hospital, advancing 
slowly, since each stage must reach perfection without short cuts. Concurrently he 
supervised the physical side of a long clinical experiment which compared the relative 
merits of radium beam and 200 kv. x-ray treatments, with the fields of both types arranged. 
to have closely matched dose distributions. At Hammersmith he built up a fine team 
and a magnificently equipped physical laboratory. 

Then in 1944 he disappeared from his usual scenes, reappearing briefly from time to 
time from Paris, or wherever else his service for the British Council, UNESCO, and UNO 
had taken him. Finally it was reported that he had settled in Houston, Texas, where he 
was developing a radiobiological research laboratory in the really big way which suited. 
his temperament. 
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Grimmett’s career might suggest the dominant, aggressive personality, yet he was 
kind, gentle, always soft-spoken, and quite imperturbable except in sympathy for others. 
As we walked down Regent Street in 1934 we passed seedy, out-at-elbow violinists playing 
in the gutter. He was filled with distress—‘‘ I know this kind of man; I have played with 
them. They are good musicians displaced by the ‘Talkies’. He declined to assist in 
the atom bomb development. ‘‘ The M.R.C. were approached’, he said, ‘‘ but I asked them 
to squash the idea”; then, with a wry grin in oblique explanation, “I don’t mind killing 
Germans in odd ones ’”’. 

Grimmett showed that, with few initial advantages other than talent, courage and 
tenacity, he could be a ‘ 20th century ’ adventurer in the best sense, and at the same time 
make substantial contributions to human welfare. All who knew him will miss him 
greatly, and feel deep sympathy with his widow. 

JOHN READ. 
Firnt, H. T., Grimmett, L. G., Rock CARLING, F., and CabE, S., 1934, Brit. Med. F., 14 April. 


Grimmett, L. G., 1937, Brit. J. Radiol., 10, 105; Ibid., 1942, 15, 144. 
Spear, F. G., and Griwme_ett, L. G., 1933, Brit. J. Radiol., 6, 387. 


JOHN R. LOOFBOUROW 


Dr. LoorBouRow, who had been a Fellow of the Physical Society since 1941, died in 
Boston, Massachusetts, on 22nd January 1951. At the time of his death he was a professor 
of biophysics at the Massachusetts Institute of Technology, where he was in addition 
Executive Officer of the Biology Department and Chairman of the Faculty. 

Loofbourow, the son of John Wilson and Henrietta Loofbourow, was born in Cincinnati, 
Ohio, on 1st November 1902. He was educated at the University of Cincinnati, where he 
obtained the degree of Bachelor of Arts in 1923. He served as an instructor in Physics there 
from 1925 to 1929, and as a research associate for the following six years. In 1936 the 
University of Dayton in Ohio awarded him the honorary degree of Doctor of Science, and 
he became professor of Biology at that institution. From 1935 to 1940 he was a research 
professor at the Institution Divi Thomae in Cincinnati, and in 1940 he was called to the 
Massachusetts Institute of Technology, where he was made an associate professor of bio- 
physics. 

During the war years Dr. Loofbourow was Executive Secretary of the Radar Division of 
the National Research Committee. He served as a special adviser to the Atomic Energy 
Commission in 1947 and 1948, in which year he was awarded the Presidential Certificate of 
Merit. 

Dr. Loofbourow was a Fellow of the American Association for the Advancement of 
Science, the American Physical Society, the American Academy of Sciences, the New York 
Academy of Sciences, the Optical Society of America, the Beta Theta Pi. He was also a 
member of the Chemical Society, the Faraday Society and the Biochemical Society of 
England, as well as the Physical Society. 

Dr. Loofbourow’s untimely death from a heart complication was tragically followed by 
that of his wife, a practising physician, from the same cause within ten days. They are 
survived by a son, John W. Loofbourow. 

Professor Loofbourow was a frequent contributor to the literature of biophysics, special- 
izing in the field of absorption spectroscopy, particularly in the ultra-violet region. At the 
time of his death he was engaged in important work on the purines and pyrimidines, and was 
engaged in developing new techniques for absorption studies at extremely low temperatures. 
With, Professor Richard C. Lord and the undersigned, he was the author of the text Practical 
Spectroscopy published in 1948. 

With his unassuming and friendly manner, breadth of interest, and complete grasp of his 
subject, Professor Loofbourow had won a host of friends among both colleagues and students. 
‘As a member of a committee on improving educational methods and as Chairman of the 
Faculty he had contributed much to the academic community. His tolerance and quiet 
competence left a deep impression on all who came in contact with him, and contributed 
greatly to the present feeling of loss among his colleagues, students, and many friends. 


GEORGE R. HARRISON. 
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EDWARD MALLETT 


EDWARD MALLETT, born in 1888, was educated at Southend High School and University 
College, London. Heserved an apprenticeship with Messrs. Crompton & Co., Chelmsford,,. 
and entered the Post Office Engineering Department through a Civil Service competition. 
in 1909. Following a short period of training, he was posted to the West Section of the 
Scotland District and remained there until the outbreak of war in 1914. He was interested 
in the Territorial Army and as a result, as soon as war was declared, was posted to the 52nd 
Lowland Divisional Signals and saw service in the Middle East area of hostilities, finally 
leaving the Service with the rank of Captain. 

After the war, Mallett was posted to the Transmission Section of the Engineer-in-Chief’s. 
Office, and here his grasp of electrical fundamentals combined with his abounding energy 
found outlet in dealing with the problems of telephone transmission, which was then entering 
into a new phase due to the introduction of the thermionic valve. There is no doubt that 
Mallett, had he remained in the Post Office service, would have had a brilliant future, but a 
friendship made during the war years turned his energies in another direction and he left 
the Post Office service in 1920. It was during this period of service in the Post Office that 
he carried out his researches into the vibrational modes of a circular diaphragm which led 
to the work for which he was eventually granted his doctorate. 

In September 1921 he joined the academic staff of the City and Guilds College, South 
Kensington, as Assistant Professor of Electrical Engineering, which post he held until 
September 1932. 

During these eleven years he was developing and expanding the ‘light current ” 
engineering section of the department. A considerable expansion was taking: place under 
the direction of Professor Fortescue, and Mallett, with his active and fertile mind, backed by 
his considerable experience in the Post Office Engineering Department, found scope for 
developing the instruction and research in electrical communication. 

He had an intuitive engineering sense and considerable academic skill, and these com- 
bined produced a lecturer of inimitable style and a research worker capable of producing 
results from the minimum of apparatus. He was extremely popular with both staff and 
students and joined in many of their social activities. 

The title of Reader was conferred on him in 1927 by the University of London, and from 
that date until his death he remained a member of the Board of Studies in Electrical 
Engineering. He was for several years Secretary to the Board of Studies, acquiring an- 
intimate knowledge of the activities of the Board and of the University, which knowledge 
served them well over the period of his long membership. 

In addition to these activities he still found time to write several papers on a variety of 
subjects of interest to the communications engineer which he presented to the Institution 
of Electrical Engineers and to the Physical Society. 

In July 1926 he was awarded the Degree of D.Sc.(Eng.) as an Internal Student of 
University College, London, the thesis title being Forced Oscillations, Electrical and Mechan- 
ical. We was frequently consulted by outside bodies and the University on electrical 
engineering matters, and in 1929 he produced his textbook on Telegraphy and Telephony 


which became, and still is, a valuable foundation for students of electrical communication. 


It was at the height of his academic activities that he left the City and Guilds College to. 
become Principal of Woolwich Polytechnic but, although he became more and more 
immersed in administrative studies, he never lost contact with the College or the University 
of London. 

As Principal of Woolwich Polytechnic he was intensely interested in the University 
aspect of the Polytechnic but at the same time took a very active interest in the large non- 
university side of Polytechnic activities. He was a keen believer in the National Certificate 
Course for part-time and evening students, and had a place in his heart for the craftsman, 
ensuring that he was provided with proper shops and tools. It was during his Principalship 
that the Sandwich Course in Engineering was started, and this eventually developed into the 
Engineering Degree Course of today. This was soon followed by Degree Courses for 
B.Sc. (Chemistry, Physics, Mathematics), and the full-time School of Art Course. 

The Polytechnic Union came into being during his time. He saw that the union of the 
diversified clubs and societies of the Polytechnic into one corporate whole could be of great 
value to the students and members, and the passing of time has proved him correct. The 
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success of the Union in all its aspects is a tribute to him: he gave it a free, unfettered 
hand, and held it in high esteem. 

Extensions to the Polytechnic were made during his time, and future developments 
were being planned by him when his fatal illness intervened. 

He was a member of the Institution of Civil Engineers, the Institution of Electrical 
Engineers, the Institute of Physics and the Physical Society (since 1921), and a Fellow of 
‘University College. 

Edward Mallett, D.Sc.Eng. (Lond.), M.Inst.C.E., M.IE.E., F.Inst.P., Fellow of 
University College, London, died on the 10th December 1950. C. R. URWIN. 


HANTARO NAGAOKA 


PROFESSOR HANTARO NAGAOKA, who had for nearly forty years occupied a worthy 
place emong the few Honorary Fellows of the Physical Society, died on 10th December 1950 
at the age of 85. In the English-speaking world his name is probably most familiar in 
connection with ‘ Nagaoka’s formula’ for the calculation of inductance, and indeed his 
work on the subject has been of great value to workers in the radio field and to physicists 
concerned with precision measurements of inductance. For example, in addition to his 
formulae he published, with S. Sakurai, tables which give, with very little computation, 
values for the self-inductance of single layer coils, and for the mutual inductance of 
coaxial circles to six-figure accuracy. "This work, however, is hardly representative of his 
main interest; he was essentially an all-round physicist in the Rayleigh, Kelvin, 
J. J. Thomson tradition, and one who could refer to physics as ‘the Study of Nature’. _ 

He was born in 1865 and was of honourable lineage. He graduated in physics at the 
‘Tokyo Imperial University in 1887, a time when, as he recalled later, the library in the 
physics department occupied but a single shelf six feet wide, and the only periodical 
available was the Philosophical Magazine. Nevertheless he must have immediately 
embarked upon experimental research in magnetism, for papers published by him in 1888, 
1889 and 1890 are discussed in some detail by Ewing in his Magnetic Induction in Iron 
and Other Metals. Among other things Nagaoka showed that in certain circumstances 
the magnetization of nickel is reversed by torsion, an effect so surprising that both Kelvin 
(then W. Thomson) and Weiss had it confirmed in their own laboratories. Sixty years 
later Nagaoka remarked that the reason why still remained quite unknown. He continued 
his work in Tokyo for six years, took his doctorate in 1893, and then proceeded to 
Germany, where he studied for three years, attending lectures by Helmholtz, Boltzmann, 
Planck and others, and working on magnetostriction. During this time he visited the 
Physikalisch-Technische Reichsanstalt and was received by Lindeck who, having recently 
translated Ewing’s book into German, knew of his work. Lindeck asked him when he had 
left Ewing’s laboratory at Cambridge, and when he learnt that the work had been done 
several years earlier in Japan was sufficiently impressed to show him the laboratory in 
detail. Nagaoka with a pleasant touch of humour said later that he realized then that 
a bit of good work is a much better introduction than a distinguished visiting card. 

Returning to Japan, he was appointed Professor of Physics at his old university, and for 
the next fifty years and more he flourished as an active and versatile leader of physics in 
Japan, occupying positions of increasing responsibility and influence, which in view of 
his own sentiment mentioned above need not be set out in detail here. 

Much of his early work was done on geophysics with Professor A. Tianakadate who 
was at that time Director of the Physical Institute of the university; a little later he returned 
to magnetism and with Professor K. Honda investigated the magnetostrictive properties 
of nickel-steels of various nickel contents; they showed that changes of volume, as well 
as of length, occur. ‘The geophysical work included the determination of g at Tokyo, 
comparisons of g at Tokyo and Potsdam, and a systematic study of the elasticity of rocks 
in relation to the propagation of seismic waves through the earth’s crust. His interest in 
geophysics remained with him until the end of his life, but round about 1900 he seems 
to have been strongly attracted by two widely different scientific developments : atomic 
physics arising from the discoveries of the electron, x-rays and radium on the one hand, 
and wireless communication on the other. His interest in the atom ultimately gave rise 
to a very considerable body of detailed work on spectroscopy, hyperfine structure and 
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the Zeeman effect, carried out with Dr. T. Takamine. His interest in ‘ wireless ’ seems 
to have led to his work on formulae for inductance calculation, already mentioned, and 
to a consideration of the effect of the ionosphere, familiar to him from its association with 
terrestrial magnetism which he had investigated earlier, on radio communication. 
He foresaw that the solar eclipse of 1914 would provide an opportunity of testing the 
Kennelly—Heaviside theory, if the ionized layer arises from the photoelectric effect of 
solar radiation. He made a calculation of the effect that should be observable on suitable 
wireless signals. Unfortunately the first world war prevented the making of the 
observations, but the detailed observations made in the U.S.A. at the eclipse of 1925 
showed that the idea was sound although, in the light of later information, Nagaoka’s 
original calculation was of course only the first rough approximation. 

It is estimated that Nagaoka published in all some 300 papers, and it will be realized 
from what has been said that they are so diverse in character that they do not lend 
themselves to brief effective description. What does stand out very clearly is his life-long 
devotion to physical science: we see him as mathematical physicist, an experimenter—he 
devised a magnetograph, a silica pendulum for gravity determinations and spectroscopic 
devices—a patient observer—he himself is said to have made fifty determinations in the 
gravity survey that he organized—and an original thinker. On this last point it must 
be noted that he had the honour of being the first to suggest a planetary structure for the 
atom. Rutherford in Radioactive Substances and their Radiations wrote ‘‘ This type of 
“Saturnian’ atom was first suggested in 1904 by Nagaoka, who investigated mathematically 
the stability of an atom composed of a central attracting mass surrounded by rings of 
rotating electrons. He showed that such a system was stable if the attractive force was 
large and the electrons in rapid motion’’. When it is remembered that it was in this 
year also that J. J. Thomson was still investigating the model consisting of electrons 
inside a positively charged sphere, the originality of the conception can be better 
appreciated. It was spectroscopy rather than radioactivity that was in Nagaoka’s mind; 
his aim was to account for the discrete spectral lines characteristic of the atom, and he 
thought that they might be associated with vibrations of the electrons in and perpendicular 
to the plane of the ring. We now know that the vital key to the problem was still missing, 
and that Bohr was to find it some years later, but Nagaoka’s contribution commands our 
respect. 

Nagaoka came to occupy an eminent position in the scientific life of Japan. He was 
elected a member of the Imperial Academy in 1905, became Director of the Institute of 
Physical and Chemical Research in Tokyo (now the Scientific and Research Institute) 
in 1926 and was among the first recipients of ‘ Bunkwa Kunsho’, a medal of great 
distinction, in 1937. He was President of the Imperial Academy in 1939 and 1948. 
He visited England on two or three occasions, was awarded an honorary degree at Cambridge, 
and read papers in London before the Physical Society. He was elected an Honorary 
Fellow of the Society in 1912. That this tribute to an outstanding fellow worker in the 
new school of physics then rapidly rising in Japan has survived the shock of two world wars 
may serve as a welcome reminder that physics still stands as a co-operative human enterprise. 

L. H. 


FRANK BENNETT YOUNG 


Ir is recorded with regret that Frank Bennett Young died on 29th January 1951. He wasa 
Fellow of the Physical Society from 1916 to the time of his death, and a founder Fellow (1920) 
and life member of the Institute of Physics. 

He was born in 1877 in Bristol where he received his early education. In 1896 he went 
to the Borough Road College, Isleworth, where he was awarded a scholarship enabling him 
to spend a year at Jena University. In 1899 he obtained the B.A. degree of London Univer- 
sity. At this time he became interested in science from a philosophical point of view and 
began to study physics, returning to Bristol to teach physics in a Higher Grade School. 
During this period he attended evening classes at University College, Bristol, before the 
foundation of the University. It was there he met Professor Chattock who had an important 
influence on his career. His enthusiasm for research in physics increased and he spent all 
his spare time (evenings, week-ends and holidays) on research under Chattock’s guidance. 
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His published papers at this time included ‘“‘ The field at the point of an electrified needle ”” 
(Phil. Mag., 1907, 13, 542) and “‘ Critical phenomena of ether ” (a study of the opalescent 
phenomena near the critical point) (Phil. Mag., 1910, 20, 793). On the foundation of the 
University of Bristol he took advantage of one of its early ordinances which permitted him 


me 


in 1910 to present these publications for the degree of B.Sc. by research. In 1915, soon after — 


the outbreak of the 1914-18 war, (Prof.) A. M. Tyndall obtained for him the job of localizing: 
shrapnel bullets in the body by x-rays at the Bristol Royal Infirmary, a problem quite beyond 
the medical man of that day. This was the first time a physicist had entered its doors to 
give service. ; 

In 1916 he joined the Admiralty, for anti-submarine research under the Board of Invention. 
and Research (W. H. Bragg and Rutherford were on the controlling committee). He worked. 
in Admiralty research establishments at Aberdour (Fife), Parkeston Quay (Harwich), 
Shandon (Gareloch) and, after the war, at Teddington (Middlesex). In the earlier days he 
was mainly concerned in research work on hydrophones and in underwater acoustics. 
Some of this work was published after the war, viz. ‘‘ On light body hydrophones and the 
directional properties of microphones ”’ (Proc. Roy. Soc. A, 1921, 100, 252) and “‘ On the 
acoustic disturbance produced by smali bodies on plane waves transmitted through water,. 
with special reference to the single-plate direction finder ’’ (Proc. Roy. Soc. A, 1921, 100, 
261). This had an important bearing on the design of directional hydrophones. At a later 
stage of the 1914-18 war he carried out important researches on the generation of electric 
currents in the sea (a) due to the presence of steel ships, and (b) due to the tidal flow of a 
conducting medium (the sea) in the earth’s field. Some of the experiments and conclusions. 
reached in the latter were published in Phil. Mag., 1920, 40, 149, in an article entitled 
“Electrical disturbances due to tides and waves’”’ by F. B. Young, H. Gerrard and W. 
Jevons. At this time he was also interested in the development of a new method of locating 
faults in submarine cables—using alternating currents and silver-chloride electrode search. 
gear. This method has since been extensively used by cable companies and by the G.P.O. 
Soon after the end of the 1914-18 war he obtained the degree of D.Sc. (Bristol 1919). 
On the removal of the Admiralty research establishment from Shandon (Gareloch) to: 
Teddington he worked on course plotters for ships and developed what is now known as the 
A.R.L. (Admiralty Research Laboratory) Course Plotter. Another line of research which 
he directed related to the preservation of the rubber insulation of cables in ships. It will 
be understood that much of his research work for the Admiralty remained unpublished on 
grounds of security. In July 1934 he became Superintendent of the Admiralty Research 
Laboratory, Teddington, a post which he filled until his retirement on reaching the age limit 
in October 1937. Whilst in this post he was awarded the O.B.E. 

In 1939 he returned to Admiralty work in Portsmouth soon after the outbreak of the 
1939-45 war, taking charge of an important section of research on the degaussing of ships 
as a protection against enemy magnetic mines. He finally retired from Admiralty service 
in 1941, taking on various unpaid wartime duties in Bristol. 


He eventually settled in Portsmouth after the war where he revived his old-time interests’ 


in music and in sketching and water colour painting, He was of a retiring disposition,very 
reliable, thorough and conscientious in all his work, and generous in giving more than their 
share of credit to others who worked with him. He will be missed by those who, like 
myself, knew him well. A. B. WOOD. 


ARNOLD SOMMERFELD 


On 26th April 1951, Arnold Sommerfeld, Honorary Fellow of the Physical Society, died 
from the consequences of a street accident at the age of eighty-two. 

He was born in 1868 at Kénigsberg, East Prussia, and attended school and university in 
his native town. In 1893 he went to Gottingen and came under the influence of Felix 
Klein. Here he had a thorough mathematical training and acquired the skill in the theory of 
differential equations, analytical functions and complex integration which he later applied 
with great success to physical problems. The direct fruit of his co-operation with Klein was 
their famous book on the spinning top (1895-1910). 

In 1895 Sommerfeld became lecturer in Géttingen, in 1897 professor at the Mining 
College in Clausthal and in 1900 he was appointed to a chair of mechanics at the Technical 
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College in Aachen. From 1906 until his death he was professor of theoretical physics in 
Minich. He declined several attractive positions, as for instance the chair of Planck in 
Berlin. 

Sommerfeld’s scientific work covers a vast field of which only some major points can 
here be mentioned. His contributions to classical physics were in hydrodynamics, by work 
on the stability of laminar motion, and in optics which he enriched by ingenious rigorous. 
solutions of some diffraction problems. But soon he directed his attention to the new ideas. 
of relativity and quantum mechanics. Combining these two subjects he developed Niels. 
Bohr’s theory of atomic spectra by taking into account the relativistic variability of mass, 
and obtained in 1915 his celebrated fine structure formula for hydrogen-like spectra. 
Sommerfeld’s fine structure constant dominates not only these phenomena but the whole of 
quantum electrodynamics. In 1917 be began to write his great book on the structure of 
atoms and spectral lines which has so immensely contributed to the understanding and 
general acceptance of quantum theory. The advent of quantum mechanics led him first to 
write a short complementary volume which was later extended to a comprehensive repre- 
sentation containing a wealth of information. "The study of wave mechanics in combination 
with Pauli’s principle led Sommerfeld to another great success, its application to the electron 
theory of metals. 

Sommerfeld was a brilliant teacher and attracted pupils from all over the world. When 
in 1938 he retired from his Chair, he began to prepare the publication of his lectures, of which 
all but one volume have appeared, and it is hoped that this last one will also be published. 
His style and power of representation is such that his books read like novels. But still 
greater was the charm of his personality which won him the devotion of his pupils, amongst 
whom were Pauli, Heisenberg, Brillouin, F. London and many other outstanding men. 

Sommerfeld was a member of many academies and learned societies, amongst them the 
Royal Society. He received the Max Planck Medal of the German Physical Society and 
many other honours. He was at one time Karl Schurz Professor at the University of 
Wisconsin and visiting professor at the Californian Institute of Technology and at several 
other American universities. MAX BORN. 


STUART ARTHUR FRANK WHITE 


STUART ARTHUR FRANK WHITE, who died on 12th January 1951, aged 80 years, was 
elected a Fellow of the Physical Society in 1895. He was for many years Professor of 
Mathematics at King’s College, London, Deen of the Faculty of Science and latterly 
Vice-Principal and Fellow of the College. He retired in 1933 after occupying bis Chair for 
thirty years. 

White’s education began at Portsmouth Grammar School from which he proceeded 
(with a scholarship or exhibition, I believe) to Wadham College, Oxford, at Michaelmas 
1887. In his student days he seems to have been as much interested in physics and 
astronomy as in mathematics, and before going to King’s College he held an appointment 
at Oxford under Professor Clifton. In 1894 he was appointed Demonstrator of Natural 
Philosophy—strangely enough, on the recommendation of the Engineering Commuttee— 
under Professor Grylls Adams, and in 1903 he succeeded W. H. Hudson in the Chair of 
mathematics. The appointment to this Chair of the demonstrator of natural philosophy, 
rather than of a mathematician vom Fach, appears to have been due to the desire of the 
Council to infuse a more vivid and practical spirit into the teaching of mathematics to the 
students of such subjects as engineering, physics, etc., whose interest lay less in mathematics. 
for its own sake than in its practical applications. For such students White was ideal. 
He was a most able teacher, lucid and always interesting. He was moreover whole-heartedly 
devoted to his teaching duties and to the progress and welfare of his students, whose affection 
and esteem he undoubtedly won. 

White was a keen golfer and skier and for a long period spent part of every Christmas. 
vacation in Switzerland. He was intensely devoted to the college he served and his wise 
counsel and guidance were greatly valued, not only by the Professorial Board of the College, 
but also by his mathematical colleagues in the university. WM. WILSON. 
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Miss JESSIE S. CHREE 


Fellows of the Physical Society will be sorry to learn of the death of Miss Jessie S. 
Chree, sister of the late Dr. Charles Chree, F.R.S. 

Dr. Chree was at one time President of the Physical Society, and the Society’s 
Charles Chree Medal and Prize was instituted in memory of him through the generosity 
of Miss Jessie Chree. 


ALLAN FERGUSON 


We record with regret the death on 9th November 1951 of Dr. Allan Ferguson who was 
the Papers Secretary of the Physical Society from 1928 to 1938 and President from 1938 
to 1941. 


An appreciation will be published with next year’s Obituary Notices. 


CONTENTS FOR SECTION B 


PAGE 
Professor 5. CHapMAN (President 1949-50). ; F ; : Frontispicee 
Mr. J. H. ANDREAE and Dr. JoHN Lamp. Ultrasonic Relaxation Processes in Pure 

Liquids : 1021 
Mr. E. A. NEPPIRAS anal Dr. B. E: NOpawES Cavitation produced by Uleeasorieee 

Theoretical Conditions for the Onset of Cavitation 2 : : . 1032 
Dr. M. BLackman. Self-Magnetic Field in High Current Deka: ; : 1039 
Dr. J. R. Tittman. Space Current Changes in Thermionic Valves following Sinail 

Pulses of Current . : 3 : ‘ ; ; ‘ : . . 1046 
Mr. D.A. Bett. Ionospheric Interaction in Disturbed Conditions. i . 1053 
Mr. K. JoHNsEN. Heavy Beam Loading in Linear Electron Accelerators : . 1062 
Mr. E. K. Inari. Currents in a High-Voltage Ion Accelerator’Tube . . . 1068 
Dr. K. J. Le Courrur. ‘The Regenerative Deflector for Synchro-cyclotrons . - 10733 
Letters to the Editor : 

Dr. E.O. Hatt. The Liiders Deformation of Mild Steel. ; . 1085 

Mr. A. N. Hunter. Ultrasonic Absorption in Supercooled Liquids , 1086 

Dr. D. G. Avery. ‘The Optical Constants of Lead Sulphide and Lead Telluride 

in the Region 0-5-3 microns . 1087 

Dr. W. EHRENBERG and Mr. R. E. Stay: Dr. W. Guaser. The Refractive 

Index in Electron Optics. 5 ; ¢ ‘ ; , . 1088 
Obituary Notices : 

LESLIE JOHN CoMRIE : ; 3 , ; ; : ; . 1090 
LEONARD GEORGE GRIMMETT . ; : : : : : : . 1092 
JoHN R. LoorBourow . ; P : . : ; , : . 1093 
Epwarp MALLETT . : : ; : : : : . : » 109s 
HAntTaRO NAGAOKA : : : : ; : : : : . 1095 
FRANK BENNETT YOUNG : : ; : ; : : : . 1096 
ARNOLD SOMMERFELD. : : : : 4 ‘ . 109% 
Stuart ARTHUR FRANK WHITE : : ‘ : : : : . 1098 
Miss JESSIE S. CHREE : : . 5 , 2 : : . 1099 
ALLAN FERGUSON (brief note: Obituary Notice next year) . : ; . LOGE 
Contents for Section A , : : : : : ’ : . . 1099 
Abstracts for Section A : : : F ; : 5 : . . 1100 
Subject Index, Section B, Vol. 64 ; : : : : 2 : . 1102 
Index of Authors (with Titles), Section B, Vol. 64 . : : : : . 1108 


Index to Reviews of Books, Section B, Vol. 64 5 ; : - : . 1114 


LLG 
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Ultrasonic Relaxation Processes in Pure Liquids, by J. H. ANDREAE and 
JoHN Lams. 


SUMMARY. The assumptions made by Hall and Herzfeld in previous theories of 
ultrasonic absorption are interpreted on the basis of a new and more general analysis. 
‘The formulae given by them are discussed, and it is suggested that the conditions implicit 
in Herzfeld’s formula should be modified. By this modification a new relation is derived, 
which provides a satisfactory explanation of the observed relaxation effects in those liquids 
which have been studied. 

A method of ‘curve fitting’ is employed to compare experimental results with 
theoretical requirements. Finally, it is concluded that the ultrasonic behaviour of acetic 
acid and propionic acid is governed by the excitation of bending vibrations of carbon— 
hydrogen bonds. Measurements in carbon disulphide indicate that in this case also a 
bending vibration is responsible. Some incomplete data for toluene suggest that here 
a stretching vibration of the carbon—hydrogen bonds in the methyl group is involved. 


Cavitation Produced by Ultrasonics: Theoretical Conditions for the Onset of 
Cavitation by E. A. Nepprras and B. E. NoiTIncK. 


ABSTRACT. An attempt is made to define theoretically the conditions for the 
appearance of cavitation in liquids subjected to alternating pressure changes. It is found 
that cavitation is restricted to a definite range of variations of the following parameters : 
(i) alternating pressure amplitude, (ii) frequency of the pressure wave, (iii) radius of the 
bubble nucleus, (iv) hydrostatic pressure. Under certain conditions, the change from 
non-cavitating to cavitating conditions is found to be exceedingly sharp, and in these 
‘cases it is shown that the threshold for cavitation can be accurately expressed in very 
simple terms. 


Self-Magnetic Field in High Current Discharges, by M. BLACKMAN. 


ABSTRACT. The theory of the self-magnetic field in high current discharges, originally 
developed by Tonks, has been investigated omitting a special assumption used in the work 
of Tonks. The influence of the parameters describing the discharge conditions has been 
considered in some detail. A comparison is also made with a simple, though more 
approximate, version of the theory due to Thonemann and Cowhig. ‘The two theories are 
shown to be in good agreement. 


Space Current Changes in Thermionic Valves following Small Pulses of Current, 
by J. R. TILLMan. 


ABSTRACT. A train of regularly spaced rectangular pulses of voltage, amplitude 

modulated at a frequency f,, and applied to the control grid of a commercial valve working 
‘conventionally, caused space current changes of up to 0-01 amp/cm® from the oxide-coated 
cathode. An unmodulated pulse train, interlacing the first with a delay of tg, sampled the 
anode current; the samples contained a component at the frequency f, of low level. The 
component decayed as tg increased, with a time constant of about 1 msec. It was sometimes 
in phase with the component obtained when ta=0, and sometimes in antiphase. The 
magnitude and phase were often dependent on the history of the valve. ‘Thus an operation 
known to poison the cathode was usually followed, temporarily, by a large antiphase 
component. 

An in-phase component denotes an enhancement of emission following a positive 
pulse; an antiphase component denotes fatigue. There is evidence that the impedance 
‘at the interface between cathode core and coating is not responsible for the effects. The 
fatigue may well be largely a surface, rather than a bulk or matrix, effect. It bears some 
relationship to effects reported elsewhere with pulses involving current densities of 
10 amp/cm?. 
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Ionospheric Interaction in Disturbed Conditions, by D. A. BELL. 


ABSTRACT. Measurements of ionospheric interaction depend on special radio 
transmissions which are only available on a limited number of occasions, and only on 
some of these occasions are ionospheric conditions steady enough for direct comparison. 
between experiment and theory. This paper describes the use of averaging methods to- 
extract some information from observations made when conditions are not steady. The 
results confirm the existence of a ‘low-frequency anomaly’ observed by Huxley, but 
show no exceptional interaction for a variation of the frequency of the disturbing wave 
around the nominal gyromagnetic resonance frequency. 


Heavy Beam Loading in Linear Electron Accelerators, by K. JOHNSEN. 


ABSTRACT. The influence of the beam current on the length, the obtainable particle 
energy and the efficiency of a linear accelerator is considered. The beam loading is especi- 
ally important if the accelerator waveguide has low losses. The efficiency, defined as the ratio: 
of H.F. power transformed into kinetic energy to the total H F. power fed to the accelerator,. 
can have a high value, 70-80%, but the particle energy obtainable drops considerably with. 
heavy beam loading. 


Currents in a High-Voltage Ion Accelerator Tube, by E. K. INALL. 


ABSTRACT. It is well known that a large electron current passes through a vacuum tube 
when used for the electrostatic acceleration of positive ions. This paper deals with experi-- 
ments which show that many of the electrons are produced by neutral particles bombarding 
the lower end of the tube. The particles are formed in the region of the ion source canal 
and have a low energy due to the extraction potential. It has been found possible to suppress. 
all the electrons produced below the bottom of the accelerator tube. 


The Regenerative Deflector for Synchro-cyclotrons, by K. J. LE CourEvr. 


ABSTRACT. ‘Tuck and Teng have proposed to extract the beam from a synchro- 
cyclotron by modifying the normal magnetic field in such a way as to build up the 
amplitude of radial oscillations until protons can escape from the magnet. In the present: 
paper the method is studied analytically and formulae are derived which express the 
performance in terms of the magnetic field perturbations : these formulae have been. 
evaluated numerically in sufficient detail to guide practical design. 
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Figure 3. Electrolytically etched (100) silver Figure 6. Specimen as for Figure (3); 
surface. [011] azimuth. +4° from [011] azimuth. 


Figure 8. Electrolytically etched (110) silver Figure 9. Electrolytically etched (flat) (111) 
surface. [110] azimuth. silver surface. [110] azimuth. 


a, 


Figure 11. Electrolytically polished (110) Figure 10. Electrolytically etched (rough) 
silver surface. [110] azimuth. (111) silver surface. [211] azimuth, 
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